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Unsteady Aerodynamic Models for Turbomachinery 
Aeroelastic and Aeroacoustic Applications 


Summary 


Theoretical analyses and computer codes have been developed for predicting compressible 
unsteady inviscid and viscous flows through blade rows of axial-flow turbomachines. Such 
analyses are needed to determine the impact of unsteady flow phenomena on the structural 
durability and noise generation characteristics of the blading. The emphasis here has been 
placed on developing analyses based on asymptotic representations of unsteady flow phenom- 
ena. Thus, high Reynolds number flows driven by small amplitude unsteady excitations in 
which viscous effects are concentrated within thin layers have been considered. The resulting 
analyses should apply in many practical situations and lead to a better understanding of the 
relevant flow physics. In addition, they will be efficient computationally, and therefore, appro- 
priate for use in aeroelastic and aeroacoustic design studies. Finally, the asymptotic analyses 
will be useful for calibrating and validating the time-accurate, nonlinear, Euler /Navier- Stokes 
analyses that are currently being developed for predicting unsteady flows through turboma- 
chinery blade rows. 

Under the present research program, the effort has been focused on formulating invis- 
cid/viscid interaction and linearized inviscid unsteady flow models, and on providing inviscid 
and viscid prediction capabilities for subsonic steady and unsteady cascade flows. In this 
report we describe the linearized, inviscid, unsteady aerodynamic, analysis, LINFLO, the 
steady, strong, inviscid/ viscid interaction analysis, SFLOW-IVI, and the unsteady viscous 
layer analysis, UNSVIS, that have been developed under this program. The LINFLO analysis 
can be applied to efficiently predict the unsteady aerodynamic response of a two-dimensional 
blade row to prescribed structural, (i.e., blade) motions and external aerodynamic (acous- 
tic, vortical and entropic) disturbances. The SFLOW-IVI analysis can be applied to predict 
steady cascade flows, at high Reynolds numbers, in which regions of strong inviscid/viscid 
interaction, including viscous layer separation, occur. The UNSVIS analysis can be applied 
to predict nonlinear unsteady flows in thin boundary layers and wakes. The capabilities of 
the three analyses are demonstrated via applications to unsteady flows through compressor 
and turbine cascades. The numerical results pertain to unsteady flows excited by prescribed 
vortical disturbances at inlet, acoustic disturbances at inlet and exit and blade bending and 
torsional vibrations. Recommendations are also given for the future research needed for ex- 
tending and improving the foregoing asymptotic analyses, and to meet the goal of providing 
an efficient inviscid/viscid interaction capability for subsonic and transonic unsteady cascade 
flows. 
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1. Introduction 


The unsteady aerodynamic analyses intended for use in predicting the aeroelastic and 
aeroacoustic responses of turbomachinery blading must be applicable over wide ranges of 
blade-row geometries, mean operating conditions, and modes and frequencies of unsteady 
excitation. In particular, these analyses must be capable of predicting unsteady pressure 
responses of blade rows to a variety of unsteady excitations. The latter include structural 
(blade) motions, variations in total temperature and total pressure (entropy and vorticity 
waves) at inlet and variations in static pressure (acoustic waves) at inlet and exit. Finally, 
because of the large number of controlling parameters involved, there is a stringent requirement 
for computational efficiency, if an analysis is to be used successfully in the blade design process. 
To sati f y this latter requirement a nu. iber of restrictive assumptions must be introduced into 
the development of appropriate unsteady aerodynamic models. The analyses described in this 
report have been developed to provide reliable and efficient theoretical prediction capabilities 
for inviscid and viscid, steady and small-disturbance unsteady, flows, at high Reynolds number, 
through two-dimensional cascades. 

1.1 Background 

The theoretical analyses that have been developed to predict the aeroelastic and aeroa- 
coustic behavior of turbomachinery blading, i.e., the onset of blade flutter, the amplitudes 
of forced blade vibration and the sound pressure levels that exist upstream and downstream 
of the blade row have, for the most part, been based on the following geometric and aerody- 
namic assumptions. The blades of an isolab 1, two-dimensional cascade are usually considered 
with the aerodynamic effects associated with neighboring structures being represented via pre- 
scribed nonuniform flow conditions at inlet and exit. The unsteady excitations are assumed to 
be periodic in time and in the blade-to-blade direction. The flow Reynolds number is taken to 
be sufficiently high so that viscous effects have a negligible impact on the unsteady pressure 
field. Finally, the unsteady excitations axe assumed to be sufficiently small so that a linearized 
treatment of the unsteady inviscid flow is justified. 

Until fairly recently, the inviscid unsteady aerodynamic analyses that have been available 
for turbomachinery aeroelastic and aeroacoustic design applications have been based on clas- 
sical linearized theory (see [Whi87] for a review). Here the steady and harmonic unsteady 
departures of the flow variables from their uniform free-stream values are regarded as small 
and of the same order of magnitude, leading to uncoupled, linear, constant coefficient, bound- 
ary value problems for the steady and the complex amplitudes of the unsteady disturbances. 
Thus, unsteady solutions based on the classical linearization apply essentially to cascades 
of unloaded flat-plate blades. Very efficient semi-analytic solution procedures have been de- 
veloped for two-dimensional attached subsonic and supersonic flows and applied with some 
success in turbomachinery aeroelastic and aeroacoustic design calculations. It should also be 
mentioned that extensive efforts (as reviewed by Namba [Nam87]) have been made to develop 
three-dimensional unsteady aerodynamic analyses, based on the classical linearization. 

Because of the limitations in physical modeling associated with the classical linearization, 
more general two-dimensional inviscid linearizations have been developed. These include the 
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effects of important design features such as real blade geometry, mean blade loading, and 
operation at transonic Mach numbers. Here, unsteady disturbances are regarded as small- 
amplitude harmonic fluctuations relative to a fully nonuniform, isentropic and irrotational, 
mean or steady background flow. The steady flow is determined as a solution of a nonlinear 
inviscid equation set, and the unsteady flow is governed by linear equations with variable 
coefficients that depend on the underlying steady flow. This type of analytical model is de- 
scribed in [Ver92, Ver93] and is often referred to as a linearized potential model; however, 
the potential-flow restriction applies only to the steady background flow. It has received con- 
siderable attention in recent years, and solution algorithms for the nonlinear steady and the 
linearized unsteady problem have reached the stage where they are being applied in turbo- 
machinery aeroelastic and aeroacoustic design studies (e.g., see [Smi90, SS90, MM94]). In 
particular, one such analysis, the linearized inviscid flow analysis, LINFLO, has been ex- 
tended under the present contract for cascade gust response predictions. In previous work, 
this analysis has been developed and applied to predict unsteady subsonic and transonic flows 
excited by blade vibrations or acoustic disturbances [VC84, Ver89a, UV91, KV93b] and, un- 
der the present contract, to predict unsteady subsonic flows excited by entropic and vortical 
gusts [VH90, HV91, VBHA91]. 

The unsteady flows of practical interest usually occur at high, but finite Reynolds number, 
so that viscous-layer displacements can have an impact on the unsteady pressure response. 
Provided that large-scale flow separations from the blade surfaces do not occur, the overall 
flow field can be divided conceptually into “inner” viscous or dissipative regions, consisting of 
thin boundary layers and wakes, and an “outer” inviscid region. Solutions to the complete flow 
problem can then be determined by an iterative procedure involving successive solutions to the 
inviscid and viscid equations. If the inviscid/ viscid interaction is “weak”, then at each step of 
the iteration process, the inviscid and viscid solutions can be determined sequentially with the 
pressure being determined by the inviscid flow. However, in most flows, strong inviscid/ viscid 
interactions occur due, for example, to boundary-layer separations, shock/boundary-layer 
interactions and trailing-edge/ near- wake interactions. For such flows the pressure must be 
determined by solving the inviscid and viscous layer equations simultaneously at each iteration 
level. 

The construction of a high Reynolds number viscous cascade solver involves first, the de- 
velopment of component flow solvers, and second, the implementation of these component 
solvers into an overall computational procedure to provide an analysis for the complete flow 
field. Solution methods for steady, subsonic and transonic, inviscid flows through cascades and 
for steady boundary-layer and wake flows have been developed to a relatively mature state. 
Methods for coupling such solutions have also been developed and applied to predict steady 
cascade flows, with strong inviscid/viscid interactions (e.g., see [HSS79, JH83, CH80, BV89, 
BHE91, BVA93a, BVA93b]), with the work in [BVA93a, BVA93b] being performed as part 
of the present effort. Similar approaches are needed for unsteady flows. Under the present 
Contract, nonlinear steady and linearized unsteady inviscid analyses have been coupled to the 
unsteady viscous-layer analysis of [PVK91], but only to predict unsteady cascade flows with 
weak inviscid/viscid interactions [VBHA91, BV93]. To date, there has been very little effort to 
couple unsteady inviscid and viscous-layer analyses to provide a strong inviscid/viscid interac- 
tion analysis for unsteady cascade flows. As steps toward this goal, the steady inviscid/viscid 
interaction analysis for cascades, SFLOW-IVI, that can provide the foundation for an un- 
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steady procedure to be developed later, and the unsteady viscous-layer analysis (UNSVIS) 
have been developed under the present research program. 

1.2 Scope of the Present Effort 

The objectives of the research program conducted under Contract NAS3-25425 are to 
provide efficient theoretical analyses for predicting compressible unsteady flows through two- 
dimensional blade rows. Such analyses are needed to understand the impact of unsteady 
aerodynamic phenomena on the aeroelastic and aeroacoustic performances of turbomachinery 
blading. For this purpose, we have developed a detailed inviscid/ viscid interaction formu- 
lation for unsteady cascade flows, which is described in § 2 of this report. We have also 
developed or extended several analyses that will form the components of an unsteady, invis- 
cid/viscid interaction, solution procedure. The latter include the linearized inviscid unsteady 
aerodynamic analysis, LINFLO, the steady inviscid/ viscid interaction analysis, SFLOW-IVI, 
and the unsteady viscous layer analysis, UNSVIS. The steady full potential analysis, SFLOW 
[HV93, HV94], which is also considered to be a component of the overall prediction scheme, 
has been developed at NASA Lewis under a separate, but related, research program. The 
present work has been directed primarily towards subsonic aeroelastic applications, however, 
for the most part, this work will apply, more generally, for predicting the aeroelastic and 
aeroacoustic responses of turbomachinery blading operating at subsonic and tr ^ sonic Mach 
numbers. 

In the first phase of this program [VH90, HV91] the linearized inviscid analyt LINFLO) 
was extended to predict the responses of a cascade to entropic and vortical excitations. A 
velocity decomposition introduced by Goldstein [Gol78, Gol79], and later modified by Atassi 
and Grzedzinski [AG89], is employed to split the linearized unsteady velocity into rotational 
and irrotational components. This decomposition leads to a very convenient descrip on of 
the linearized unsteady perturbation — one in which closed form solutions can be determined 
for the entropy and vorticity or rotational velocity fluctuations in terms of the drift and 
stream functions of the underlying steady flow. Numerical field methods are required only 
to determine the unsteady potential and from this, the unsteady pressure. The potential is 
governed by an inhomogeneous wave equation in which the source term depends upon the 
rotational velocity field. The potential fluctuations are determined numerically on an H - type 
mesh in which the streamlines of the steady background flow are used as mesh lines. Numerical 
solutions are reported in [VH90, HV91, VBHA91] for several configurations including flat- 
plate cascades, a compressor exit guide vane, a high-speed compressor cascade, and a turbine 
cascade. The LINFLO analysis is described in § 3 of this report, with particular emphasis on 
its capabilities for predicting unsteady flows excited by vortical gusts. 

An efficient steady analysis for predicting strong inviscid/ viscid interaction phenomena; 
such as, viscous-layer separation, shock/boundary-layer interaction and trailing-edge/near- 
wake interaction, in turbomachinery blade passages is needed as paxt of a comprehensive 
analytical blade design prediction system. Such an analysis, called SFLOW-IVI, has been 
reported in [BVA93a, BVA93b] and is described in § 4. Here, the flow in the outer or in- 
viscid region is governed by the full-potential equation and that in the inner viscous region, 
by Prandtl’s viscous-layer equations. The non-hierarchical nature of strong interactions is 
taken into account in the semi-inverse iteration procedure used to couple the two solutions. 
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The steady, full-potential analysis, SFLOW, [HV93, HV94] is employed in this procedure 
to determine inviscid solutions. SFLOW was constructed, for use with the LINFLO analy- 
sis, to provide comprehensive and compatible, steady and unsteady, inviscid, flow prediction 
capabilities for cascades. 

In the IVI procedure, which is referred to as SFLOW-IVI, viscous effects are incorpo- 
rated by adjusting the blade and wake surface boundary conditions in SFLOW to account 
for the effects of viscous displacement. Inviscid solutions are determined using an implicit, 
least-squares, finite-difference approximation, viscous-layer solutions using an inverse, finite- 
difference, space-marching method which is applied along the blade surfaces and wake stream- 
lines. The inviscid and viscid solutions are coupled using a semi-inverse, global iteration proce- 
dure, which permits the prediction of boundary-layer separation and other strong-interaction 
phenomena, and allows nonlinear changes to the inviscid flow, due to viscous effects, to be 
evaluated. Results are presented for three cascades, with a range of inlet flow conditions 
considered for one of them, including conditions leading to large-scale flow separation. Com- 
parisons with Navier-Stokes solutions and experimental data are also given. 

Finally the UNSVIS analysis of [PVK91] was extended so that unsteady viscous effects 
in the vicinity of leading-edge stagnation points and in blade wakes could be predicted 
[VBHA91, BV93]. The nonlinear unsteady flow in a viscous layer is described by Prandtl’s 
equations. As in the SFLOW-IVI analysis, algebraic models are used to account for the effects 
of transition and turbulence. The viscous-layer equations are solved in terms of Levy-Lees 
type variables using a finite-difference technique in which solutions are advanced in time and 
in the streamwise direction. Numerical solutions are determined by marching implicitly, first 
in time and then in the streamwise direction, over several periods of unsteady excitation, 
from an initial steady solution, and from an approximate time-dependent, upstream flow so- 
lution. This analysis, rather than one in which the results of separate nonlinear steady and 
linearized unsteady viscous-layer analyses are superposed, allows an assessment to be made 
of the relative importance of nonlinear unsteady effects in viscous regions. 

Under the present effort, a similarity analysis was developed to predict unsteady viscous 
compressible flow in the vicinity of a moving leading-edge stagnation point and incorporated 
into the UNSVIS code. The stagnation region analysis provides the “initial” upstream flow 
information needed to advance or march a viscous-layer calculation downstream along the 
blade surfaces and into the wake. In addition, the wake analysis, used previously in UNSVIS, 
was extended so that the changes or jumps in the inviscid velocity that occur across vortex- 
sheet unsteady wakes could be properly accommodated. The linearized inviscid analysis, 
LINFLO, and the nonlinear viscous-layer analysis, UNSVIS, were also coupled to provide a 
weak viscid/inviscid interaction solution capability for unsteady cascade flows. The UNSVIS 
analysis is described in § 5, where it is also applied to study the viscous-layer responses of 
an unstaggered flat-plate cascade to pressure or acoustic excitations originating upstream and 
downstream of the blade row. Finally, the coupled LINFLO/UNSVIS analysis is applied to a 
turbine cascade subjected to a pressure excitation from upstream to demonstrate the current 
weak inviscid /viscid interaction solution capability on a realistic cascade configuration. 

The component analyses, described in this report, are important in their own right. They 
can be applied to improve our understanding of the complex steady and unsteady flow pro- 
cesses that occur in turbomachine cascades, and to provide useful aeroelastic and aeroacoustic 
design information. Hopefully, along with the general inviscid/viscid interaction and inviscid 
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small-disturbance formulations, presented in § 2, they will contribute to the future develop- 
ment of a comprehensive inviscid/ viscid interaction analysis for unsteady cascade flows. Such 
an analysis will provide useful aeroelastic and aeroacoutic response information for a wide 
range of blade-row geometries and operating conditions. It will also be of help in calibrat- 
ing the time-accurate, nonlinear, Euler and Navier-Stokes analyses that are currently being 
developed for predicting turbomachinery unsteady flow fields. 
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2. Physical Problem and Mathematical Models 


2.1 Unsteady Flow through a Two-Dimensional Cascade 

We consider time-dependent flow, at high Reynolds number (Re) and with negligible body 
forces, of a perfect gas with constant specific heats and constant Prandtl number (Pr) through 
a two-dimensional cascade, such as the one shown in Figure 2 . 1 . The unsteady fluctuations 
in the flow arise from one or more of the following sources: blade motions, upstream total 
temperature and total pressure disturbances, and upstream and/or downstream static pres- 
sure disturbances that carry energy toward the blade row. These excitations are assumed to 
be of small amplitude, periodic in time, and periodic in the blade-to-blade or cascade “cir- 
cumferential” direction. Also, blade shape and orientation relative to the inlet freestream 
direction, the inlet to exit mean static pressure ratio and the amplitudes, modes, frequencies 
and wave numbers of the unsteady excitations are such that viscous effects are confined within 
thin layers, that lie along the blade surfaces and extend downstream from the blade trailing 
edges. We should note that the aerodynamic models to be developed below, apply to subsonic 
flows and to transonic flows containing normal shocks; however, the applications presented 
throughout this report are restricted to subsonic flows. 

In the following discussion all physical variables are dimensionless. Vector quantities are in 
boldface type, and a tilde over a dependent variable indicates time dependence. Lengths have 
been scaled with respect to blade chord (/>*), time with respect to the ratio of blade chord 
to upstream freestream flow speed (Vj^J, density with respect to the upstream freestream 
density (pl^), velocity with respect to the upstream freestream flow speed, and stress with 
respect to the product of the upstream freestream density and the square of the upstream 
freestream speed (pLoo^-oo 2 )- Here the superscript * denotes a dimensional quantity and the 
subscript — oo refers to the prescribed freestream conditions far upstream. The scalings for 
the remaining variables can be determined from the equations that follow, which have the 
same forms as their dimensional counterparts. 

We will analyze the unsteady flow in a blade-row fixed coordinate frame in terms of the 
Cartesian ( x,y ) or (£,77) coordinates and the time t. Here, for example, £ and 77 measure 
distances in the cascade axial and circumferential directions, respectively. To describe flows 
in which the fluid domain varies with time it is useful to consider two sets of independent 
variables, say (x,<) and (x,t). The position vector x(x,t) = x + 72.(x,t) describes the instan- 
taneous location of a moving point, say V, x refers to the reference or steady-state position 
of V, and 7?.(x, t) is the displacement of V from its reference position. The displacement 
field, n, is usually prescribed so that the solution domain moves with solid boundaries and 
is stationary far from the blade row. If we set 72. = 0, then the position vector x, like x, 
describes a stationary point in the blade-row fixed reference frame. 

The mean or steady-state positions of the blade chord lines coincide with the line segments 
77 = f tan 0 + mG, 0 < f < cos 0 , m = 0, ±1, ± 2 , ... , where m is a blade number index, 0 is 
the cascade stagger angle, and G is the cascade gap vector which is directed along the 77-axis 
with magnitude equal to the blade spacing. The blade motions are prescribed as functions of 
x and t, i.e., 

72. Bm (x + 772 G,t) = Re{RB(x)exp[7(wt -I- ttict)]}, for x€B (or x € B) . (2.1) 
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Figure 2.1: Two-dimensional compressor cascade. 

Here 7Z.B m is the displacement of a point on a moving blade surface (B m ) relative to its mean 
or steady-state position (B TO ); Rb is the complex amplitude of the blade displacement; a is 
the phase angle between the motions of adjacent blades; Re{ } denotes the real part of { }; 
and B denotes the reference (m = 0) blade surface. In the present report we will use the 
notations 7?-B m or x 6 B m , to indicate the displacement of a point on the mth blade 
surface. Similar notations, i.e., 'R. Wm or TZ, x € W m , and 'R-sh m ,„ or x € Sh m>n , will be 
used to represent inviscid wake and shock displacements. 

For aeroelastic and aeroacoustic applications we axe usually interested in a restricted class 
of unsteady flows; those in which the unsteady fluctuations can be regarded as perturbations 
of a background flow that is steady in the blade-row frame of reference. Thus, we consider 
situations in which the background flows far upstream (say £ < £_) and far downstream (£ > 
£+) from the blade-row consist of at most a small steady perturbation from a uniform flow. In 
this case any arbitrary unsteady aerodynamic excitation of small amplitude can be represented 
approximately as the sum of independent entropic, vortical, and acoustic disturbances that 
travel towards the blade row, as indicated in Figure 2.2. 

The entropic, 5_oo(x, t), vortical, ^'_ 00 (x, t), and acoustic, p/, T<x >(x, f), excitations, where 
the subscripts — oo and +oo refer to the far upstream and fax downstream flow regions, 
respectively, are also prescribed functions of x and t. However, these functions must be 
solutions of the fluid-dynamic field equations that describe disturbances that travel towards 
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S_oo(x,f) 



C-oo(*,<) 


Pl,+ oo(x, t ) 




Figure 2.2: Unsteady excitations: blade motion, incident vortical and entropic disturbances 
from upstream, and incident acoustic disturbances from upstream and downstream. 


the blade row. Since the unsteady aerodynamic excitations are periodic in fj and t, it is often 
useful to represent them via Fourier series. In this case the complex amplitudes, temporal 
frequencies and circumferential wave numbers of the various excitations are the prescribed 
quantities. In the present study, the external aerodynamic excitations are small-amplitude, 
harmonic perturbations of a uniform freestream flow, and therefore must be of the form 

s(x,t) = S_oo(x- V-oot) = Rets.*, exp[i(rc_oo -x + ud)]} ... , £ < £_ , (2.2) 

C(x,t) = C-oo(x-V_ 00 t) = Re{C_ 00 exp[*(K_ 00 -x + u;t)]}... , £ < £- , (2.3) 

and 

p/(x, t) = p/ )T oo(x, t) = R e{p/, Too expf-/?:,:^ + i(/c Too • x + u;f)]} . . . , £ > • (2.4) 

Here s_oo, C-oo and Pi , too are the complex amplitudes of the entropic, vortical and acoustic 
excitations, u is the temporal frequency, k Too is the wave number, with k^oo = crG~ 
and the fi Too are attenuation constants. It follows from the governing equations that the 
temporal frequency and wave number of an entropic or vortical excitation are related by 
a? = -k_oo • V_oo, where V_oo is the uniform relative inlet velocity, but more complicated 
relationships exist between u ; and k t<X) for pressure excitations [Ver89b]. 
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2.2 Reynolds-Averaged Navier-Stokes Equations 


The field equations that govern the fluid motion follow from the mass, momentum and 
energy conservation laws. After taking ensemble averages of the resulting equations, the 
so-called Reynolds-averaged Navier-Stokes equations for the (statistical) mean values of the 
flow variables can be determined. The conservation law forms of these equations axe usually 
preferred in modern numerical simulations of steady and unsteady flows. However, for the 
present purpose of developing computationally efficient aerodynamic models, we will consider 
the conservative form of the mass equation and convective forms of the momentum and energy 
equations, i.e., 

+ V x • {pV) = 0 (2.5) 

X 

_ D V - 

+ V X P = V x • (17 - pv' ® v') (2.6) 


dp 

dt 


and 


_DH t _ dP 
P Dt dt 


— V x • (IT • V — Q — ph' T \') 


(2.7) 


Here, D/Dt = d/dt | x + V ■ V x is the material or convective derivative operator; p, V, P 
and Ht axe the ensemble or Reynold’s-averaged values of the time-dependent fluid density, 
velocity, pressure, and specific tot al ent halpy, respectively; J7 is the viscous stress tensor; Q 
is the heat flux vector; v' <g> v' and h' T \' axe statistical correlations that account for the effects 
of random turbulent fluctuations on the ensemble- averaged flow; and ® denotes the tensor or 
dyadic product of two vectors. To obtain a closure of the mean flow equations, the turbulent 
correlations can be related to gradients in the ensemble-averaged vaxiables via algebraic eddy 
viscosity models. Note that, as is the usual practice in aerodynamic calculations, we have 
only included selected turbulent correlations in equations (2.6) and (2.7). For example, we 
could have also included density-velocity correlations in the mass, momentum and energy 
equations. We are assuming, therefore, that the correlations omitted have a negligible impact 
on the deterministic mean flows under consideration. 

For a Newtonian fluid with zero coefficient of bulk viscosity, we can write 


n = -piRe)- 1 [(2V X • V/3)I -V X ®V-(V X ® V c ] . (2.8) 


Here, p is the coefficient of shear viscosity, Re = (Tl^VZ^L* / is the flow Reynolds number, 
I is the unit tensor, and the subscript C denotes the conjugate dyadic. Assuming Fourier’s 
Law for the conduction of heat, the heat flux Q is related to the gradient of the temperature 
f by 

Q = -j kPr-'Re-'V^f , (2.9) 

where k is the coefficient of thermal conductivity, Pr = p m Cp/k* is the flow Prandtl number 
and Cp is the dimensional value of the specific heat at constant pressure. The coefficients p. and 
Jc in the foregoing equations are usually assumed to be known functions of the temperature. 
For example, the following empirical laws relating the molecular viscosity and the thermal 
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conductivity to the temperature, 


A 



t \ 3/2 r-oo + r c 

T. oo) T + Tc ’ 


( 2 . 10 ) 


are often used in Navier-Stokes calculations. Equation (2.10) is a form of Sutherland’s equa- 
tion. Here T_ «, = CpTl^/iV*^) 2 is the non-dimensional reference temperature, and T c is a 
constant, which for air has a dimensional value, T£, of 110°K [Sch60). 

The pressure, P, and temperature, T, can be expressed in terms of the dependent variables 
p, V and Hr- In particular, for a thermally and calorically perfect gas 

P = 7 _1 (7 “ 1 )p(Ht ~ V 2 /2) = 7' 1 (7 - 1 )~PH ( 2 . 11 ) 


and 

T = (H t ~ V 2 /2) = H , (2.12) 

where 7 is the specific heat ratio (constant pressure to constant volume) and H is the specific 
enthalpy of the fluid. Although enthalpy and temperature are different fluid dynamic proper- 
ties, the relationship H = T applies here because of the scalings used to non-dimensionalize 
the various flow variables. 

We will also find it useful to introduce the fundamental thermodynamic identity 

TdS = dH — p~ l dP, (2.13) 

where S is the specific entropy of the fluid. It then follows from (2.11) through (2.13) that 

PTe-^ = . ( 2 - 14 ) 

where the subscript Ref refers to a reference thermodynamic state, which is usually taken to 
be the free-stream state far upstream of the blade row, i.e., at £ = £_<»• We can make use of 
(2.13) and the equation of motion (2.6) to re-write the energy equation (2.7) in the form 


~pT^- = T - V x • (Q + P~h' T V) + (V • V x )(^v' ® v') , (2.15) 

ly L 

where T = tl : V x <g> V is the viscous dissipation and the symbol : indicates the scalar 
product of two tensors. 


Boundary Conditions 

For application to turbomachinery unsteady flows the foregoing field equations must be 
supplemented by conditions at the vibrating blade surfaces and conditions at the inflow and 
outflow boundaries. Since transient unsteady aerodynamic behaviors are usually not consid- 
ered, a precise knowledge of the initial state of the fluid is not required. The no-slip condition, 

\ = 71 for x € B m (or x € B m ) , (2.16) 
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where 'R-B m is prescribed, applies at blade surfaces. In addition, either the heat flux Q • n or 
the temperature T must be prescribed at such surfaces, i.e., 

Q-n = Q w (x,t) or T = T w (x,t), x € B m (or x € B m ) (2.17) 

We also require conditions on the flow far upstream and far downstream from the blade 
row, i.e., at the inflow and outflow boundaries of the computational domain. Typically the 
circumferentially- and temporally-averaged values of the total pressure, total temperature 
and the inlet flow angle axe specified at the inflow boundary. At the outflow boundary, the 
“circumferentially” and temporally averaged pressure is specified. In addition, total pressure 
and total temperature fluctuations at inlet and pressure fluctuations at inlet and exit, that 
carry energy towards the blade row, must be specified. Total pressure and total temperature 
fluctuations at exit and unsteady pressure disturbances at inlet and exit, that carry energy 
away from the blade row, are determined as part of the unsteady solution. 

2.3 Inviscid/ Viscid Interaction Model 

The Reynolds-averaged Navier-Stokes equations must be considered if viscous effects are 
expected to be important throughout the fluid domain. However, for most flows of practical 
interest the Reynolds number (Re) is usually sufficiently high so that such effects are concen- 
trated in relatively thin layers across which the flow properties vary rapidly. Provided that 
large scale flow separations do not occur, these layers generally lie adjacent to the blade sur- 
faces (boundary layers) and extend downstream from the blade trailing edges (wakes). Thus, 
if we assume that the Reynolds number is high and the flow remains essentially “attached” to 
the blade surfaces, we can apply an inviscid/viscid interaction (IVI) analysis to determine the 
unsteady flow field. The terminology “inviscid/viscid interaction” refers to all flow situations 
in which viscous layers have a significant influence on the pressure field. Weak interactions 
are defined as those in which viscous effects on the pressure are small, or more specifically, on 
the order of the viscous-layer displacement thickness, 8 ~ 0(Re~ 1 ^ 2 ). If viscous effects on the 
pressure disturbance are larger than this, i.e., > 0(8), the interaction is classified as strong 
[Mel80]. 

In an IVI analysis separate sets of approximate equations, i.e., reduced forms of the Navier- 
Stokes equations, are constructed, using the method of matched asymptotic expansions, to 
describe the flows in “outer” inviscid and “inner” viscous-layer regions. This approach offers 
the potential for providing efficient predictions of the effects of viscous layer displacement and 
curvature on the unsteady aerodynamic behaviors of blade rows. The governing equations can 
be derived with respect to moving, shear-layer surfaces S m , each of which is contained entirely 
within a viscous layer. These surfaces are usually taken to coincide with the suction and 
pressure surfaces of the blades and lie entirely within the viscous wake, e.g., see Figure 2.3. 
The curvatures of the reference shear layer surfaces are assumed to be 0(1) or less, and here 
we will assume that their motions, which depend on the blade motions, are also of 0(1) or 
less. Ultimately, we will assume that the motions of the blades and their wakes, and hence, 
those of the reference shear-layer surfaces are of small amplitude. The inviscid and viscous 
layer equations must be solved simultaneously, subject to appropriate matching conditions 
along the blades and wakes, to determir ' the flows in the inviscid and viscous-layer regions. 
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Figure 2.3: Two-dimensional steady flow, at high Reynolds number, over a blade surface. 

In the limit as Re — ► oo, viscous effects become negligible, and it is sufficient to consider only 
the outer or inviscid flow. 

To arrive at the governing equations for the flows in the inner and outer regions the 
dependent fluid variables are expanded in power series in the small parameter, 8. For example, 
the series expansions for the pressure have the form 

P = P* + ~6Pf + 6P* + ... and P = P° + 8P° + IP? + . . . (2.18) 

where the superscripts X and O refer to the inner and outer regions respectively. The lead 

terms in the inner-region expansions for the components of the fluid velocity, V, tangential 

and normal to the shear layer surface S, i.e., in the r— and n— directions, are assumed to 

be of 0(1) and of 0(5), respectively. The lead terms in the remaining inner region and in 

all outer region series expansions for the deterministi c uns teady flow variables are assumed to 

be of 0(1). The turbulent correlations v'<g> v' and h' T \' are assumed to be 0(8) in viscous 

~2 

regions and of higher order, i.e., 0(8 ) or higher, in the inviscid region. 

Scalings for the independent variables are also introduced. The thickness of the viscous 
layers is assumed to be 0(8), hence, in the inner region, distances, n, normal to the shear 
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- .... . . - _1 
layer surfaces are of 0(8) and partial derivatives in the normal direction are of order 6 , i.e., 

n ■ V x « 2- « 0(1 *) . (2.19) 

on 

Partial derivatives along the surface, r • V x , and local time derivatives, d/dt | x , are assumed 
to be of 0(1) in both the inner and outer regions. In the outer region, normal derivatives, 
n • V x , are also assumed to be of 0(1). 

Approximate field equations that describe the flows in the inner and outer regions are 
determined by simply substituting the series expansions for the dependent variables and the 
scalings for the independent variables into the Reynolds-averaged Navier-Stokes equations and 
equating terms of like power in 6. When expressed in terms of the original variables, P, V, 
etc., the zeroth-order inviscid flow is governed by the Euler equations, i.e., equations (2.5), 
(2.6) and (2.7) or (2.15) with right-hand-sides set equal to zero; the zeroth-order viscous flows, 
by Prandtl’s viscous-layer equations, which will be given below. The inviscid flow in the outer 
region is determined as a solution of the Euler equations subject to flow tangency conditions 
at the blade surfaces and jump conditions at shocks and at blade wakes. The blade and 
wake conditions contain terms that account for the effects of viscous-layer displacement and 
curvature, which become negligible as Re — * oo. The viscous flows in the inner regions are 
determined as solutions to Prandtl’s equations subject to edge conditions that depend on the 
behaviors of the inviscid flow variables along the blade and wake surfaces. The specific forms 
of the inviscid and viscid blade and wake conditions follow from an asymptotic matching of 
the outer inviscid and the inner viscous-layer solutions, e.g., see [Mel80]. 

Inviscid Region 


The field equations that govern continuous, inviscid, fluid motion axe determined from the 
Reynolds averaged Navier-Stokes equations, (2.5), (2.6) and, in this case, (2.15), and are given 
by 


d 2 

at 


+ V x • (V~p) = 0 

(2.20) 

X 


DV 


Dt +V * P = 0 

(2.21) 

DS 


Dt=° 

(2.22) 


For turbomachinery applications, we require solutions of these equations subject to boundary 
conditions at moving blade surfaces, jump conditions at moving wake and shock surfaces, and 
appropriate conditions far from the blade row. 

The inviscid solution for the normal component of the fluid velocity at a moving blade 
surface must match the viscous solution for this velocity at the outer edge of the viscous 
layer. This is equivalent to the condition that the inviscid flow must be tangential to the 
blade and wake displacement surfaces. After carrying out the asymptotic matching process 
and neglecting terms of second and higher order in 6, we find that the normal component of 
the inviscid fluid velocity must satisfy the condition 


(V — TZ ) . n = K 1 [dikh/dt | x + dCP'VrJydr) + (h • V x )l , x € B m (2.23) 
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at the mth (m = 0, ±1,±2, ...) moving blade surface. In equation (2.23), p e and V T>e are 
the inviscid values of the fluid density and streamwise velocity at the moving blade surface, 
B m , or the viscous density and streamwise velocity at the edge (e) of the viscous layer, r is 
the distance measured along this surface downstream from the leading-edge, and n is a unit 
vector normal to B m and pointing into the fluid. 

Two types of terms arise from the wake matching conditions, one due to the displacement 
thickness effect and the other, to the wake curvature effect. The first leads to the requirement 
that the inviscid solution for the normal component of the fluid velocity must be discontinuous 
across a wake with jump given by 

[V] • n = <?; X [d(p.3)/&| x + d(}'V r J)/dr]) + (h • V x )(l) , x € W m . (2.24) 

Here [ J and ( ) denote the difference (upper minus lower) and the sum (upper plus lower), 

respectively, across a wake, and n is an “upward” pointing, unit vector, normal to the moving 
reference wake surface, W. Note that (6) = 8 W is the displacement thickness of the complete 
wake. The wake curvature effect gives rise to a pressure difference across the wake. The 
requirement that the outer inviscid flow match this pressure difference leads to the condition 

[P] = {pAAZ + h [HVr,e -hr)- n • dh/d f] ) , X € W m (2.25) 

where (0) = 6w is the momentum thickness of the complete wake and k = r • dn/dr\ w is 
the curvature of the reference wake surface, which is taken as positive when the curvature 
is concave upwards. Since there is some ambiguity in establishing the shape and location of 
a viscous wake the right-hand-side of (2.25) is usually ignored in inviscid/viscid interaction 
calculations. Consequently, the pressure jump across the wake is usually set equal to zero. 

In deriving the surface conditions (2.23)-(2.25), we have assumed that the viscous layer is 
thin, i.e., $ « 1, and therefore, that the effects of terms of second and higher order in the 
viscous displacement thickness are negligible. In the inviscid limit Re — ► oo, the thicknesses 
of the viscous layers, and hence, the right-hand-sides of (2.23)— (2.25) become zero. 

Jump conditions must also be imposed at inviscid shock discontinuities. Such conditions 
are obtained from the integral conservation laws by considering a control volume that contains 
a segment of a shock surface, and taking limits, first, as the lateral extent of this volume, nor- 
mal to the surface segment, approaches zero, and, then, as the area of the segment approaches 
zero. The resulting jump conditions for conserving mass momentum and energy at a shock 
are given by 

[AT/] | = 0, Af/|[V] + [P]n = 0 , and + [PV] • n = 0 , x 6 Sh m<n (2.26) 

respectively. Here | J denotes the jump in a flow quantity as experienced by an observer in 
moving across the shock Sh miTl in the n direction, 

Mf = p(Y -h) n , x € Sh m<n (2.27) 

is the fluid mass flux through the shock, J E T = H T - P/p is the total specific internal energy 
of the fluid and the subscripts m,n refer to the nth shock associated with the mth blade. 
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In the inviscid limit, the conditions (2.26) also apply across the vortex-sheet “viscous” 
layers. In this case, since J V T ] ± 0, and the inviscid jump conditions (2.26) indicate that, for 
vortex sheets, AT/ = 0, [ jPJ = 0, and | Vj-n = 0. At shocks, 5h m ,„, the mass flux is generally 
nonzero (i.e., M/ ^ 0). Hence, it follows from (2.26) that the component of fluid velocity 
tangent to a shock surface, V • r, must be continuous across the shock. The remaining jump 
conditions, along with the thermodynamic equations of state, are then required to determine 
the shock velocity, 'R-sh m , n ? and the changes in the normal component of the fluid velocity and 
the thermodynamic properties of the fluid as it passes through the shock. 

The far-field conditions used in the inviscid approximation are the same as those indicated 
in § 2.2 for Navier-Stokes simulations. In particular, averaged values of the inlet total pressure, 
total temperature and flow angle and the exit static pressure are specified along with the 
entropic and vortical fluctuations at inlet and the static pressure disturbances at inlet and 
exit that carry energy towards the blade row. In addition, disturbances generated within the 
solution domain must be allowed to pass through the inflow and outflow boundaries without 
reflection. 


Viscous Layers 


The flows in the viscous layers are governed by Prandtl’s equations and are subject to 
no-slip and prescribed heat-flux or wall temperature conditions, cf. (2.16) and (2.17), at the 
moving blade surfaces. In addition, the streamwise velocity and the thermodynamic properties 
of the fluid at the edges of the viscous layers are determined by the values of the corresponding 
inviscid quantities at the blade surfaces and along the reference wake (shear-layer) surfaces. 
We describe the flow in these layers in terms of curvilinear coordinates r and n which measure 
distance along and normal to, respectively, a moving, reference, shear-layer surface S. The 
unit vectors r(f,f) and n(f , t), where f measures distance along the mean position of the 
shear layer surface, are tangent and normal to S. 

The field equations that govern the lead terms, Vj 0 , Vf • n, P£, etc., in the inner-region 
series expansions, are 'etermined in a manner similar to that used for the inviscid region. In 
particular, the series expansions for the inner-region dependent variables and the scalings for 
the independent variables are substituted into the Reynolds-averaged Navier-Stokes equations, 
(2.5), (2.6) and (2.7), and only terms of 0(1) are retained. Then, in terms of the original 
variables, the zeroth-order equations for the two-dimensional flows in the thin boundary layers 
that lie along the upper and lower surfaces of each blade and in the thin wake that extends 
downstream from the blade trailing edge have the form 


d£ 

dt 


+ = 0 , 


(2.28) 


and 



+ ^-(R«)- 


d_ 

dn 



+ = 0 


(2.29) 


-_DH t dP 
9 Dt~ Hi 


r^Tn 


(£ - (/>r)-'t) + (Pr)-'% dHT 


dn 


d ,~jt 


+^(M't<) = 0 . (2.30) 
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Here 


Dim = | 


„ 5 T ~. d d 

+ (! + «") ’'% + '% “ft 


+ ^l + ^ 4 - <2 - 3i) 


H T = H + V 2 /2 & H + V? /2, and K and V n are the fluid velocity components in the positive 
r- and n-directions, respectively. The viscous layer equations (2.28)-(2.30) are based upon 

the approximations 6 << 1, V • n ~ 0(6) and d/dn ~ 0(6 ), where 6 is the viscous-layer 

displacement thickness. These equations are written relative to a fixed reference frame, but the 
tangential and normal coordinates are m easur ed relative to moving blade and wake surfaces. 

The turbulent correlations v' T v' n and h' T v' n , appearing in equations (2.29) and (2.30), are 
related to gradients of the mean-flow variables, using Prandtl’s mixing-length hypothesis; i.e., 
we set 


pvy n = -e(Re) 


dn 


and phjv' n = — (i?e) 


-l 


~ 8Ht /- ~ \fr 

+ (f " e/,)Vr 


dVr 

dn 


(2.32) 


Here e and are the eddy viscosity and eddy diffusivity, respectively, and play roles similar 
to their molecular counterparts. The eddy diffusivity is related to the eddy viscosity, i.e., 
tfj = ( Prj)~ 1 e, through the introduction of the turbulent Prandtl number, Prj , which has 
a value near one for the flows of interest in turbomachinery applications. Values for e are 
provided via the use of an algebraic turbulence model; e.g., see [CS74] or [BL78]. 

It follows from the momentum equation (2.6) and the viscous layer approximations, that 
the component of the pressure gradient normal to the shear layer surface is given by 


n • V X P = 


dP_ 

dn 


- DV ~_, t72 , 

-pn ■ + - - • « PkV t + 


(2.33) 


Since the right-hand-side of (2.33) is of 0(1) and the thickness of the viscous layer is 0(6), the 
change in pressure across this layer must be of 0(6). Hence, in the zeroth-order approximation, 
the pressure in the viscous layer can be set equal to the inviscid pressure at the shear layer 
surface. Therefore, as a consequence of the high Re assumption, the pressure in thin viscous 
layers is a function only of r and t. Also, the pressure within and the flow properties at the 
edge of each viscous layer are equal to the inviscid values of these variables at the reference 
shear layer surface. 

The behaviors of the remaining flow variables in the viscous layer are determined by solving 
the mass (2.28), momentum (2.29), and energy (2.30) equations subject to appropriate surface 
and edge conditions. Also, since the system of viscous-layer field equations is parabolic in time 
(t) and in the streamwise (r— ) direction, the density, p, surface velocity, V T , and total enthalpy, 
H t , must be known for all time at some upstream location, and these variables, along with 
the normal velocity, must be known throughout the solution domain at some initial time. 

A no-slip condition and either a prescribed temperature or heat flux condition must be 
imposed at the solid blade surfaces, i.e., 

V = 'ft , xGB m or n = 0,T< tte , 

and _ . (2.34) 

H t = H T , w (r,t) or ^ - V T ^ = -k~\Re • Pr)~ 1 Q w for x € B m , 
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and a condition on the normal component of the fluid velocity, i.e., 

(V-£).n = 0 for x 6 W m or n = 0 , r > r TE (2.35) 


must be enforced at the reference wake surfaces. At the edge(s) of the viscous layer we require 
that 

V T -+ V Tie (T,t) and Hr — ► #r,e(M) for n — » oo , (2.36) 


where the limit n — * oo refers to the edges of the upper and lower surface boundary layers 
and the upper and lower edges of a wake. Here the subscripts w and e denote the values 
of the fluid properties at a solid wall and at the edge of the viscous layer, respectively, and 
the subscript TE refers to a blade trailing-edge. The fluid velocity, V T>e , and total enthalpy, 
Hr,e, at the edges of the viscous layers are determined by the inviscid solution along the blade 
and reference wake surfaces. Note that the exact locations of the reference wake surfaces 
are unknown a priori ; however, to within lowest order, the wake boundary conditions can be 
referenced to any arbitrary surface emanating from a blade trailing edge and lying within the 
actual viscous wake [Vel80]. 

The displacement, 6, and momentum, 0, thicknesses of a viscous layer are needed to de- 
termine the effects of viscous displacement and wake curvature, cf. (2.23)-(2.25), on the outer 
inviscid flow. For the two-dimensional unsteady flows, considered herein, these parameters 
are defined by 


and 


Hr,t)=r(l-j^L) dn 

J ° V P e (Vr,e-n T )J 

~8{T,t ) = r ^ ( 1 - Y r ~ *7 ) dn , 

Jo p e (Vr, e -n T )\ v T , e -nJ ’ 


(2.37) 

(2.38) 


where the integrations are carried out 


along lines normal to a moving blade or wake surface. 


2.4 Linearized Inviscid Unsteady Aerodynamic Model 

Even with the simplifications resulting from the assumption of a high Reynolds number, 
attached flow, the unsteady aerodynamic problem still involves prohibitive computing times 
and computational resources. Because of this, the traditional approach has been to assume 
that the unsteady excitations are of small-amplitude, harmonic in time, and exhibit a phase- 
lagged blade-to-blade periodicity, cf. (2.1)-(2.4). These assumptions lead to time-independent, 
linear equations for the inviscid unsteady perturbation, that can be solved over a single, 
extended, blade-passage region, for the complex amplitudes of the first-order unsteady flow 
variables. Because they are governed by linear equations, the first-order unsteady motions that 
arise from various Fourier modes of unsteady excitation are not coupled. Therefore, solutions 
for arbitrary excitations can then be obtained by a superposition of solutions for fundamental 
harmonic excitations. These features can lead to useful and efficient theoretical predictions 
of the unsteady aerodynamic response information needed in turbomachinery aeroelastic and 
aeroacoustic design studies. 

Thus, we regard the unsteady flow in the inviscid region as a small perturbation of a non- 
linear mean or steady background flow. Since the unsteady excitations are of small amplitude 
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[i.e., of O(e) << 1], the unsteady part of the inviscid flow can be approximated as a first-order 
(in e) perturbation of an underlying nonlinear background flow, that is steady in the blade- 
row frame of reference. In this case, the first-order unsteady fluid motions are governed by 
linear equations with variable coefficients that depend upon the steady background flow. Also, 
since the unsteady excitations are harmonic in time and the equations that govern the first- 
order flow properties are linear, these properties will have a harmonic time-dependence. We 
will take advantage of this feature by introducing complex representations for the first-order 
flow variables; thereby removing explicit physical time dependence from the linear unsteady 
problem. 

To determine the inviscid small-disturbance equations we first expand the dependent flow 
variables into asymptotic series of the form 

V[x(x, t ), t] = V(x) -|- v[x(x, t), t) + . . . = V(x) + #e[v(x) exp (iwf )] + . . . . (2.39) 

Here V(x) ~ 0( 1) and v[x(x,t),t] ~ O(e) are the velocities in the steady background flow 
at x and the first-order unsteady flow at x = x + i?e{R(x) exp (fort)}, respectively, v is the 
complex amplitude of the linearized unsteady velocity, u; is the temporal frequency of the 
unsteady motion, |R| ~ O(e), and the dots refer to higher order terms. 

The asymptotic expansion (2.39), with a field-point displacement field 7?.(x, t ) that satisfies 
the condition H{x,t) = n Bm {x,t) for x € has been employed in recent linearized Euler 
analyses [HC93a, HC93b, KK93, MV94]. It leads to a rather complicated set of unsteady 
field equations, but allows blade-surface boundary conditions to be imposed directly at the 
blade surfaces. However, since the blade motions are of small amplitude, we can set x = x, in 
(2.39), i.e., K = 0, which then becomes identical to the asymptotic expansion used in earlier 
linearized analyses [Whi87, Ver92, Ver93]. The latter expansion will be applied herein. It 
leads to a more convenient set of linearized unsteady field equations and precludes the need 
to prescribe a displacement field 72.(x,f) over the solution domain. However, Taylor series 
expansions, e.g., 

Vs = [v + (72. • V)V + . . -] s (2-40) 

must be applied to refer surface information from moving blade, wake and shock surfaces 
(5) to the mean positions (S) of such surfaces. Thus, in the present linearization the spatial 
coordinates x and x, refer to the same fixed location in the cascade frame of reference. 

The field equations that govern the zeroth-order steady and the first-order unsteady flows 
are obtained by substituting the asymptotic expansions for the flow variables [e.g., (2.39)] 
into the governing nonlinear field equations; equating terms of like power in e; and neglecting 
terms of second and higher order in e. Conditions on the unsteady perturbation at the 
mean blade, wake and shock surfaces are obtained by substituting the asymptotic (2.39) and 
Taylor (2.40) series expansions, and the relations between the unit tangent and normal vectors 
at corresponding points on the instantaneous (5) and mean surface (5) positions, i.e., 

t s m = r Sm + ^n Sm • n s m + - • • and n 5m = n Sm - ^n Sm • r Sm + - • • , (2.41) 

into the full time-dependent surface conditions, subtracting out the corresponding zeroth- 
order conditions and neglecting terms of higher than first order in e. These procedures lead to 
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nonlinear and linear variable-coefficient equations, respectively, for the zeroth- and first-order 
flows. The variable coefficients that appear in the linearized unsteady equations depend upon 
the steady background flow. 

Note also, that as a consequence of the assumptions regarding cascade geometry, the inlet 
and exit mean-flow conditions, and the temporal and circumferential behaviors of the unsteady 
excitations, the steady background flow will be periodic from blade-to- blade and the first-order 
unsteady flow will exhibit a phase-lagged, blade-to-blade periodicity. Thus, for example, we 
can write 

V(x + mG) = V(x) and v(x + raG) = v(x)exp(irocr) . (2.42) 

Such conditions allow numerical resolutions of the steady and linearized unsteady flows to be 
limited to a single extended blade-passage region. 

Steady Background Flow 

The field equations for the steady background flow are [cf. (2.20), (2.21) and (2.22)] are 

V, • (jSV) = 0 , (2.43) 

?(V ■ V,)V + V S P = 0 , (2.44) 

and 

V • v x s = 0 . (2.45) 

Here p , V, P and S are the density, velocity, pressure and entropy, respectively, in the steady 
background flow. The steady pressure can be determined in terms of the dependent variables 
p and S , using the steady form of (2.14), i.e., 

P = [Pp“' ,, exp(-75)] Ref p Y exp(75) (2.46) 

Surface conditions for the steady background flow follow from (2.23) through (2.25) and 
are imposed at the reference or mean positions of the blade, wake and shock surfaces. Mean 
blade positions are prescribed, but the mean wake and shock locations must be determined 
as part of the steady flow solution. 

The conditions 

V ■ n = p~ l d{p < .V f ^S)ldf , for x € B m , (2.47) 

apply at the mean blade surfaces. In addition, the zeroth-order normal velocity and pressure 
must satisfy the conditions 

[VI • n = (p- l d{p e VfJ)ldf) and [PJ = K(p e V? e (S + 0)) , for x 6 W m , (2.48) 

along blade wakes. Note that in the limit Re — ► oo, 5 — » 0 and 6 — > 0; therefore, the fluid 
normal velocity approaches zero at blade surfaces and the normal velocity and pressure become 
continuous across blade wakes. 

Finally, the requirements of mass, momentum and energy conservation provide the follow- 
ing conditions on the mean flow vaxiables at the mean shock locations, x € Sh m>n : 

[M f ] = 0, M / [V] + [P]n = 0 , and M } \Et\ + [PV] • n = = 0 , (2.49) 
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where Mj = pV ■ h. It follows from (2.49) that the tangential component of the fluid veloc- 
ity, V T , must be continuous across shocks. In principle, the shock-jump conditions must be 
imposed in nonlinear inviscid calculations, but the usual practice is to capture shocks by solv- 
ing conservative forms of the inviscid field equations across shocks using special differencing 
techniques. 

Steady-state, non-reflecting, inflow and outflow, boundary conditions (e.g., see [Gil90]) 
must be imposed on the steady background flow far upstream and far downstream from the 
blade row. In particular, averaged values of the inlet total temperature, total pressure and 
flow angle and the exit static pressure are specified, and circumferential harmonics of the mean 
flow variables, caused by the interaction of the mean flow with the blade row, are determined 
as part of the mean-flow solution. In addition, blade-to-blade periodicity conditions, cf. (2.42), 
can be imposed to restrict the steady flow solution domain to a single extended blade passage 
region. Since the inlet and exit conditions are imposed at finite distances from the blade row, 
say at £ = £_ and = <f + , respectively, the numerical solution domain is further restricted to 
a single extended blade passage region of finite axial extent. 

Linearized Unsteady Flow 

The differential equations that govern the first-order or linearized inviscid unsteady flow 
in continuous regions can be written as 

iujp + V* • (pV + pv) = 0 , (2.50) 

+ p{V • V X )V + = 0 (2.51) 

and 

^ + (v ■ V X )S = 0 . (2.52) 

Here p , v, p and s are the complex amplitudesof the time-dependent first-order density, ve- 
locity, pressure and entropy, respectively, and Dj Dt = iu + V • V is a convective derivative 
operator based on the mean flow velocity. To complete this system we require an additional 
equation that relates the first-order density, pressure and entropy. This is obtained by ex- 
panding the thermodynamic relation (2.14) and using (2.46) to obtain 

p = p[7 -1 p/P — s] = A~ 2 p — ps = 0 , (2.53) 

where A is the speed of sound propagation in the steady background flow. 

The thermodynamic relation (2.53) can be used to eliminate the density p from the lin- 
earized momentum (2.51) and continuity (2.50) equations, thereby expressing these equations 
in useful alternative forms. After performing the necessary algebra, we find that the linearized 
momentum and continuity equations can be written as 

— (v - sV/2) + [(v - sV/2) • V X ]V + V*{p/p) = V(v • V*S)/2 + p-'pV*S , (2.54) 

and 

^-[p/(M 2 )-s] + r I V i -(pv)=0, (2.55) 


Dv , 
« +(v 


V S )V 
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respectively. The linearized unsteady flow can then be determined by solving the system 
of convection equations (2.52), (2.54) and (2.55), subject to appropriate surface and far-field 
conditions, for the complex amplitudes, s, v and p, of the first-order unsteady entropy, velocity 
and pressure. Note that, if the mean flow is isentropic, the right-hand side of (2.54) is zero. 

Conditions on the inviscid unsteady perturbation at the reference or mean positions of the 
blade, wake and shock surfaces are obtained by substituting the asymptotic (e.g., (2.39)) and 
Taylor (e.g., (2.40)) series expansions and the surface vector relations (2.41) into the full time- 
dependent surface conditions (2.23)-(2.26), subtracting out the corresponding zeroth-order 
conditions (2.47)-(2.49), and neglecting terms of higher than first order in e. We have also 
made use of the steady field equations (2.43)-(2.45) and jump conditions (2.49) to write the 
first-order surface conditions in more convenient forms. We find that the linearized unsteady 
flow tangency condition can be expressed as 

[v — twR — VfdRj df + (R • V*)V] • ri = — p + t • dR/df) V • h + iuip~ 1 (p8 + p8 ) 


+p~ 1 r • V* [(pvl + p\8 + pV8 + 81 l • V*(pV)) • r] , x on B m , 

(2.56) 

where <5(x, t ) = £(x) + iie{<$(x) exp(iu;i)} + . . . and 6 and 6 ~ C?(|R|<5) are the steady and the 
complex amplitude of the first-harmonic components of the viscous displacement thickness, 
respectively, and R(x), x 6 B m , is the complex amplitude of the unsteady blade displacement. 
As a convenience, we have omitted the subscript e on the densities and velocities appearing 
in the surface conditions (2.56) and those given below. But, it is to be understood that the 
values of these inviscid fluid properties at a blade or wake surface are equal to their viscous 
values at the edges of the corresponding viscous layers. 

The linearized conditions on the jumps in the normal velocity and the pressure across 
wakes follow from (2.24) and (2.25), and the corresponding zeroth-order conditions, and have 
the form 

[vj-n - [ [V^]5R/df + (R- Vx)[Vj ] • n = fv] • n + Rr9[Vfi]/5f 


- d{RnlV?])/df - Rn[(V • V) lnpj - *£*[%] 

= 1(PP~ 1 + t ■ df ) V] • n + (iup~ l {pi + p8)) (2.57) 

Q 

+ (p- 1 + />V6 + pV8 + 81 l • V*(pV)) ■ r]) , x on W m . 

and 

[pj + (R-V JC )|[P] = |[p] + fi ? 5[Pl/5f-i? fi [p(V.V x )V].n 

= *(pV f \6 + e)) + K(V^(p + R-Vxp)V f + 2p(v f + (R-V 1t )V f )](8 + d)) (2.58) 

+ (pVf(8 + 0) (nVf - iutdFU/df)) , x on W m , 

where k = t dh/df and k are the steady and the complex amplitude of the first-order 
unsteady wake curvatures. 
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As do the steady, the first-order unsteady blade and wake conditions also simplify con- 
siderably in the inviscid limit Re — > oo. In particular, the right-hand-sides of (2.56)-(2.58) 
become zero. Also, for Re — > oo the left-hand-sides of the first-order wake conditions can be 
simplified by making use of the inviscid forms of steady wake-jump conditions and the field 
equations (2.43) and (2.44). After performing the necessary algebra, we find that 

I- v}.h = d(RnlV f l)/d? + R fi lV f d(\np)/dTl X € W m , (2.59) 

and 

M = -*RnlpV f % X € W m . (2.60) 

Equations (2.57) and (2.58) provide two independent relations for determining the jump in the 
linearized unsteady normal velocity and pressure across each wake. However, since the wake 
normal displacement, R • n, x € W m , is unknown a priori , these relations are not sufficient 
to determine [p] and |[v] • n, unless the steady tangential velocity, V • r, and density, p, are 
continuous across wakes. 

The linearized equations that ensure that mass, momentum and energy are conserved 
across shock discontinuities are 

[m/1 = [m /cl = [pv + pV - ia>pRj • n - d(Rnlp]V f ) / df = 0, x on Sh m , n , (2.61) 


M } [v + [< d(RnVn)/dT + i?nC] r J + ra/cIKln + [pJ n 

(2.62) 

-kRnlp(V* - V?)]h - RnVfipdVn/dfjh = 0 , x on Sh m , n 

and 

MAer + p/pj +mfAE T ] + K-Pv! • n - d(R*V f \P\)ldf - R*Vf\pdE T ldf\ = 0, x on Sh m , n , 

(2.63) 

where £ = V x x V and e T = e + V-v. Equations (2.61), (2.62) and (2.63), with the first-order 
mass fluxes, m/ c and m / defined by 

m/ c = (pv + pV — iwpR) • n — d(pRnV f )/df , 
and (2.64) 

m/ = m/ c — RpRnVn + Rfd(pVn)/df 

respectively, are the relations needed for determining the jumps in the first-order fluid prop- 
erties across moving shocks and the normal component of the shock displacement R • n, x € 
Sh mtn . These conditions along with the zeroth-order conditions (2.49) must be enforced to 
ensure that mass, momentum and energy are conserved to within first-order across moving 
shocks. Again, however, there is one more unknown associated with the unsteady shock- 
jump conditions, than there are independent equations. Therefore, these conditions are not 
sufficient for determining the relevant unsteady shock information. 

The first-order density and specific total internal energy can be eliminated from the first- 
order shock jump conditions by applying the thermodynamic relations (2.53) and 

e = p -1 [7 -1 p + (7 - = (7 /») -1 P + 7 -1 (7 - 1) _1 A 2 s (2.65) 
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to (2.61)-(2.64). By so doing, we would retain the convention, adopted in the derivation of 
the field er ations, of regarding pressure, entropy and velocity as the dependent variables of 
the lineariz i unsteady flow problem. However, since the jump conditions, derived above, are 
not sufficient for fitting wakes and shocks into an unsteady solution, based on the linearized 
Euler equations, additional information will be required. 

The foregoing equations provide linearized surface conditions, in which viscous displace- 
ment and wake curvature effects are taken into account. Such conditions are needed to account 
for viscous effects in a linearized inviscid analysis of the unsteady perturbation of a nonlin- 
ear steady background flow. The flow tangency condition (2.56) applies at the mean blade 
positions, and also on the upper and lower sides of the mean wakes; the jump conditions 
on normal velocity and pressure (2.57) and (2.58) also apply at mean wake positions; and 
the mass, momentum and energy conservation conditions (2.61), (2.62) and (2.63), respec- 
tively, apply at the mean shock positions. Unfortunately, these surface conditions are quite 
complicated and, to date, they have not been fully incorporated into a linearized unsteady 
aerodynamic analysis. Indeed, as presently posed, the wake and shock conditions are not 
sufficient to determine the jumps in the flow variables across a wake or shock surface and the 
surface normal displacement. Inviscid forms of the flow tangency condition (2.56) have been 
used successfully in linearized Euler calculations in which wake and shock effects are captured 
[HC93a, KK93, MV94], but there have not been any attempts to include the viscous terms 
on the right- hand- side of (2.56) in such calculations. 

In addition to the foregoing surface conditions, phase-lagged periodicity [cf. (2.42)] and 
far-field conditions must be imposed on the linearized unsteady flow. The latter must allow 
for the prescription of incoming entropy, vorticity and pressure disturbances at the inflow 
boundary and incoming pressure disturbances at the outflow boundary of the computational 
domain. In addition, unsteady disturbances coming from within the solution domain must 
pass through the computational inflow and outflow boundaries without distortion or reflection. 

It should be noted that to be precise, we have presented the foregoing steady and linearized 
unsteady equations in terms of the independent variable x, the mean-surface coordinates t 
and n, and the mean-surface, unit vectors r and n. However, since the displacement field 
TZ. = 0 and, as a result, the surface coordinates and unit vectors in the resulting steady 
and linearized unsteady equations always apply to the mean surface locations, we could have 
omitted the overbars in presenting the governing equations. To simplify the nomenclature, we 
will adopt the latter strategy in describing the LINFLO and SFLOW-IVI analyses in § 3 and 
§ 4, respectively. 

2.5 Discussion 

We have presented an inviscid/viscid interaction model for two-dimensional unsteady flows, 
occurring at high Reynolds numbers, in which the unsteadiness is driven by excitations of small 
amplitude. Although we have not yet developed solution procedures for the complete model, 
we have developed and evaluated solution procedures for several of the components needed 
for a complete inviscid/viscid interaction analysis of unsteady cascade flows. In particular, we 
have developed efficient flow solvers for linearized inviscid unsteady flows, for steady flows with 
strong inviscid/viscid interactions, and for unsteady flows with weak inviscid/viscid interac- 
tions. In constructing these analyses we have restricted our consideration to flows in which 
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any shocks that might occur are of weak to moderate strength and in which the free-stream 
flow conditions far upstream of blade rows are uniform. In such cases, the steady background 
flows in the inviscid regions can be regarded as isentropic and irrotational. For potential 
mean flows, the inviscid wake- and shock-jump conditions become well-posed, in that, there 
are a sufficient number of conditions to determine the flow properties at wakes and shocks. In 
addition, the potential mean flow assumption leads to two-dimensional, steady and unsteady, 
aerodynamic analyses that are very efficient computationally. 

In the following sections of this report we will describe the component unsteady aerody- 
namic analyses mentioned above. In particular, in § 3 we will describe the linearized inviscid 
analysis LINFLO, which applies to flows in which the unsteadiness can be regarded as a small 
perturbation of an isentropic and irrotational mean or steady background flow. The SFLOW- 
IVI analysis for steady flows with strong inviscid/ viscid interactions will be described in § 4. 
Finally, the unsteady viscous layer analysis UNSVIS, which can be used in conjunction with 
LINFLO to predict unsteady flows with weak inviscid/viscid interactions, will be presented 
in § 5. 

To demonstrate these analyses, we will apply them to three of the cascades studied in 
previous investigations [VH90, VBHA91, BVA93a] — a compressor exit guide vane (EGV), a 
high speed compressor (HSC) cascade, known as the Tenth Standard Cascade Configuration 
[FS83, FV93], and a turbine cascade, known as the Fourth Standard Cascade Configuration 
[FS83]. The blades of the EGV and HSC cascades are constructed by superposing the thickness 
distribution of a modified NACA four-digit series airfoil on a circular-arc camber line. The 
thickness distribution is given by 

T(x) = tfr[2.969x 1/2 - 1.26x - 3.516x 2 + 2.843x 3 - 1.036x 4 ], 0 < x < 1 . (2.66) 

where Ht is the nominal blade thickness. The coefficient of the x 4 term in (2.66) differs from 
that used in the standard NACA airfoil definition (i.e., —1.015) so that the example blades 
close in wedge-shaped trailing edges. The camber distribution is given by 

C(x) = H c - R + [R 2 - (x - 0.5) 2 ] 1/2 , 0 < x < 1 , (2.67) 

where H c (> 0) is the height at midchord and R = (2 H c )~ l {Hl + 0.25) is the radius of the 
circular-arc camber line. Thus, the surface coordinates for the reference blade are given by 

[X, Y] b,± = [iT 0.5T(x) sin 0 , C (x) ± 0.5T(x) cos 6] , 0 < x < 1 , (2.68) 

where 9 = tan _1 (dC/dx). 

The blades of the compressor exit guide vane (EGV) are constructed by setting Ht = 0.12 
and H c = 0.13. This cascade has a stagger angle © of 15 deg, a blade spacing G of 0.6 
and operates at an inlet Mach number and inlet flow angle of 0.3 and 40 deg, respectively. 
The blades of the high speed compressor cascade are constructed by setting Ht = 0.06 and 
He = 0.05. This cascade operates at a high-subsonic inlet condition, i.e., M- = 0.7 and 
fl-oo = 55 deg, and has a blade spacing and a stagger angle of unity and 45 deg, respectively. 

As a representative turbine configuration we have selected the Fourth Standard Config- 
uration, but for the present study we have modified the blade profiles, defined in [FS83], so 
that our example blades close in sharp, wedge-shaped, trailing edges. The turbine cascade 
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operates at an inlet , xh number of 0.19 and an inlet flow angle of 45 deg. The blade spacing 
G is 0.76, and the stagger angle 0 is 56.6 deg. In addition to the foregoing “real” blade 
cascades, we will also consider unsteady flows through flat-plate cascades. Here, the blade 
mean positions are aligned with the mean inlet flow direction, i.e., 0 = fl_oo; therefore, the 
local steady Mach number, M = M_ and flow angle, fl = fl_oo? are constants throughout 
the flat-plate flow fields. 
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3. The Linearized Inviscid Analysis: LINFLO 


In this Section, we will describe the linearized inviscid unsteady aerodynamic analysis, 
called LINFLO, in which unsteady disturbances are regarded as small perturbations of a 
potential steady background flow. The latter assumption leads to very efficient unsteady 
aerodynamic response predictions for realistic cascades and mean-flow operating conditions. 
The LINFLO analysis applies to unsteady flows excited by prescribed blade vibrations and/or 
external entropic, vortical and acoustic excitations. However, in the following discussion we 
will emphasize the application of LINFLO to unsteady flows excited by vortical gusts, as this 
was the major capability developed under contract NAS3— 25425. 

For more details on the development and application of LINFLO to subsonic, transonic 
and supersonic unsteady flows excited by prescribed blade motions (the flutter problem) we 
refer the reader to [VC84, Ver89a, UV91, MVF94]. Applications of LINFLO to unsteady 
flows excited by prescribed acoustic disturbances and to controlling the noise produced by 
wake/blade row interactions can be found in [KV93b] and [KV93a, KV94], respectively. Fi- 
nally the LINFLO analysis has been applied recently to help validate modern, time-accurate, 
Euler and Navier-Stokes analyses of unsteady cascade flows. These applications are described 
in [AV94] for unsteady subsonic and transonic flows excited by blade vibrations, and in [DV94] 
for subsonic flows excited by vortical and acoustic disturbances. 

3.1 Unsteady Perturbations of a Potential Mean Flow 

We consider inviscid unsteady flows through two-dimensional cascades in which the un- 
steadiness is caused by excitations of small-amplitude. We assume that the steady background 
flow fax upstream of the blade row is at most a small, isentropic and irrotational, perturbation 
of a uniform freestream, and any shocks that occur axe at most of weak to moderate strength. 
Because of these assumptions, the steady background flows will be isentropic (with 5 = 0) 
and irrotational; i.e., V x V = 0. Thus, we can set V = and (V ■ V)V = V(V$) 2 /2, 
where $ is the steady velocity potential. We also assume that the steady flow far downstream 
is a small perturbation from a uniform stream. In this case, analytical representations for 
the entropic, vortical and acoustic excitations at inlet and the acoustic excitations at exit can 
be provided, cf. (2.2)-(2.4). Also, closed form solutions can be determined to describe the 
unsteady entropic and vortical perturbations throughout the flow field. Therefore, numerical 
field methods are required only for determining the first-order unsteady pressure fluctuations. 

Note that in the present development 72 = 0 throughout the field; therefore, the spatial 
coordinates x and x are identical, both referring to a fixed point in the blade-row frame of 
reference. Thus, as a convenience, in this and the following chapter, we will simply use V to 
indicate the gradient operator with respect to x or x, 72 to denote a surface (blade, wake, 
or shock) motion and r and n to denote unit tangent and normal vectors, respectively, at a 
mean surface location. 

The overall procedure for determining linearized unsteady solutions is first to determine 
the steady flow for a given cascade configuration, characterized by G, 0, the blade shape 
and the inlet and exit flow conditions; and then, to determine the unsteady flow through this 
cascade for an unsteady excitation at a prescribed amplitude, frequency, and circumferential 
wave number, n n = cr/G, or interblade phase angle, cr. Because of the assumptions used in 
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developing the linearized unsteady aerodynamic model the solution domains for the nonlinear 
steady and the linearized unsteady problems can be restricted to a single extended blade- 
passage region of finite extent in the axial flow direction. 

The Steady Background Flow 


For a potential steady background flow the mass conservation equation (2.43) reduces to 


V • (pV$) = 0 


(3.1) 


and the momentum equation (2.44) can be integrated and combined with the thermodynamic 
relations for a perfect gas to yield the following (Bernoulli) relations for the mean flow variables 


= p = h - 1 )Ml„T 


= 1 - - 1 ] = 


2 + (7 — Wl<* 

2 + (t-1 )M ! ' 


(3.2) 


Here M and A are the local Mach number and speed of sound propagation, respectively, in 
the mean or steady background flow and 7 is the specific heat ratio of the fluid. 

Surface conditions for the inviscid zeroth-order or steady background flow apply at the 
mean positions, B m , W m and Sh m>n , of the blade, wake and shock surfaces. Since, by assump- 
tion, the flow remains attached to the blade surfaces, a flow tangency condition, cf. (2.48), 
which for inviscid flow has the form 


V$ • n = 0 , for x € B m , (3.3) 

applies at such surfaces. In addition, the steady normal velocity component and pressure must 
be continuous across the blade wakes (Wm), i.e., 

[VJ • n = 0 and [P] - 0 , for x G W m , (3.4) 

respectively. Since entropy and vorticity changes across shocks (Sh TOt „) are regarded as negligi- 
ble, the mass, momentum and energy conservation laws cannot all be enforced at such surfaces. 
The usual practice is to require only that mass and tangential momentum be conserved, i.e., 

[pVJ ■ n = 0 , and V] ■ r = 0 , for x € Sh m , n , (3.5) 

respectively. 

Numerical procedures for determining two-dimensional, steady, potential flows through 
cascades have been developed extensively, e.g., see [WN85, Cas83, HV93, HV94], particularly 
for flows with subsonic relative inlet and exit Mach numbers (i.e., M T<X> < 1). In such calcu- 
lations far-field boundary conditions are imposed at axial stations placed at finite distances 
upstream and downstream (i.e., at £ = £ T ) from the blade row, where linearized solutions 
describing the behavior of the steady potential can be matched to a nonlinear near-field so- 
lution. In addition, a Kutta condition is usually imposed at blade trailing edges in lieu of 
prescribing an exit freestream flow property. Finally, the usual practice is to solve the con- 
servative form of the mass-balance equation (3.1) throughout the entire fluid domain. Thus, 
the shock- and wake-jump conditions are not imposed explicitly. Instead, shock phenomena 
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Figure 3.1: Mach number contours and blade surface Mach number distributions for steady 
flow at M_oo = 0.3 and fl_oo = 40 deg through the EGV cascade. 

are captured through the use of special differencing techniques. Wake conditions are satisfied 
implicitly, because, in a potential flow, the fluid properties are continuous and differentiable 
across wakes. If needed, mean shock and wake (i.e., the downstream stagnation streamlines) 
locations can be determined a posteriori from the resulting steady flow solutions. 

Steady flows through the EGV and turbine cascades, that will be used here to demonstrate 
the vortical gust prediction capabilities of the LINFLO analysis, are illustrated in Figures 3.1 
and 3.2, respectively. These solutions have been determined using the methods of [Cas83]. 
In each case a Kutta condition has been applied at blade trailing edges; therefore, only inlet 
uniform flow information, e.g., M_ and fl_oo, have been specified for the steady calculations. 
The predicted steady Mach number field and Mach number distribution along a blade surface 
for the steady flow at M = 0.3 and = 40 deg through the EGV are shown in Figure 3.1. 
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Figure 3.2: Mach number contours and blade surface Mach number distributions for steady 
flow at M_oo = 0.19 and (l-oo = 45 deg through the turbine cascade. 

The calculated exit Mach number and exit flow angle for this flow are 0.226 and —7.4 deg, 
respectively, and the mean lift force, F y , acting on each blade is 0.360. The predicted steady 
Mach number contours and blade-surface Mach number distributions for the steady flow at 
M - oo = 0.19 and Q_oo = 45 deg through the turbine are shown in Figure 3.2. For the turbine 
cascade, the calculated exit Mach number and flow angle are 0.49 and 72.0 deg, respectively, 
and the mean lift acting on each blade is —2.11. 
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The Linearized Unsteady Equations 


The field equations that govern the first-order unsteady perturbation of an isentropic and 
irrotational steady flow can be expressed as a system of coupled differential equations for the 
complex amplitudes of the first-order entropy (s), velocity (v) and pressure (p), respectively 
(see [Gol78, Ver87]). In general, we require a solution to this system subject to prescribed flow 
conditions far upstream and far downstream from the blade row, flow tangency conditions at 
moving blade surfaces, and, since the inviscid field equations apply only in continuous regions 
of the flow, jump conditions at moving shocks and blade wakes. 

As indicated by Goldstein [Gol78, Gol79] the system of field equations that governs the 
linearized unsteady flow can be cast into a very convenient form by decomposing the unsteady 
velocity into rotational (vr) and irrotational (V<£) parts. The rotational velocity, vr, is taken 
to be divergence-free far upstream of the blade row, i.e., V • vr = 0 for £ < and the 
unsteady pressure depends only upon the potential 4> through the relation p = — pD<f>/ Dt , 
where D/Dt = iuj + • V is a convective derivative operator based on the mean flow 

velocity. The unsteady vorticity is given by £ = V x vr. For realistic blade profiles Vr ■ n and 
hence, V<£-n, will be singular along the mean blade and wake surfaces. Therefore, Atassi and 
Grzedzinski [AG89] introduced a modified form of the Goldstein decomposition, by setting 

v = v» + V<f> — Vfl + V<£« + V<£ , (3.6) 


to facilitate the numerical resolution of the velocity potential, <f>. Here, <f>* is a convected or 
pressure-less potential (i.e., D<f>./Dt = 0) which satisfies the condition V</>„ • n = —vr ■ n at 
blade and wake mean positions. 

The system of field equations that governs the linearized unsteady flow variables, s, vr and 
<t> is determined by substituting S = 0, V = and the velocity decomposition, Eq. (3.6), 
into the linearized Euler equations (2.52), (2.54) and (2.55) to obtain 



(3.7) 


JUv* - sV*/2) + [(v* - sV$/ 2) • V]V$ = 0 (3.8) 

and 

~ f "' v ’ = p v ' Wvb + = > r ' v ■ <^ v -> ■ (3 - 9) 

These equations are coupled only sequentially; hence, they can be solved in order to determine 
the complex amplitudes, s, vr and <f>, of the unsteady entropy, rotational velocity and velocity 
potential, respectively. Moreover, closed form solutions [Gol78, HV91] can be determined 
for the entropy and rotational velocity fluctuations in terms of the prescribed entropy and 
rotational velocity fluctuations at inlet. The velocity potential, <f ) , is governed by convected 
wave equation with source term p -1 V • (pv„) and depends, therefore, upon the rotational 
velocity, the convected potential, the conditions imposed at wake and shock surfaces, and far 
upstream and far downstream of the blade row. Phase-lagged periodicity conditions [cf. (2.42)], 

</>(x + mG) = <^>(x) exp(imcr) (3.10) 
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are also imposed on the entropy, rotational velocity and velocity potential to reduce the 
computational domain to a single, extended, blade-passage region. 

Surface Conditions 

As a consequence of the small unsteady-disturbance approximation, conditions on the lin- 
earized unsteady perturbation at moving blade, shock and wake surfaces can be transferred to 
the mean positions of these surfaces (see [Ver87]), with the mean wake, i.e., the downstream 
stagnation streamlines, and shock locations being determined from the nonlinear steady solu- 
tion. Thus, the (inviscid) flow tangency condition, cf. (2.56), can be written as 

v • n = V<£ • n = [iwR + (V$ • r)(r • V)R - (R • V) V$] • n, x € B m (3-11) 


where the complex amplitude of the blade displacement, R(x), x € B m , is a prescribed 
quantity. In addition, since the irrotational steady velocity and pressure are continuous and 
have continuous derivatives across the mean-flow downstream stagnation streamlines, the wake 
conditions, (2.57) and (2.58), reduce to 

[v • n] = [V^] • n = 0 and [p] = [Dj/Dt] = 0 , x € W m . (3.12) 


Finally, if we neglect changes in entropy and rotational velocity across shocks, the con- 
servation laws for mass and tangential momentum yield the following linearized shock-jump 
condition for a shock that terminates in the fluid 


[p(V^ + Vfi - A" 2 ^V$)J • n = -[p][io;+ (V$ • r)r - V][([V$] • n)" 1 ^] ] 
- (V* • n) -1 [^>] r • V([p]V* • r) , x <E Sh TO , n . 


(3.13) 


Equation (3.13) provides a relation for determining the jump in the unsteady potential, [<£]], 
at the mean position of a shock. The shock displacement normal to the mean shock locus 
is then given by R • n = — ([V$J • n) _ 1 [<^J. Because the mean flow is potential, the mean 
flow variables are continuous across blade wakes, and we have assumed that s and V/* are 
continuous across shocks, the foregoing wake- and shock-jump conditions provide sufficient 
information for determining the jumps in the unsteady potential, and hence, the remaining 
unsteady flow variables, across wake and shock surfaces. 


Far-Field Behavior 


We have assumed that the mean or steady flow is at most a small (i.e., of 0{t)) perturbation 
from a uniform stream both fax upstream (£ < £_) and far downstream (£ > £ + ) from the 
blade row. Therefore, in these regions, the first-order (in e) unsteady field equations can be 
reduced to constant coefficient equations for which analytical solutions can be determined (see 
[Ver89b]). 

For example, it follows after replacing by in (3.7) and (3.8), that far upstream 
of the blade row, the complex amplitudes of the entropy and rotational velocity fluctuations 
have the form (see also (2.2) and (2.3)) 

s(x) = s_oo exp(i/c_ 00 • x), £ < i- , (3.14) 
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and 

Vfi(x) = v R ,_oo exp(iK_ 0O • x), £ < Z- . (3.15) 

Here, s_oo and vr,_ <*, are the complex amplitudes (at x = 0) and K - oo is the wave number of 
the entropic and rotational velocity disturbances. The wave number vector has a component 
«jj,-oo = K -oo ’ e n = crG~ l in the cascade “circumferential” or 77 -direction and a component 
Kt = K-oo • ej = —uiVZ^o = — u? in the inlet freestream direction. Therefore, we can write 


K—oo — — (w sec fl-co “I - ffG tan Q — oo)^£ T e^ 

= —ojex + (wtanf2_oo 4- crG~ l secfI_oo)e;v , 


(3.16) 


where, and e„ are unit vectors pointing in the axial (£-) and circumferential ( 77 -) directions, 
respectively, and ej and e^r axe unit vectors pointing in the inlet freestream and normal to 
the inlet freestream directions, such that e z = e T x e N points out from the page. 

The complex amplitudes, s-00 and vr,-<x> • &Ni of the entropy and the rotational velocity 
component normal to the freestream direction are prescribed quantities. The component of 
the gust velocity in the inlet freestream direction, vr • e T , is determined by the divergence- 
free or orthogonality condition, iK -00 • Vfi.-oo = 0. The complex amplitude of the vorticity far 
upstream of the blade row is given by 

C -00 = (V X Vr)_oo = iK -00 X Vh,-oo = iK-oo v R,-°° ‘ ejv/«r,-oo , (3-17) 


and, since the vectors K -00 and Vr-oo are orthogonal, we can write 

VR- 00 = *(«-oo x C-oo)/l*-oo | 2 • 


(3.18) 


The entropy and vorticity, or rotational velocity at inlet are essentially prescribed functions; 
at exit, these quantities are determined as part of the unsteady solution 

The velocity potential fluctuations in the far upstream and far downstream regions depend 
upon the prescribed acoustic excitation as well as the acoustic response of the cascade. To 
describe these fluctuations, we can set 

<t>(x) = <f> E {x) + <M X ) for f > & , (3.19) 


where the potential component 4>e accounts for acoustic excitations, i.e., pressure disturbances 
that either attenuate as they approach the blade row or propagate and carry energy towards 
the blade row. In particular, it follows from (3.9) with V3> = V Too that for an acoustic 
excitation at temporal frequency u>, 4>e has the form (see [Ver89b]) 

4>e(x) = 4>l^oo exp[/S^oo( + O • x], £ > £* > ( 3 - 20 ) 


where 

4 >I , qpoo — P^loifl^ooV^oo COS fl^oo — *[w + ( K T<» ’ ( 3 . 21 ) 

the complex amplitudes p/ iT<x > are prescribed, = &/G, and ft^oo and /c^oo depend 

upon the inlet/exit freestream conditions, and the temporal frequency (u?) and interblade 
phase angle (<r) of the acoustic excitation. The potential component 4>r is associated with the 
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acoustic and vortical response of the blade row and therefore, must be determined as part of 
the unsteady solution. 

Usually only acoustic excitations that are of propagating type are considered. For subsonic 
inlet and exit conditions (M Tco < 1) the velocity potential corresponding to a propagating 
acoustic excitation at temporal frequency u> > 0 has the form (3.20) with /? Too = 0, and 

K rj,:foo — ~ — 4.T°° > (3.22) 

where the k^oo are the circumferential wave numbers at which cut-off or acoustic resonance 
occurs. These wave numbers axe given by 

43oo = wU“^M T oo (1 - Af^ 0O )" 1 (M Too sin n ?00 ± a/i - cos 2 fl Too ) . (3.23) 

It follows from (3.9), that the axial wave number of the propagating acoustic excitation is 

K Z&oo — "F|d^:oo| A/^oo > (3.24) 

where 

MtooI = 1(1 — -^TOO COS 2 fl Too ) 1 K 2 — •Af^ 00 6^ 00 | 1 / 2 , (3.25) 

and 

^oo = (wV^, + K v sin n Too )/ (1 - M cos 2 H Too ) . (3.26) 

Analytic solutions to (3.9), with V4> = V TOO , for the far-field potential component <j>n 
which satisfy the requirements that acoustic response disturbances either attenuate with in- 
creasing axial distance from the blade row or propagate carrying energy away from or parallel 
to the blade row, and that vorticity must be convected downstream are given in [Ver89b]. 
These solutions contain arbitrary constants that are determined by matching the far-field 
analytic solutions for the velocity potential to a near field numerical solution. 

3.2 Linearized Unsteady Solutions 

At this point we have presented the linearized equations that govern a general inviscid 
unsteady perturbation of a potential steady background flow, and we have indicated the 
quantities, i.e., u ;, cr, Rb, ^-oo, "Vr,-oo'^n and pi ,^ 0 0 that must be prescribed to determine the 
unsteady perturbation. We proceed to describe the procedures currently used in the LINFLO 
analysis to determine the unsteady entropy, rotational velocity and velocity potential. Once 
the unsteady potential is determined, it is a simple matter to determine response information 
needed for aeroelastic and aeroacoustic design applications, i.e., the unsteady pressures acting 
on the blades and in the far upstream and downstream regions of the flow. 

Entropy and Rotational Velocity 

Closed form solutions for the linearized entropy and rotational velocity fluctuations can be 
determined in terms of independent variables that describe the steady background flow [Gol78, 
HV91]. For this purpose we introduce the Lagrangian coordinate vector 

X = Ae r + Ve N . (3.27) 
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Figure 3.3: Drift and stream function contours for steady flow at M_oo — 0.3 and O_oo — 
40 deg through the EGV cascade. 

Here, 

A (x) = x_-er+/ V _1 dT* and 'I'(x) = x_ -e^+ / p(e z x V)-dr (3.28) 

v yx_+[>t(x)-'J'(x_)]e N ^ x - 

are the drift and stream functions, respectively, of the steady background flow. Also, x_ is the 
position vector to the point of intersection (£_ , 77 — ) of the reference blade stagnation streamline 
and the axial line £ = £_, dr# is a differential element of arc length along a streamline, and 
dr is a differential vector tangent to the path of integration. 

Drift and stream function contours for the EGV operating at Af_oo = 0.3, D_oo = 40 deg 
and for the turbine operating at M _«> = 0.19, fi_oo = 45 deg are shown in Figures 3.3 and 
3.4, respectively. These results have been calculated from the steady solutions depicted in 
Figures 3.1 and 3.2. Far upstream of each blade row, the drift function contours are parallel 
to each other and perpendicular to the inlet freestream direction, but as the flow proceeds 
downstream, gradients in the steady velocity field produce distortions. For the EGV , these 
distortions are mild over most of the blade passage, but severe in the immediate vicinity 
of the blade and wake surfaces, where the drift function contours stretch downstream from 
the leading-edge stagnation point. In addition, to the severe distortions near the blades and 
wakes, the drift function contours for the turbine are highly stretched within a blade and wake 
passage, because of the high steady velocities at mid passage. 

To determine the closed form solutions for the entropy and rotational velocity, we first 
note that DX/Dt = (V • V)X = V.*, = e T and X -► x a s£-> -oo. The solution to the 
entropy transport equation (3.7), which satisfies the far upstream condition (3.14), is then 
given by 

s(x) = 5_oo eXp( 2 K-oo • X) , (3.29) 
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Figure 3.4: Drift and stream function contours for steady flow at M_oo = 0.19 and = 
45 deg through the turbine cascade. 

and the solution to the rotational velocity transport equation (3.8), which satisfies the far 
upstream condition (3.15), is 

vr(x) = [V(X • A- oo) + s-coVS^exp^K-oo • X) , (3.30) 

where 

^ — oo VR,— oo 5— ooV— qq/2 . (3.31) 

It follows from (3.30) that 

C = V x r R — V(i#c_oo -X) x [V(A_oo -X) + exp(i<c_oo -X) , (3.32) 

and C — ► C-oo exp(i/c_ OQ • x) = x va i _ 00 exp(i#c_ 00 • x) as X -» x. 

If the steady background flow stagnates in the vicinity of a blade leading edge, as it will 
for realistic configurations, the drift function will have a logarithmic singularity at the mean 
blade and wake surfaces, i.e., A — * aolnn as n — ► 0, where n is the normal distance from 
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the surface and a 0 is a real constant. As a result, the exponential function exp(i/c — CO X) 
will be indeterminate, and the normal component of the rotational velocity will be singular, 
i.e., Vr ■ n —* din -1 exp(iK._ 00 • X), where a\ is a complex constant, at such surfaces. We can 
remove this singular behavior from the blade and wake surface conditions that are used to 
determine the unsteady potential, <f>, by prescribing a convected potential of the form [AG89] 


<t > » = 


-iu> 1 A-oo ■ V-oo + exp(i#c_oo • X) , 


(3.33) 


where 


f(*) = 


u;~ 1 (/c_ 00 x A- oo) • e^Gcos 0_ e 
27r(l — ia.o<jj) 


sin 


2 tt[^(x) - <g(x_)] ‘ 

GCOS fi_oo 


(3.34) 


is a complex function that depends upon, among other things, the behavior of the mean flow 
in the vicinity of a leading-edge stagnation point. This choice of </>„ ensures that v» • n = 
(v/i + V^*) • n = 0 at blade and wake mean positions. 

After combining (3.30), (3.33) and (3.34), we find that the complex amplitude of the 
source- term velocity, v» = \r + V</>,, is given by 


v» = 


• X) + 


f dF _ 1( 

. 3 *-" (K -= 


>A-oo) • e z I ViP + s_ooV$/2 


exp(i/c_oo • X) . 

(3.35) 

It follows from (3.34) and (3.35) that v* behaves like s_ 00 V$exp(iK_ 00 -X)/2 in the immediate 
vicinity of the mean blade and wake surfaces, i.e., asn-+0. Thus, v„ • n = 0, but, if S— oo 4- 
the tangential component of the source-term velocity will be indeterminate at such surfaces. 
It would be useful in future work to construct a convected potential, </>„, that also removes 
this indeterminacy, thereby allowing more accurate numerical resolutions of unsteady flows 
excited by entropic disturbances. 

The velocities vr and v« depend upon A and ^ and the first partial derivative of these 
functions. Therefore, the complex amplitudes of the unsteady vorticity, £ = V x v R = 
V x v», and the source term, p -1 V • (pv„), in (3.9) depend also upon the second partial 
derivatives of A and Thus, an accurate solution for the nonlinear steady background flow 
is a critical prerequisite for determining the unsteady effects associated with entropic and 


vortical excitations. 

The complex amplitudes of the entropy, rotational velocity, vorticity, and source term 
velocity are readily determined once the values of the drift and stream functions and their 
spatial derivatives are specified over the single extended blade-passage solution domain. For 
this purpose it is convenient to use an /f-grid in which one set of mesh lines are the streamlines 
of the steady background flow, for resolving unsteady flows excited by entropic and vortical 
gusts. An H- grid which covers the solution domain, i.e., one which is bounded by the upstream 
and downstream axial lines ( = and two neighboring mean-flow stagnation streamlines, 
is appropriate. The locations of the latter are determined a posteriori from the solution for 
the nonlinear steady background flow. Once the boundaries of the H - grid are established, 
the locations of the interior grid points can be determined using an elliptic grid generation 
technique, as described in [VH90, HV91]. 

Because a streamline mesh is used, the drift function can be evaluated at each point in the 
computational domain by a straightforward numerical integration of (3.28). The procedure 
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Figure 3.5: Contours of the in-phase component of the unsteady vorticity for the EG V cascade 
subjected to vortical gusts with vr-oq • e n -- (1,0) and ut = 5. 

used in [VH90, HV91] is simply to specify the drift function along the far upstream bound- 
ary £ = £_, and then to evaluate this function along each streamline using a second-order 
accurate difference approximation. The derivatives of the drift and stream functions at a 
given grid point are determined using the finite difference operators developed by Caspar and 
Verdon [CV81]. Because the drift function is singular at blade and wake surfaces, one-sided 
difference approximations are used to evaluate the derivatives of this function at points on 
the mesh streamlines adjacent to these surfaces. 

Calculated vorticity and source term fields for unsteady flows through the example EGV 
and turbine cascades excited by vortical gusts with vr-oo • £n = (1,0), uj = 5, and a = 
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Figure 3.6: Contours of the in-phase component of the source term for the EGV cascade 
subjected to vortical gusts with vr,-<x> • &n = (1,0) and u = 5. 

— 7 r, — 2x and — 3x are shown in Figures 3.5 through 3.8. These results were determined 
by performing the unsteady calculations on a (155 x 40) streamline H-mesh. Contours of 
the in-phase component or real part of the unsteady vorticity and source term are shown in 
Figures 3.5 and 3.6, respectively, for the EGV operating at M_ = 0.3, = 40 deg. The 

wave-number magnitudes, |*C-oo|, associated with the gusts at a — — x, — 2x and — 3x are 
5.65, 10.71 and 17.08, respectively, and the arguments relative to the axial flow direction, 
a_oo = tan -1 (/c, i _ 00 /K{,-oo)f are -112.2 deg, -77.8 deg and -67.0 deg, respectively. The 
vortical gusts are distorted as they axe convected by the nonuniform mean flow through the 
EGV. The vorticity contours in Figure 3.5 and the source term contours in Figure 3.6 indicate 
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Figure 3.7: Contours of the in-phase component of the unsteady vorticity for the turbine 
cascade subjected to vortical gusts with Vr-oo • = (1,0) and a; = 5. 

that this distortion increases in severity, i.e., the vorticity and source-term contours are more 
severely stretched and re-oriented within the blade and wake passages with increasing values 
of |a|. The results in Figure 3.6 reveal the rather strong variations in the source term that 
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Figure 3.8: Contours of the in-phase component of the source term for the turbine cascade 
subjected to vortical gusts with vr,-<x> • ejv = (1,0) and u> = 5. 

occur over the extended blade passage solution domain, particularly for the gust at a = — 3tt. 

Similar results for the turbine cascade operating at M_ <» = 0.19, fl-oo = 45 deg are shown 
in Figures 3.7 and 3.8 respectively. The wave number magnitudes, |k_oo|, and arguments with 
respect to the axial-flow direction, a_ <», for the vortical excitations at a — —ir, — 2tt and -37 r 
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axe 5.07, 8.35 and 13.50 and —125.4 deg, —81.8 deg and —66.7 deg, respectively. As indicated 
in Figures 3.7 and 3.8, the vortical gusts are highly distorted as they axe convected past the 
thick, highly cambered turbine blades. The unsteady vorticity and source term contours for 
the gusts at a = — and a = — 37r are quite different from those for the gust at cr = — 7 r. 
For a — — 7r the rectilinear vorticity contours far upstream of the blade row evolve into bowed 
shapes as the gust is carried through the blade row by the mean flow. The vorticity contours, 
within the passage and far downstream of the blade row, for the gusts at o = — 2tt and 
cr = — 37r are close to being straight lines. These lie at substantially different orientations than 
the contours upstream of the blade row. The source term contours in Figure 3.8 are severely 
distorted for the turbine blade row from mid-blade passage to the downstream boundary of 
the solution domain, particularly for the vortical gusts at a = —27 r and a = — 37 t. Also, the 
source terms associated with the gusts at a = — 27r and a = — 37r have very large gradients 
within the blade passage and downstream of the blade row. These features make it difficult 
to determine an accurate numerical resolution of the unsteady potential, for these turbine 
unsteady flows. 


Velocity Potential 

The unsteady potential ( <f > ) is determined as a solution of the field equation (3.9) subject 
to conditions at the mean blade, wake and shock surfaces, and conditions in the far field. Flow 
tangency [cf. (3.11)] applies at the blade surfaces, the fluid pressure and normal velocity must 
be continuous [cf. (3.12)] across blade wakes, and mass and tangential momentum must be 
conserved [cf. (3.13)] across shocks. The velocity potential in the far field is given by (3.19); 
the potential due to an acoustic excitation at frequency w and circumferential wave number 
= cr/G, by (3.20). 

A numerical resolution of the linear, variable-coefficient, boundary- value problem for <j> 
is required over a single, extended, blade-passage region of finite axial extent. The field 
equation must be solved in continuous regions of the flow subject to the surface and far 
field conditions on the unsteady potential. In particular, the near-field numerical solution for 
the potential must be matched to far-field analytical solutions at finite distances (£ = £ T ) 
upstream and downstream from the blade row. Numerical methods for determining (f> for 
isentropic and irrotational (i.e., 5 = ^ = 0) unsteady, subsonic, transonic, and supersonic 
flows have been reported in [CV81, VC82, VC84, UV91, MVF94]. Such solutions apply to 
unsteady flows excited by prescribed blade motions and/or acoustic excitations at inlet and 
exit. Numerical solution procedures for unsteady flows excited by entropic and/or vortical 
gusts [VH90, HV91, VBHA91] have been developed and implemented only for subsonic flows. 
The development of such procedures for transonic and supersonic flows remains, therefore, as 
a subject for future research. 

Because of the stringent and conflicting requirements placed on computational meshes for 
cascade flows, a composite-mesh (see Figure 3.9), which is constructed by overlaying a polar- 
type local mesh on an if-type cascade mesh, has been adopted for determining the unsteady 
potential. The H mesh is used to resolve unsteady phenomena over the entire solution domain; 
the local surface-fitted mesh, to resolve phenomena in the vicinities of rounded blade leading 
edges and/or normal shocks. The cascade mesh facilitates the imposition of the phase-lagged, 
periodicity conditions [cf. (3.10)] and the matching of the analytic and numerical unsteady 
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Figure 3.9: Extended blade-passage solution domain and composite mesh used in LINFLO 
unsteady transonic calculations. 

solutions at the far upstream (£ = £_) and far downstream (£ = £+) boundaries of the 
numerical solution domain. Use of this mesh alone is often sufficient for resolving unsteady 
subsonic flows, and this has been the strategy applied for calculating the unsteady subsonic 
solutions for cascade/vortical gust interactions presented in this report. The local mesh allows 
an accurate modeling of the unsteady flow in the vicinities of blade leading edges and normal 
shocks. It is constructed so that two “radial” lines coincide with the predicted mean shock 
locus to provide upstream and downstream shock mesh lines for the accurate imposition of 
unsteady shock-jump conditions. 

Since the cascade and local body-fitted meshes differ topologically, a zonal solution pro- 
cedure for overlapping meshes has been adopted in [UV91] for determining the unsteady 
potential. In the region of intersection between the two meshes, i.e., the region covered by the 
local mesh, certain cascade mesh points are eliminated depending upon their location within 
the local mesh domain. The discrete equations are written separately for the cascade and 
local meshes and coupled implicitly through special interface conditions, resulting in a single 
composite system of finite-difference equations that describe the unsteady flow over the entire 
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Figure 3.10: Contours of the in-phase component of the unsteady potential for the EGV 
cascade subjected to vortical gusts with vr-oo • e w = (1,0) and u> = 5. 

solution domain. 

The finite-difference model used to approximate the unsteady equations on the cascade 
and local meshes has been described in detail in [CV81]. Algebraic approximations to the 
various linear operators, which make up the unsteady boundary-value problem, are obtained 
using an implicit, least-squares, interpolation procedure that is applicable on arbitrary grids. 
This procedure employs a nine point “centered” difference star at subsonic field points, and 
a twelve point difference star at supersonic points. At a blade boundary point a nine point 
one-sided difference star is used on the cascade mesh, whereas nine- or six-point one-sided 
stars are used on the local mesh. Normal shocks axe fitted in the local-mesh calculation 
by approximating the shock-jump condition (3.13) using one-sided difference expressions to 
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Figure 3.11: Contours of the in-phase component of the unsteady potential for the turbine 
cascade subjected to vortical gusts with • e jv = (1,0) and to = 5. 

evaluate the normal derivatives of the unsteady potential on the upstream (supersonic) and 
downstream (subsonic) sides of the shock. At those points along the shock mesh lines at which 
the steady flow is continuous (i.e., at points lying beyond the end of the shock), the condition 
[<£] = 0 is imposed. 
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The systems of linear algebraic equations that approximate the unsteady boundary- value 
problem on the cascade and local meshes are block-tridiagonal for subsonic flow and, because 
shocks are fitted, block-pentadiagonal for transonic flow. A subsonic solution on the //-mesh 
alone is determined using a direct block inversion scheme. Composite (cascade/local) mesh 
solutions axe determined using a different scheme. Because of the cascade/local mesh coupling 
conditions, the composite system of discrete equations contains a sparse coefficient matrix of 
large bandwidth. Consequently, special storage and inversion techniques must be applied to 
achieve an efficient solution. Once the composite system of unsteady equations is cast into an 
appropriate format, it can be solved using Gaussian elimination [UV91]. 

The calculated unsteady potential ( <j > ) fields associated with the interactions of vortical 
gusts with the example EGV and turbine configurations are depicted in Figures 3.10 and 3.11, 
respectively. In particular, contours of the in-phase components of the unsteady potential, 
Re{<f > }, are given for vortical gusts at \r,-oo * e n = (1,0), u> = 5, and a = —7 r, — 2tt and 
—37 r. Vorticity and source-term fields for the unsteady flows through the EGV are shown in 
Figures 3.5 and 3.6, respectively; those for the flows through the turbine, in Figures 3.7 and 
3.8. The potential solutions depicted in Figures 3.10 and 3.11 were determined on 155 x 40 
streamline meshes and show variations over a blade passage that are associated primarily with 
the source term on the right-hand-side of (3.9). Once the unsteady potential is determined, 
the complete, linearized, inviscid, unsteady flow problem is solved. 

3.3 The Inviscid Response 

At this point we have provided a linearized unsteady aerodynamic formulation that de- 
scribes the general first-order fluid-dynamic perturbation of an isentropic and irrotational 
mean or steady background flow. We have also outlined the solution procedures used in the 
LINFLO analysis to determine the unsteady entropy, rotational velocity and velocity poten- 
tial. Solutions to the linearized unsteady problem are required to determine the aerodynamic 
response information needed for aeroacoustic and aeroelastic applications, e.g., the unsteady 
pressure fields far upstream and fax downstream of the blade row, and the unsteady pres- 
sures acting at the moving blade surfaces. We refer the reader to [Ver89a, Ver92] for detailed 
derivations of other local and global unsteady aerodynamic response parameters that are used 
in aeroelastic investigations. 

Approximate solutions for the full, nonlinear, time-dependent flow properties are con- 
structed by superposing the results for the steady and the linearized unsteady flow properties, 
e.g., we can write 

P(x, t ) = P(x) + p(x, t ) + J2 Psh m ,n (x, t) + . . . , (3.36) 

where p(x) = Re{p(x) exp(iu>f)}. The first two terms on the right-hand-side of Eq. (3.36) 
account for the steady and the first-harmonic contributions to the time-dependent fluid pres- 
sure acting at the point x. The third term accounts for the anharmonic contribution to 
the time-dependent pressure caused by the motions of normal shocks, and is determined by 
analytically continuing the steady solution from the mean to the instantaneous shock lo- 
cations [Wil79, Ver89a]. Thus, the first-order anharmonic response depends on the jumps, 
e.g. in the steady values of the fluid properties at the mean shock locations and the 
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unsteady shock displacements normal to the mean shock loci, TZ- • n, x € S /i m , n • The regions 
of anharmonicity are confined to thin strips that contain the mean shock loci. 

The steady background flow is described by the velocity potential, $, which is determined 
in terms of a prescribed inlet Mach number (M_ <»), flow angle (fi-oo)? and cascade geometry. 
The steady velocity, V = V$, is determined from this potential, as are the steady values of 
the thermodynamic properties of the fluid [cf. (3.2)]. The total enthalpy, Ht — T + V 2 / 2, in 
the steady background flow is constant and is given by 

Ht = Ht,-oo = (7 ~ + 2 (3.37) 

The linearized unsteady flow is described by the time-dependent, first-order, dependent 
flow variables 5 , v„ and 72-sh • n, which are determined in terms of prescribed values of the 
frequency u>, interblade phase angle a, and the complex amplitudes, Rb, -s_oo, C -00 an< l Pi,t- 00 ■> 
of the imposed unsteady excitation. The first-harmonic unsteady velocity is v = v„ + V<£, 
and the first-harmonic thermodynamic properties of the fluid can be determined in terms of 
the independent variables s and <f>. 

It follows from the equation of state for a perfect gas [cf. (2.11)] that 

p = 7 -1 ( 7 — 1)(M + pH) , (3.38) 

where h is the linearized unsteady component of the fluid specific enthalpy. After combining 
(3.38) with the fundamental thermodynamic identity relating the first-order entropy, pressure 
and density [cf. (2.53)] and recalling that p = -pDj>/Dt = -p(d/dt + • V)£, we find that 

the first-harmonic flow properties are related by 

p/p + S = 7 - 1 p/P = (7 - 1 F'ik/H -s) = —A~ 2 D<i>/Dt . (3.39) 

In addition, the first- harmonic total enthalpy, hj, is given by 

h T = Ts — D~4>!Dt + V* • [v. + V# + . . . . (3.40) 

The Inviscid Flow Along Moving Blade and Wake Surfaces 

The values of the inviscid flow quantities along the moving blade and wake surfaces (5) 
must be available to perform an analysis of the unsteady flows in the viscous boundary layers 
and wakes. These can be determined in terms of their values at the mean surface locations 
(S) and the blade and wake displacements via Taylor series expansions, e.g., 

Ps=Ps+Ps+K-'VP\s + 'EPsh m . n + ... , (3-41) 

where 72-, x € S is the surface displacement vector. In addition, relations between the unit 
tangent and normal vectors at the instantaneous and mean surface positions, cf. (2.41), are 
needed to determine the velocity components at the moving surface from information available 
at the mean surface. Here, the subscript S refers to points x$ = Xs+72-s on the moving surface, 
and the subscript S, to the corresponding points Xs on the mean surface. 

If we now restrict our consideration to continuous, i.e., shock-free steady background flows, 
the unsteady velocity at a moving blade or wake surface is given by 

V 5 = [V + v + {K ■ V)V] S + . . . = [V$ + v, + Vj> A (K • V)V$] S + . . . (3.42) 
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Expressions for the tangential and normal components of the fluid velocity at a moving blade 
and wake surface follow from (2.41), (3.42) and the conditions • ns = 0 and v» • n$ = 0. 
We find that 

Vs • Ts = [V$ + v» + V<^ + (72. • V)V4>]s • ts + — (3.43) 

and 

fi'T? 

Vs • ns = [V<?!> + (72- • V) V $]s • ns + ... = • ns . (3.44) 

The thermodynamic properties at the moving blade and wake surfaces are determined using 
Taylor series expansions along with the thermodynamic property relations for the steady (3.2) 
and linearized unsteady (3.39) flows. After performing the necessary algebra, we find that for 
shock-free flows 



P s 


7 ]Ps 
P s 




( 


D(f> 

Dt 


+ n • v[(v$) 2 /2] 


+ ... . (3.45) 
s 


Finally, the total enthalpy of the fluid at a moving blade or wake is given by 


Ht,s — [Ht + hj 4- TZ. ■ ^ Hx ] s + — ■ (3.46) 

Since the total enthalpy is constant in the steady background flow, i.e., Hj = Hx,~ and 
• v» = s( V 4> • t) 2 / 2 at blade and wake mean positions, we can write 


H T ,s = H t ,-oo + [Ts - [Dj>/Dt\ + • (V^ + v.)]s + . . . 

= [H Tt - oo(1 + s) - D<j>/Dt + • V^] s + . . . . 


(3.47) 


Expressions for the anhaxmonic values of the time-dependent flow properties in the field 
and at moving blade and wake surfaces can be determined following the methods of [Wil79, 
Ver89a]. However, for the subsonic inviscid and viscous-layer flows that are of primary concern 
in the present effort, only the steady and first-harmonic values of the time-dependent inviscid 
flow properties are required. 


3.4 Numerical Results: Cascade/ Vortical Gust Interactions 

Response predictions will be given below for two-dimensional, unsteady flows excited by 
prescribed vortical gusts. The intent is to demonstrate further the capabilities added to 
the linearized inviscid unsteady flow analysis, LINFLO, under Contract NAS3-25425. This 
analysis will be applied to predict the unsteady pressure responses of the EGV cascade, and 
a flat-plate counterpart to harmonic, in time and space, vortical excitations at inlet. The 
two cascades operate at a mean inlet Mach number (M_ <») and flow angle (fl-oo) of 0.3 and 
40 deg, respectively, and have a blade spacing ( G ) of 0.6. The blades of the EGV cascade 
are staggered at 0 = 15 deg; those of the flat-plate cascade axe aligned with the mean inlet 
flow direction, therefore, 0 = fl_oo = 40 deg. The predicted steady Mach number field and 
blade-surface Mach number distribution for the steady flow through the EGV are shown in 
Figure 3.1; the drift (A) and stream ('P) function contours, in Figure 3.3. As noted previously, 
the calculated exit Mach number, exit flow angle and mean lift force, F y , acting on each EGV 
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blade are 0.226, —7.4 deg and 0.360, respectively. The steady background flow through the 
unloaded ( F y = 0) flat plate cascade is uniform with M = 0.3 and 17 = 40 deg. 

The unsteady flows are excited by vortical disturbances that axe convected through the 
blade row by the steady background flow. A single harmonic component of such a disturbance 
is described by specifiying the temporal frequency, u>, interblade phase angle, cr, and either the 
complex amplitude, C-ooi of the vortical excitation or, as is done here, the complex amplitude 
of the normal component of the rotational velocity, i.e., Vr^oo • ejv, that is associated with 
this vortical excitation. For a pure vortical excitation, the rotational velocity and vorticity 
axe given by (3.30) and (3.32), respectively, with $_«> = 0 and A- = vr.oo. The wave 
number k_<x, of this excitation is given in terms of the frequency o>, interblade phase angle 
cr, blade spacing G and inlet flow angle, 17_ o©, in equation (3.16). For the example flat-plate 
cascade and far upstream of the EGV cascade, v/j(x) = vr-<x> exp(tfc_oo -x) and the complex 
amplitude of the gust vorticity is £_oo = in 2 _ 0 0 vr-oo • ejv//cx,-oo- 

A vortical gust at a given frequency and interblade phase angle can be regarded as a 
harmonic component of a wake excitation from an upstream blade row. For example, con- 
sider an upstream row that has the circumferential blade (wake) spacing Gexc and moves 
at circumferential velocity Vexc^ where Vgxc is the wheel speed, relative to the reference 
row, in our case the EGV or corresponding flat-plate blade row. The circumferential wave 
number and temporal frequency of the first-harmonic component, i.e., the component at the 
blade passing frequency (BPF), of the wake excitation are /c,_oo = — 2x/Gexc — cr/G and 
u> = 27tVexc/Gexc = — ctVexc/G. Thus, a vortical gust at a = —tt and u = 5 represents 
the first haxmonic component of a wake excitation from an upstream blade row having twice 
the gap or one-half the number of blades (JVexc = -Nref/2) as the reference row. The tem- 
poral frequency = 5 corresponds to an upstream flow coefficient, i.e., ratio of mean flow 
axial speed to wheel speed, Vloo cos O-oo/Pexc, of 0.802. In this case the vortical gusts at 
a = — 7 rn, u> = 5n, n = 2,3, . . . , i.e., at n x BPF, describe the higher harmonic content of 
the wake excitation. The strengths of the various harmonics axe determined by their complex 
amplitudes, VR, n _oo, -ejv. Alternatively, the gusts at o = — 7rn, u = 5n, n = 2, 3,... , can 
be thought of as the first-harmonics of wake excitations associated with upstream blade rows 
having blade spacings of Gexc = 2 G/n or blade counts of Nexc — «Aref/2. 

In general, for the example EGV and flat-plate cascades operating at an upstream flow 
coefficient of 0.802, the interblade phase angle and temporal frequency of the mth harmonic 
of a wake excitation from upstream are 

a = Tr(2mG / Gexc) — and u> = 7rG _1 Vw(2.mG / Gexc ) = 5n (3.48) 

Thus, the coefficient n = 2m G/ Gexc = 2mA r E xc/AREF varies inversely with blade spacing, 
or directly with the number of blades in the upstream row, and directly with m. 

We will examine the response of the EGV and the corresponding flat-plate cascade to 
vortical gusts with vr-oo • £n = (1,0), co = 5n and a = — tt n. The wave-numbers, <», 
associated with these gusts have magnitudes |fC_oo| = 5.65n and arguments relative to the 
axial flow direction, a = tan -1 (/c„ ,_ oo/k« ,-oo), of -112.2 deg. The complex amplitudes of the 
vortical disturbances at inlet are C-oo = (0,6.384n). It should be noted that, with the conven- 
tions adopted herein, u> and a must be of opposite sign to model a realistic cascade/ vortical 
gust interaction. We have chosen to set u > 0, and therefore, a must be taken to be less than 
zero. 
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We will determine the linearized unsteady pressure responses caused by the interactions 
of vortical gusts at w = 5n, a = -mr n = 1,2, and 3 with the EGV and flat-plate cascades. 
We will also examine the behaviors of the unsteady aerodynamic lift (/ y ) and moment (m) 
acting on the reference (m = 0) blades of the EGV and flat-plate cascades for 0.1 < n < 5. It 
is helpful in interpreting these results to recall that for subsonic flows an unsteady excitation 
that produces an acoustic response in which all pressure waves attenuate with increasing axial 
distance from the blade row is termed subresonant. A superresonant (m, n) excitation is one 
that produces an acoustic response in which m and n pressure waves persist far upstream and 
far downstream, respectively, and travel away from the blade row. Acoustic resonance occurs 
if at least one acoustic response wave persists far upstream or fax downstream, but travels 
only circumferentially, i.e., along the blade row, relative to an observer moving at the mean 
flow velocity. 

The linearized unsteady solutions reported below have been determined by performing 
LINFLO calculations on a global if -type mesh consisting of 40 mean-flow streamlines and 
155 “axial” lines, and extending one axial chord upstream and downstream of the blade row. 
The mesh streamlines were packed near the blade and wake surfaces; the axial lines, near the 
blade edges. The global mesh analysis is very efficient; for example, a complete unsteady flow 
calculation on a 155 x 40 global mesh requires about 2 minutes of CPU time on an IBM 370 
work station. Most of the information determined from such a calculation can be saved and 
re-used in subsequent unsteady flow calculations for different frequencies, interblade phase 
angles, and modes of excitation, which then require only about 30 seconds of CPU time for 
each unsteady case. It should be noted that LINFLO composite-mesh calculations typically 
require four (4) times the CPU time required for a global mesh calculation. 

As part of the present study, response predictions, i.e., unsteady pressure difference distri- 
butions [Ap(x) = p(xg) — p(xj)] and unsteady lifts and moments, have also been determined 
for the flat-plate cascades using the classical linearized analysis of Smith [Smi72]. Recall that 
in the LINFLO formulation the unsteady potential equation (3.9) contains the source term 
p -1 V • (pv„), and the normal derivative of the potential, • n, is zero at the mean blade 
surfaces. For flat-plate blades aligned with the inlet flow direction the source term reduces 
to V 2 ^,. In the classical linearization the potential equation is homogeneous, and the normal 
derivative of the potential cancels the normal component of the gust velocity at the mean 
blade surfaces. Thus, a comparison between the LINFLO and classical linear solutions is 
very meaningful, because it reveals the ability of the numerical approximations used in the 
former analysis, to accurately account for the strong source term variations that occur over 
the solution domain. 

Results for the unsteady flows axising from the interactions of vortical gusts at vr-oo-gn = 
(1,0), a; = 5n, and a = im , with the EGV and flat-plate cascades are shown in Figures 3.12 
through 3.18. Contours of the in-phase component or real part of the unsteady vorticity 
for 07 — 5 and o = —i r and those for the in-phase components of the unsteady pressure 
responses to the vortical gusts at n = 1, 2 and 3 are shown in Figure 3.12 for the EGV and in 
Figure 3.13 for the flat-plate cascade. Surface pressure responses for the gusts at n = 1, 2 and 
3 are shown in Figures 3.14-3.16, and the behaviors of the unsteady lift and moment versus 
n for 0.1 < n < 5 are presented in Figures 3.17 and 3.18, respectively. 

The contours of the real part of the unsteady vorticity for the interaction of a vortical 
gust at a; = 5, <7 = — 7r, i.e., at say, the blade passing frequency (BPF), of a wake excitation, 
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Figure 3.12: Contours of the in-phase components of the unsteady vorticity and the unsteady 
pressure responses for the EGV cascade subjected to vortical gusts with vr-oo ’ &N = (1,0). 

and the contours of the real parts of the unsteady pressure responses of the EGV to vortical 
excitations at BPF, twice BPF and three times BPF are shown in Figure 3.12. Although, 
the complex amplitudes and wave number magnitudes |/c_oo| of the vortical gusts vary 
linearly with n, the vorticity contours associated with the excitations at 2BPF and 3BPF 
have very similar shapes to those shown in Figure 3.12 for the excitation at BPF. The vortical 
gusts are convected through the EGV blade row and distorted from the rectilinear pattern 
that exists far upstream by the nonuniform steady background flow. 

The interactions of the vortical gusts with the EGV blading produce unsteady pressure or 
acoustic responses. The vortical gust at u = 5 and a = —% (n = 1) produces a subresonant 
acoustic response, i.e., all pressure response waves attenuate with increasing axial distance 
from the blade row. The vortical excitation at 2BPF (n = 2) produces a superresonant (1,1) 
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Figure 3.13: Contours of the in-phase components of the unsteady vorticity and the u eady 
pressure responses for the flat-plate cascade subjected to vortical gusts with = (1, 0). 

response in which a pressure response disturbance, at u; = 10 and a\ = a + 2v = K Vtl G = 0 
propagates away from the blade row in both the far upstream and far downstream directions. 
The axial wave numbers, 1 ,^ 00 , of these far upstream and fax downstream pressure distur- 
bances are 3.895 and —2.440, respectively, and their amplitudes, |pi, TO o| are 0.119 and 0.120. 
The vortical excitation at 3BPF (n=3) produces a superresonant (1,0) response in which one. 
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acoustic response wave at u> = 15, 02 = <r + An = K Vt 2 G = n propagates far upstream of 
the blade row and all other waves attenuate with increasing distance from the EGV blade 
row. The propagating wave has a wave number vector with magnitude |/c 2 ,— oo| = 6.325 and 
argument, 02,-00 = 55.9 deg relative to the axial flow or positive ^-direction, and an amplitude 
IP 2 - 00 I = 0.070. 

Note that spurious pressure behaviors, as indicated in Figure 3.12, are predicted at the 
upstream and downstream boundaries of the solution domain for the vortical gust at n = 2 
and at the downstream boundary for the gust at n = 3. These behaviors can be attributed to 
the nonuniform spacings between the mesh streamlines in the global H-meshes used for the 
LINFLO calculations. The mesh streamlines must be packed near the blade and wake surfaces 
to resolve the vortical flows near these surfaces, but this packing leads to a loss in accuracy in 
matching the numerical near-field and analytical far-field solutions for the unsteady potential 
at the inlet and exit boundaries, and thus, in predicting the complex amplitudes of the far-field 
acoustic response waves. 

The vortical excitation at u = 5 and a = — ir and the unsteady pressure responses of the 
flat-plate blade row to vortical excitations at u> = 5 n and a = —n n are shown in Figure 3.13 
for n = 1,2 and 3. Note that for flat-plate cascades having blades aligned with the mean 
flow direction, the vortical gusts are convected through the blade row, without distortion, 
by the uniform steady background flow. Thus the rectilinear vorticity pattern that exists 
far upstream of the blade row is maintained throughout the entire flow field. For the flat- 
plate cascade, the vortical gust at BPF (n=l) produces a subresonant pressure response. 
The gusts at 2BPF and 3BPF produce superresonant (1,1) responses. At 2BPF acoustic 
waves at u = 10 and oj = a + 2n = = 0 propagate away from the blade row in 

the negative and positive axial directions. For the upstream propagating wave, k^i,_oo = 
3.895 and |pi,-oo| = 0.225; for the downstream propagating wave, k^i.+oo = -2.439 and 
|pi,+oo| = 0.167. The gust at 3BPF produces upstream and downstream propagating acoustic 
response waves at w = 15, <r 2 = a + 4n = k Vi 2 G = n. For the upstream propagating 
wave |k 2> — 00 1 = 6.325, a 2 ,-oo = 55.9 deg and |p 2 ,_oo| = 0.122, and for the relatively strong 
downstream propagating wave, |k 2 ,+oo| = 5.308, a 2 ,+oo = 99.5 deg and |p 2 ,+oo| = 0.424. 

The acoustic response behaviors of the EGV and the flat plate cascades to the excitation 
at uj = 15, a = — 3n differ in the far-downstream region. The reason for this is that different 
mean flow conditions prevail far downstream of the two blade rows, which, for the excitation 
at 3BPF, lead to very different far-downstream acoustic environments. Recall that the EGV 
operates at an exit Mach number and flow angle of 0.226 and —7.4 deg; the flat plate cascade, 
at M+ 00 = M = 0.3 and ft +00 = ft = 40 deg. A vortical excitation at u> = 15 and a = —Zn, 
is superresonant (1,0) for the EGV cascade; therefore, a propagating acoustic wave persists 
far upstream, but all acoustic response waves attenuate with increasing distance downstream. 
This same excitation produces a superresonant (1,1) acoustic response for the flat-plate blade 
row, with a relatively strong propagating acoustic response wave in the far downstream field. 

Surface pressure responses at the reference ( m = 0) blades of the EGV and flat-plate 
cascades are depicted in Figures 3.14-3.16. These results indicate the effects of gust distortion, 
due to nonuniform mean flow phenomena, on local unsteady blade loading. The surface 
pressure distributions along the EGV and flat-plate blades for the vortical gusts at BPF 
(Figure 3.14) and 2BPF (Figures 3.15) show somewhat similar qualitative behaviors, but the 
EGV results show a more wave-like character and there are important quantitative differences 
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Figure 3.14: Unsteady surface pressure responses of the EGV and corresponding flat-plate 
cascades to a vortical gust with vr-oo • = (1,0), u = 5 and a = — ir. 


between the unsteady loads acting on the two cascades. In contrast, because of the different 
acoustic response environments that exist far downstream of the two blade rows, the surface 
pressure responses of the EGV and flat-plate blades to the vortical gust at 3BPF (Figure 3.16) 
bear little qualitative resemblance. 
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Flat Plate 



Figure 3.15: Unsteady surface pressure responses of the EGV and corresponding flat-plate 
cascades to a vortical gust with v«_ • e s = (1,0), u> = 10 and a = — 2k. 


The complex amplitudes of the unsteady lift and moment acting on the reference EGV 
and flat blade blades for the unsteady flows considered in Figures 3.14 through 3.16 are listed 
below. The unsteady moment is taken about blade midchord. 

Vortical Excitation EGV Cascade Flat Plate Cascade 

vr-oo • e^r = (1,0) 

w=5, <7 = — 7T fy = 0.827,-1.100 fy = 0.732,-1.184 

m = -0.014, 0.142 m = -0.074, 0.197 

u = 10, a = — 2 tt f y = 0.122, -0.141 /„ = 0.223, -0.268 

m = -0.036, 0.053 m = -0.108, 0.115 

w = 15, cr = -3 tt f y = -0.010, 0.059 f y = -0.406, 0.074 

m = 0.013, 0.020 m = -0.066, 0.133 
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Figure 3.16: Unsteady surface pressure responses of the EGV and corresponding flat-plate 
cascades to a vortical gust with v^-oo • e^v = (1,0), u> = 15 and o = — 3;r. 

These results indicate large differences between the global unsteady airloads, due to vortical 
gusts, that act on the EGV and the flat plate blades. 

We should note that the blade pressure-difference distributions and the aerodynamic lifts 
and moments predicted by the LINFLO and Smith (CLT) analyses for the flat-plate flows 
considered in Figures 3.14—3.16 are in excellent agreement. Indeed, the pressure- difference 
curves predicted by the two analyses are almost coincident for the gusts at n = 1 and n = 2, 
and show only slight differences for the gust at n = 3. 

Predictions for the EGV and flat-plate lift and moment responses for vortical excitations 
at v fi _ oo • ejv = (1,0), o> = 5n and a = —irn are given in Figures 3.17 and 3.18. Here we 
consider the behaviors of the unsteady aerodynamic lift ( f y ) and moment (m) about midchord 
that act on the reference blades of the EGV and flat-plate cascades versus n, as n varies from 
0.1 to 5. LINFLO predictions are given for the EGV and both LINFLO and classical linear 
theory predictions are given for the flat-plate cascade. The LINFLO and classical linear 
theory predictions for the lift and moment responses are in excellent agreement. The lift 
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Figure 3.17: Unsteady lift versus n for the EGV and flat plate cascades subjected to vortical 
gusts at vh,-oo • e n = (1,0), u> = 5n, a — —ir n. 

forces (Figure 3.17) acting the EGV and flat-plate blades differ, but not substantially, over 
the entire range 0.1 < n < 5; however, there are more significant differences between the 
moments (Figure 3.18) acting on the blades of the two cascades. 
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Figure 3.18: Unsteady moment about midchord versus n for the EGV and the corresponding 
flat plate cascades subjected to vortical excitations at vr^-oo ■ eyv, = 5n and a = —irn. 

It is helpful for interpreting the lift and moment responses indicated in Figures 3.17 and 
3.18, to describe the features of the far-field acoustic response for the two cascades. The 
vortical excitations at n = 1.463, 2.688, 2.926 and 4.389 produce an acoustic response dis- 
turbance at a cut-off or acoustic resonance condition fax upstream of the EGV and both far 
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upstream and far downstream of the flat plate cascade. The excitation at n = 2.688 pro- 
duces a resonant acoustic response disturbance that travel along the blade row in the positive 
77-direction; those at n =1.463, 2.926 and 4.389 produce resonant disturbances that travel 
down the blade row, i.e., in the negative 77-direction. Vortical excitations at n =1.558, 2.865, 
3.115 and 4.673 produce a cut-off or resonant acoustic reponse disturbance far downstream 
of the EGV. The resonant response disturbance at n =2.865 travels upward along the blade 
row; those at n =1.558, 3.115 and 4.673 travel downward. As is readily determined from the 
curves in Figures 3.17 and 3.18, the unsteady aerodynamic lift and moment undergo abrupt 
and often large changes in the vicinity of an acoustic resonance condition. Such a condition 
also represents a boundary for different types of unsteady response behaviors. 

For example, the responses of the flat-plate cascade to the vortical excitations in the range 
0.1 < n < 5.0 are subresonant for n < 1.463 and 2.688 < n < 2.926; superresonant (1,1) 
for 1.463 < n < 2.688 and 2.926 < n < 4.389; and superresonant (2,2) for n > 4.389. 
Because the inlet and exit free-stream conditions for the EGV differ, this blade row has a 
wider variety of far-field acoustic response behaviors than its flat plate counterpart. For 
example, the acoustic responses of the EGV are superresonant (0,1) for vortical excitations 
at 2.688 <7 1 < 2.865; subresonant for 2.865 < n < 2.926; superresonant (1,0) for 2.926 < n < 
3.115 and superresonant (1,1) for 3.115 < n < 4.389. The different types of far-field acoustic 
response behavior seem to have a strong impact on the unsteady blade loads. In particular, the 
flat plate cascade has large amplitude moment responses to vortical excitations in the range 
2.688 < n < 2.926 which produce subresonant acoustic responses. The moment responses of 
the EGV to vortical excitations in this range are much smaller in magnitude. The acoustic 
responses of the EGV are superresonant (0,1) for 2.688 < n < 2.865 and subresonant for 
2.865 <n< 2.926. 

Discussion 

This completes our description of the LINFLO analysis and of the capabilities of this analy- 
sis, developed under Contract NAS3— 25425, to describe cascade/vortical gust interactions. In 
subsequent sections of this report we will describe a steady inviscid/ viscid interaction analysis 
and an unsteady viscous layer analysis that are being developed for use in conjunction with 
LINFLO to provide efficient prediction capabilities for unsteady viscous flows. We will also 
present LINFLO results for unsteady flows caused by prescribed blade motions and acoustic 
excitations. 

The LINFLO analysis represents a powerful analytical capability for efficiently predict- 
ing the unsteady aerodynamic information needed in blade row aeroelastic and aeroacous- 
tic response studies, in which many controlling parameters are involved. In recent stud- 
ies [DV94, AV94] LINFLO has been applied to assist in calibrating modern, time- accurate, 
Euler/Navier-Stokes analyses. Very good agreement between LINFLO predictions and non- 
linear Euler and Navier-Stokes predictions has been determined for unsteady subsonic flows 
excited by small-amplitude, vortical or acoustic disturbances [DV94] and for unsteady subsonic 
and transonic flows excited by prescribed blade vibrations [AV94]. Based upon the limited 
ranges of parametric studies that were conducted using the nonlinear analyses, it was also 
found that unsteady pressure responses are linear over surprisingly wide ranges of excitation 
amplitudes. Such results provide evidence on the usefulness of a linearized inviscid, unsteady, 
aerodynamic analysis for aeroelastic and aeroacoustic design studies. 
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4. The Steady Inviscid /Viscid Interaction Analysis: SFLOW-IVI 


Efficient analyses for predicting the effects of strong inviscid/viscid interaction (I VI) phe- 
nomena, due, for example, to viscous-layer separations, shock/boundary- layer interactions and 
trailing-edge/near-wake interactions, on the aerodynamic, aeroelastic, and aeroacoustic per- 
formances of turbomachinery blading are needed as part of a comprehensive analytical design 
prediction system. The focus here will be on the development of an accurate and efficient IVI 
analysis for steady cascade flows, that can provide the foundation for an unsteady procedure 
to be developed later. The steady analysis is described below and is being developed as part of 
an overall research program, which has the goal of providing reliable and efficient theoretical 
prediction methods for steady and unsteady viscous flows, at high Reynolds numbers, through 
subsonic and transonic cascades. 

The present steady analysis uses an inviscid/viscid interaction (IVI) approach, in which 
the flow in the outer inviscid region is assumed to be potential, and the flows in the inner or 
viscous-layer regions are governed by Prandtl’s viscous-layer equations. The inviscid equations 
are solved using a Newton iteration procedure and an implicit, least-squares, finite-difference 
approximation, similar to that used in LINFLO. The viscous-layer equations are solved using 
a modified Levy- Lees transformation and an inverse, finite-difference, space-marching method, 
which is applied along the blade surfaces and the wake streamlines. Complete details on the 
inviscid and viscous solution procedures can be found in [HV93, HV94] and [BVA93a], respec- 
tively. The inviscid and viscid solutions are coupled via a semi-inverse global iteration pro- 
cedure that permits the prediction of boundary-layer separation and other strong-interaction 
phenomena. Numerical results will be presented below for the three example cascades de- 
scribed in § 2.5 with a range of inlet flow conditions considered for one of them, including 
conditions leading to large-scale flow separations. Comparisons with Navier-Stokes solutions 
and experimental data will also be given. 

4.1 General Concepts 

For the flows of .Tactical interest in either internal or external aerodynamics, the Reynolds 
number is usually sufficiently high so that the flow past an airfoil or blade can be divided into 
two regions: an “inner” dissipative region consisting of boundary layers and wakes, and an 
“outer” inviscid region. The principal interaction between the flows in the viscid and inviscid 
regions arises from the displacement thickness effect which leads to thickened semi-infinite 
equivalent bodies v . corresponding changes in surface pressures. 

If the interact' 3 “weak,” then the complete flow problem can be solved sequentially. 
In this case the pr ire distributions along the blades and wakes are first determined by a 
pure inviscid solution. These distributions are then imposed when solving the viscous-layer 
equations to determine viscous displacement thickness, £(r), distributions along the blades 
and wakes. The latter are then used to obtain a new inviscid solution, that acco un ts for 
viscous displacement effects. The resulting changes in the blade pressure distributions and 
in the downstream freestream flow properties (e.g., Mach number and flow angle) can then 
be calculated. It is sometimes possible to continue this sequential solution procedure until a 
converged solution for the entire flow is achieved. 
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Flows over airfoils or blades, however, involve both a weak overall interaction arising from 
standard displacement thickness and wake curvature effects, and local strong-displacement 
interactions caused, for example, by viscous-layer separations, shock/boundary-layer interac- 
tions, and t railing-edge/ near- wake interactions. Viscous displacements in a strong-interaction 
region can cause substantial changes in the local inviscid pressure and, in some cases, in the 
overall pressure field. The concept of an inner viscous region and an outer inviscid region 
still applies, but the classical hierarchical structure of the flow breaks down. In particular, for 
flows with strong interactions the inviscid/viscid hierarchy changes from “direct”, in which 
the pressure is determined by the inviscid flow, to “interactive” in which the pressure is de- 
termined by the mutual interaction between the flows in the inviscid and the viscous-layer 
regions. This change must be accommodated within an inviscid/viscid interaction solution 
procedure. 

The construction of an IVI analysis involves, first, the development of component, i.e., 
inviscid and viscous-layer analyses, and, second, the implementation of these components into 
a strong-interaction iteration procedure to produce a solution for the complete flow field. 
Solution methods for steady subsonic and transonic inviscid flows through cascades (e.g., 
see [Cas83, HV94]) and for steady boundary-layer and wake flows (e.g., [VV85] and [BV87] ) 
have been developed to a relatively mature state. Methods for coupling such solutions have 
also been developed and assessed through a number of model problem studies (e.g., see [Vel80, 
VV85, BV87]). Inviscid/viscid interaction procedures for predicting steady flows in cascades 
have also been developed [HSS79, JH83, CH80, BV89, BHE91] and applied over a wide range 
of inlet flow conditions, including conditions leading to stall [BHE91]. 

In the present approach we consider high Reynolds number (Re — p*_ x V'^L* j ji *_ x ) steady 
flow, with negligible body forces, of a perfect gas with constant specific heats and Prandtl 
number through a two-dimensional cascade (cf. Figures 2.1 and 2.3). In particular, we will 
restrict our consideration to adiabatic flows at unit Prandtl number, however, it is a relatively 
simple matter to extend the analysis to heat conducting flows at arbitrary Prandtl numbers. 
The flow in the outer inviscid region is assumed to be isentropic and irrotational and hence, 
governed by the full-potential equation; those in the inner viscous-layer regions, by Prandtl’s 
viscous-layer equations. The non-hierarchical nature of strong interactions is taken into ac- 
count by the procedure used to couple the inviscid and viscid solutions. In addition, an inverse 
viscous-layer calculation, in which the displacement thickness is specified instead of the pres- 
sure, is employed to permit viscous solutions to be continued through local strong-interaction 
regions, including regions of separated flow. In regions of the flow where the inviscid/viscid 
interactions are weak, the pressure, as determined from the inviscid solution, can be imposed, 
instead of the displacement thickness, to obtain a direct viscous-layer solution. 

In the present approach, the steady, cascade, full-potential analysis (SFLOW) developed 
by Hoyniak and Verdon [HV93, HV94] is employed to determine the flow in the inviscid 
region. SFLOW has been constructed for use with the linearized inviscid unsteady flow anal- 
ysis, LINFLO, to provide compatible and comprehensive steady and unsteady inviscid flow 
prediction capabilities for cascades. In the IVI calculation procedure (which is referred to as 
SFLOW-IVI), viscous effects are incorporated by modifying the blade and wake boundary 
conditions used in SFLOW to account for the effects of viscous displacement. The nonlinear 
inviscid analysis, coupled with the IVI iteration procedure, allows nonlinear changes to the 
inviscid base flow due to viscous effects to be evaluated. The ability to treat such nonlinear 
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perturbations is especially important in transonic flows in which shock positions and strengths 
axe significantly altered by viscous displacements. Althoug the analysis described below is 
presently restricted to subsonic flows, it can be extended to treat transonic flows in the future. 

4.2 Inviscid/ Viscid Interaction Analysis 

Inviscid Region 

The flow in the inviscid region is determined as a solution of the field equations (3.1) 
and (3.2) subject to flow tangency conditions at the blade surfaces, jump conditions on the 
normal velocity and pressure across the blade wakes, and prescribed uniform flow conditions 
far upstream of the blade row. Blade-to-blade periodicity conditions [cf. (2.42)] are applied 
upstream and downstream of the blade row to limit the computational domain to a single, 
extended, blade- passage. Kutta conditions at the blade trailing edges and a global mass 
conservation condition that relates the flows at the inlet and exit boundaries and accounts for 
the blockage effects of the viscous layers are enforced in lieu of specifying the uniform flow far 
downstream of the blade row. For the flows considered here, the inlet and exit velocities are 
subsonic. 

The field equation (3.1) can be written in the form 

A 2 V 2 $ - • V(V$) 2 /2 = 0 (4.1) 

where the speed of sound propagation A is given in terms of the potential $ in (3.2). The 
SFLOW inviscid analysis is based on this non-conservative form of the full potential equation. 

The specific forms of the blade and wake conditions for the potential steady flow in the 
inviscid region can be determined by setting V = V4> in (2.47) and (2.48). Thus, at each 
blade surface 

• n = p~ l d{p e V Ttt 8)ldr , for x G B m , (4.2) 

where p e and V T<e = d$/dr | e are the inviscid density and velocity at this surface, or the 
viscous density and streamwise velocity component at the edge of the viscous layer, and 6 is 
the boundary-layer displacement thickness (see Figure 2.3). The quantities r and n denote 
the arc distance along the blade (positive in the downstream direction and zero at the leading- 
edge stagnation point) and the local unit normal vector directed outward from the surface, 
respectively. At the blade wakes the inviscid solution for the normal component of the fluid 
velocity and the pressure must be discontinuous with jumps given by 

[V$]]-n = (p~ 1 d{p e Vr,tS)ldr) and [PJ = R(p e V* e (8 + 0)) , for x € W m , (4.3) 

where n is the upward pointing unit normal vector to the reference wake streamline (i.e., W 
in Figure 2.3), (8) = 8w and {&) = 6w are the displacement and momentum thicknesses of 
the complete wake, and « is the curvature of the wake which is taken as positive when the 
reference wake streamline is concave upwards. 

As mentioned in § 2.3, a complication arises in that the location of the reference wake 
streamline is unknown a priori ; however, to within lowest order, the wake conditions can be 
referenced to any arbitrary curve emanating from the trailing edge and lying within the actual 
viscous wake [Vel80]. Usually, wake curvature effects axe regarded as negligible. In this case 


62 



the inviscid pressure, density and tangential velocity are continuous across blade wakes, and 
the wake conditions (4.3) reduce to 

[V$] • n = pd(p e V r , e Sw)/dT and [P] = 0, x G W m . (4.4) 

For steady flows, the reference wake streamline is taken to be the aft stagnation streamline as 
determined by a pure inviscid solution. This is adequate except in extreme cases where the 
location of the stagnation streamline is significantly altered by viscous effects. In this case, it 
is possible to periodically update the location of the wake streamline during the calculation, 
although this has not been done in the present study. 

The foregoing boundary- value problem is solved using a Newton iteration procedure. Thus, 
we set 

$n+l=$n + ^», »* = 0,1,2,... (4.5) 

where is the estimate to the final solution at the nth iteration level and <£ n is a correction 
to this estimate. The quantity is determined at each step of the iteration process by solving 
linear equations that are derived by substituting (4.5) into (4.1)— (4.3) and neglecting terms 
that are of second and higher order in 4> n . For example, the linear field equation for the 
correction, <f > n , is 

- (7 - ^ - V(V<S) ! -V} /2 = -^V 2 4>„ + V4>„ ■ V(V*„) 2 /2 , 

(4.6) 

where A n , the nth estimate to the speed of sound propagation, is determined from (3.2) with 
$ replaced by $ n and D/Dt - V$ n ■ V. The right-hand side of (4.6) is known and the linear 
operator on the left-hand side is readily derived from the linear unsteady operator appearing 
on the left-hand side of (3.9). The various equations used to determine <f> n are approximated 
using the same implicit, least-squares, finite- difference approximations as those used in the 
LINFLO analysis. The resulting linear system of algebraic equations that approximate the 
linear boundary value problem for 4> n is solved by direct matrix inversion, as in LINFLO, using 
lower-upper decomposition and Gaussian elimination. The Newton iterations are continued 
until 1 1 <fi n 1 1 < e, where || || denotes a prescribed norm and e << 1 is a user specified tolerance 

level. This nonlinear steady analysis, called SFLOW, is described in detail in [HV93, HV94]. 

The present inviscid solutions were obtained on a “streamline” type H-mesh, rather than 
on the “sheared” H-mesh described in [VC84] and [HV93]. The SFLOW analysis was modified 
by Hoyniak to use the streamline H-mesh developed by Hall and Verdon [HV91]. Thus, prior 
to initiating an I VI calculation, a pure inviscid solution is first obtained on a sheared H-mesh. 
The resulting solution is then used to generate a streamline H-mesh, in which one set of mesh 
lines corresponds to the streamlines of the inviscid flow, and the second set consists of lines 
that are “nearly” orthogonal to the first set. The principal advantage of the streamline H- 
mesh over the sheared mesh is an improved resolution of the flow in the vicinities of blade 
leading edges. 

An alternative to this procedure is available in SFLOW and was used for one of the cases 
described in this report, i.e., the turbine cascade. In the turbine case, a useful streamline 
H-mesh could not be determined from the solution on a sheared H-mesh, because the latter 
provided an inadequate resolution of the flows at blade leading-edges. To remedy this, the 
initial inviscid solution was obtained on a composite mesh constructed by overlaying a local, 


63 


surface-fitted, C-mesh on a global sheared H-mesh. A detailed description of this procedure, 
as applied to linearized unsteady flows, can be found in [UV91]. 

Viscous Layers 

The flows in the inner or viscous regions are assumed to be governed by Prandtl’s viscous- 
layer equations [cf. (2.28)-(2.32)]. For steady flows ( d /dt\ x = 0 and P = P, p = p, etc.) the 
continuity and streamwise momentum equations have the form 



We assume that the flow in the viscous layer is adiabatic and occurs at unit Prandtl number. 
In this case the energy equation, cf. (2.30), reduces to the requirement that the total enthalpy 
of the fluid, Ht « T + V T 2 /2, must be constant across the viscous layer. 

In equation (4.8) the subscript e refers to fluid properties at the edge of the viscous layer, 
and the effective viscosity, p, c s, is defined to be 


Ptff = P + e , (4.9) 

where p is the molecular viscosity, which is assumed to be a function of temperature alone, 
e oc 7xr is the turbulent eddy viscosity, and 7 j r is the streamwise intermittency factor. In the 
present study, the moleculax viscosity, //, is determined by Sutherland’s equation (2.10); the 
eddy viscosity in blade boundary layers, using the Cebeci-Smith model [CS74], as modified 
in [BV87] to account for flow separation; and the eddy viscosity in wakes, using the model 
of Chang, et al. [CBCW86]. The specific turbulence model, used in the present study, is 
described in detail in § 5.1. For the present IVI calculations, instantaneous transition, i.e., 
7 jr changes abruptly from 0 to 1, is assumed to occur at specified locations along the blade 
surfaces. 

The foregoing field equations govern the flow in the viscous layers along the upper and 
lower surfaces of the blades and in the blade wakes. They are solved subject to conditions at 
the edges of the viscous layers, on the blade surfaces, and along the reference wake streamlines, 
i.e., 

V r — * K, e for n — > oo , r > 0 , (4-10) 

Vr = V n = 0 for n = 0, 0 < T < (4-11) 

and 

K = 0 for n = 0, r > r* £ , (4.12) 

respectively, where t? e are the trailing-edge values of the upper- (+) and lower-surface (— ) arc- 
length coordinates measured from the leading-edge stagnation point. The condition expressed 
by (4.10) is also applied along a wake streamline for n — ► — oo. Equations (4.11) and (4.12) 
imply that the curve n = 0 corresponds to the blade surfaces and reference wake streamlines, 
respectively. 
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The displacement (<5) and momentum (0) thicknesses of the viscous layers are needed to 
determine the effects of viscous displacement and wake curvature on the outer inviscid flow. 
For steady flows, these quantities are determined by [see also (2.37) and (2.38)] 


and 



(4.13) 


(4.14) 


where the zero lower bound on the integrals is replaced with — oo when determining the 
displacement and momentum thicknesses of a wake. 


Levy-Lees Transformation 

To facilitate the numerical resolution of the viscous-layer equations we introduce a modified 
version [BV87] of the Levy- Lees transformation [Blo70]. Thus, we define new independent 

var 1 y| |^| pq 

i = l p e V T , e ii c f[,edT and rj = (Re/2£) ll2 V r , e f pdn , (4-15) 

Jo JO 

and new dependent variables 

F = V T /V T , e and / = (Re/ 20 1/2 ® - (4-16) 

Here £ and 77 are scaled streamwise and normal spatial coordinates, F is the ratio of the local 
to the edge value of the streamwise velocity, and / is a scaled stream function for the flow in 
the viscous layer. The quantity /x e ff,e('r) in (4.15) is the effective viscosity at the edge of the 
viscous layer. It is used in the definition of £ to maintain a nearly constant value of 77 at the 
edge of a turbulent viscous layer [VWV82]. The variable in (4.16) is the stream function of 
the flow; therefore, dV/dn = pV T and d^/dr = -pV n . The Levy-Lees transformation permits 
the leading-edge, stagnation- point, laminar, similarity solution to be easily recovered and leads 
to a reduced truncation error in the numerical approximation to the viscous-layer equations 
relative to that associated with an analysis based on the use of primitive flow variables. 

After applying the transformation relations (4.15) and (4.16), we find that the continuity 
(4.7) and momentum (4.8) equations can be written as 

F=?f (4.17) 

a Tj 


and 




+ /?(0-F 2 ) = 0. 


(4.18) 


Here l = pfi e fi/(pefie({,e), 0 = h/P, and /? is a pressure gradient parameter, defined by 


df/ T , e 

K,e ^ ' 


(4.19) 
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Since the total enthalpy Ht ~ T + V ? /2 and the pressure are constant across the viscous 
layer, it follows from the Bernoulli relations (3.2) and the equation of state for a thermally 
perfect gas, cf. (2.11), that 


0 = T/T, = (PIP,)-' = 1 + 1 - F 2 ) , (4.20) 

where T/T e and p/p e are the ratios of the local to the edge values of the fluid temperature 
(or enthalpy) and density. 

The following boundary conditions are applied. The stream function is constant along 
the blades and the reference wake streamlines (i.e., at t) = 0). Therefore, without loss in 
generality, we can set 

/ = 0 at T) = 0 . (4.21) 

The no-slip condition applies at a blade surface, i.e., 

F = 0 at Tj = 0 for £ < £ T e , (4.22) 

where 8te is the trailing-edge value of £. At the edges of the blade-surface boundary layers 

F(r} e ) = 1 for ( < £ te , (4.23) 

since V T — ► V T>e as r) — ► rj e . This condition forces the flow variables to approach their appro- 
priate edge (inviscid) values as rj — *• T) e . 

The edge conditions for a wake are more complicated because of the jumps in the inviscid 
flow variables, associated with wake curvature [cf. (4.3)]. We indicate the upper- and lower- 
edge values of t) and the edge velocity, V Tje , by the superscripts + and — . If we use upper 
surface flow variables to define the Levy-Lees transformation for the wake calculation, e.g., if 
we set F = V T /V+ e , then 

F(vt)= 1 and F(r,:) = V-JV ; ► for { > ( TE . (4.24) 

If the curvature effect is negligible, i.e., if fPJ 0 for x € It', the wake. edge conditions become 

Htf) = l. 

Solution Procedure 

The viscous-layer equations are parabolic in the ^-direction and therefore require initial 
conditions. These are provided by determining a similarity solution which holds in the vicinity 
of a leading-edge stagnation point, i.e., near £ = 0. Such a solution is obtained by solving 
the Levy Lees equations with (3 = 1. Solutions for the flows along a blade surface and its 
wake are then obtained using space- marching in the downstream direction. As discussed in 
§ 4.1, a complete IVI calculation requires the ability to solve the viscous-layer equations in 
both the “direct” mode, in which the pressure gradient parameter j3 is specified and the 
displacement thickness 8 is determined, and in the “inverse” mode, in which 8 is specified 
and /3 is determined. In the present study the equations are solved in the direct mode over 
a forward part of the blade, that includes the leading edge stagnation point, and in the 
inverse mode downstream of an axial station, whose location is either specified in advance or 
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determined during the calculation to ensure that the inverse mode is initiated upstream of a 
strong interaction region. A wake is calculated entirely in the inverse mode. For a calculation 
in the inverse mode, the quantities /? and V T>e are unknown. Thus, a supplemental equation 
relating these two variables is needed. This is obtained by discretizing (4.19), which defines 
j3 in terms of V T<e . 

In the direct viscous-layer calculation, the value of ft is determined by the inviscid analysis 
and the displacement thickness is obtained from the viscous analysis. In the inverse procedure, 
the displacement thickness is specified, and the edge values of the flow variables, V T>e , M e , etc., 
are obtained as part of the viscous-layer solution. This is accomplished via the introduction 
of a “mass deficit parameter”, m = p e V r>e 8. An expression relating the value of / at the edge 
of the viscous layer to m is derived by integrating equation (4.17) across the boundary layer 
and employing the definitions of 8, v and 6 , cf. equations (4.13), (4.15) and (4.20). We find 

that rvc _ 

f{Ve) = Odv- ( Re/20 1/2 m . (4.25) 

Jo 

Equation (4.25) is used at blade surfaces to impose the specified value of m through the 
corresponding value of / at the outer edge of the viscous layer, i.e., at v = 7e- The integral 
term in (4.25) is determined from the previous global IVI iteration, and is therefore specified 
or lagged in the current iteration. 

An expression for the difference between the stream function values at the upper and lower 
edges of a wake is obtained by integrating (4.17) across the entire wake (i.e., from v e to 7+ )> 
to obtain 

M) - fiv7) = f + 0d v - (Re) 1 ' 2 {m + /(2e) l/2 + fh-/(2C) 1/2 ] • (4-26) 

Equation (4.26) is used to impose the mass deficit parameters, = pf in the wake. 

If wake curvature effects are regarded as negligible the last term in (4.26) can be replaced by 
rriw / (2£w) 1//2 , where rn\v = pV T>e 8wj (2£w and is the continuation of the upper surface 

£- variable into the wake. The integral term in (4.26) is determined from the previous iterative 
solution. 

The discretized field equations, boundary conditions and auxiliary conditions, equations 
(4.17)-(4.26), are quasi-linearized and the resulting tridiagonal system of algebraic equations 
is solved at each streamwise (i.e., r- or £-) station of the computational domain, using a 
fixed-point iteration to update the nonlinear terms. The inversion algorithm used in the wake 
is modified to account for the application of one boundary condition (4.21) at »? = 0 and 
two others at the upper and lower edges of the viscous layer, as well as to account for the 
application of the jump condition (4.26) on the values of / at the upper and lower edges of 
the viscous wake. Finally, the so-called FLARE approximation, which prevents instabilities in 
the viscous-layer solution, due to flow reversal, is applied by turning off the convective terms 
in the momentum equation wherever F < 0. Further details on the viscous- layer numerical 
analysis can be found in [BVA93a]. 

Inviscid /Viscid Iteration Procedure 

The present IVI approach determines the complete flow field by iteratively updating the 
mass deficit parameter, m(r), which affects both the inviscid and viscous solutions through 
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Figure 4.1: Semi-inverse inviscid/ viscid iteration procedure. 


their respective boundary conditions. For an arbitrary m distribution, two different surface 
and wake streamline velocity distributions generally result: one, V Tyei [r), from the inviscid 
calculation and one, K, ev ( T )> f rom the viscous-layer calculation. The objective is to deter- 
mine a converged inviscid/viscid interaction solution by finding the mass deficit parameter 
distribution that minimizes the differences between the V T<ej and V Tytv distributions. 

In this investigation the “semi-inverse” iteration procedure of Carter [Car79], which is 
illustrated schematically in Figure 4.1, has been used to update m at every streamwise mesh 
station on the blade and wake surfaces. Thus, we set 

m"« = m”[l + u - 1)] , (4.27) 

where the superscript n is the global iteration count and u is a relaxation parameter. The 


68 














solution is considered to be converged when 


max IK.ev - V T ,, 


e/, 


\/Vr, eii <Z,i 




(4.28) 


where the value of e is specified by the user and IE is the number of streamwise mesh stations. 
Equation (4.28) is applied on both blade surfaces and along the wake. The viscous-layer 
solution is obtained at the locations corresponding to the intersections of the inviscid mesh 
lines with the blades and the reference wake streamlines. This avoids the need for interpolation 
between different inviscid and viscid streamwise mesh. 

During the global iterations, the independent variable £ is updated using equation (4.15), 
where the current values of the variable appearing in the integrand are applied. Because a 
major objective of this study has been to develop an efficient analysis, various techniques for 
accelerating convergence were examined. We found that one of the most effective approaches 
for reducing the CPU time needed to obtain a converged IVI solution is to use the largest 
value of the relaxation parameter, u>, for which the iterative procedure remains stable. It was 
also observed that the inviscid velocity distribution, V T , ei , changes relatively little between the 
initial purely inviscid solution and the final IVI solution. This is in contrast to the significantly 
larger changes observed in the viscous velocity distribution, V T , ev • This observation prompted 
the introduction of a sub-iteration loop in which the viscous equations are solved repeatedly 
during each global iteration. Thus, equation (4.27) is applied N v times during a single global 
iteration, with V T , e/ being frozen at its most recent value and K,e v being re-calculated during 
each sub-iteration by solving the viscous-layer equations using the latest m distribution. The 
value of Nv is a user-specified input, and the standard iteration procedure is recovered if 
JVy = 1 . This strategy is only effective in reducing the total CPU time if the number of 
global iterations needed to obtain a converged IVI solution, Nq, can be reduced enough 
to more than balance the increased computational effort needed for the additional viscous 
calculations performed at each global iteration level. Of the three cascades examined in § 4.3, 
the sub-iteration procedure was of benefit only for the turbine. The dashed lines in Figure 4.1 
correspond to the viscous subiteration technique described above. 


4.3 Numerical Examples: Steady Flows with Strong Inviscid/Viscid Interactions 

The foregoing inviscid/ viscid interaction analysis has been applied to the cascade configu- 
rations described in § 2.5; in particular, the compressor exit guide vane (EGV), the high-speed 
compressor (HSC) cascade known as the Tenth Standard Cascade, and the turbine cascade, 
which is a modified version of the Fourth Standard Configuration. Surface pressure coefficient, 
Cp = (P — P_qo)/2, or Mach number, M, displacement thickness, 6, and surface shear stress, 
t w = Re~ x JidV r ldn\ n= o, distributions will be presented, as functions of chordwise distance 
x, for the three cascades. IVI solutions for the compressor cascades will be evaluated via 
comparisons with Navier-Stokes solutions; the solution for the turbine against experimental 
measurements. In addition, predicted values of the total pressure loss, the exit flow angle, 
and separation point location will be presented for the EGV operating over a wide range of 
inlet flow conditions. These predictions were obtained using the mixing analysis of Stewart 
[Ste55]. Finally, the performance of the SFLOW-IVI analysis, i.e., its efficiency and conver- 
gence properties, will be discussed. 


69 


Figure 4.2: The EGV cascade and streamline H-mesh. 


In the calculations described below, the SFLOW-IVI analysis was applied using the norm 
(4.28) and a convergence tolerance, e, of 0.001. The inviscid H-meshes used for the two 
compressor cascades consisted of 90 axial and 31 circumferential lines, with 24 axial lines 
upstream of the blade leading edges, 41 lines intersecting the blade surfaces and 25 lines aft 
of the trailing edges. The meshes used for the viscous-layer analyses employed a total of 81 
and 25 streamwise grid lines along the blade and wake, respectively, with 71 normal grid 
lines across a surface boundary layer and 141 normal grid lines across the wake. The inviscid 
mesh used for the turbine cascade had 150 axial and 31 circumferential lines, with 39 points 
upstream of the leading edge, 51 points along each blade surface and 60 points along the wake. 
A total of 101 surface and 25 wake stations were used in the turbine viscous-layer analysis. 
The normal mesh had the same dimensions as those used for the compressor cascades. For 
the cases considered in this study, the wake curvature effect was assumed to be negligible; 
thus, {P} w was set equal to zero, [cf. (4.4)]. 
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Compressor Exit Guide Vane (EGV) 


The EGV cascade consists of highly cambered, modified NACA 0012 airfoils. It has a 
stagger angle, 0, of 15 deg, a gap-chord ratio, G, of 0.6 and operates at a prescribed inlet 
Mach number, M_ and inlet flow angle, O_oo, of 0.3 and 40 deg, respectively. Calculations 
were performed for a purely inviscid flow, and for viscous flows at Reynolds numbers of 10 5 
and 10 6 . Instantaneous transition from laminar to turbulent flow was assumed to occur at one 
percent of the arc distance measured along the blade surfaces from the leading-edge stagnation 
point to the trailing edge on both the suction and pressure surfaces of the blades. A streamline 
H-mesh is depicted in Figure 4.2, where three adjacent EGV blade passages are shown. For the 
purpose of illustration, the mesh shown in this figure has approximately one-half the number 
of axial and circumferential grid lines as were used for the actual calculations. 

Results of the inviscid and IVI calculations are shown in Figure 4.3. The blade and wake, 
pressure coefficient and displacement thickness distributions are shown in Figures 4.3a and 
4.3b, respectively; the blade-surface shear-stress distributions, in Figure 4.3c. The expected 
approach of the viscous solutions to the inviscid solution as Re is increased is evident in 
the pressure coefficient predictions. The rate of growth of the suction-surface displacement 
thickness increases dramatically with increasing chordwise distance, x , as the viscous-layer 
separation point is approached. As shown in Figure 4.3c, suction- surface separation bubbles 
(fu, < 0) exist and span approximately 14 percent of chord for Re = 10 6 and about 24 
percent of chord for Re = 10 5 . The decrease in the extent of the separation bubble as Re is 
increased is consistent with the behavior expected for turbulent flows. Note that the suction- 
surface pressure distributions in Figure 4.4a flatten after the flow separates from the blade, 
but adverse pressure gradients still persist within the separation regions, because the suction 
surface pressures rise to meet those on the pressure surface as the trailing edge is approached. 

The blade-surface pressure coefficient, displacement thickness and shear-stress distribu- 
tions, as predicted by the SFLOW-IVI analysis for Re = 10 6 , are compared in Figure 4.4 with 
results obtained using the Navier-Stokes analysis of Dorney, et al. [DDE92]. This Navier- 
Stokes analysis uses the Baldwin- Lomax turbulence model [BL78], which is very similar to 
the Cebeci-Smith model used in SFLOW-IVI. Good agreement between the results of the two 
analyses, particularly for the shear-stress distributions, has been obtained over most of the 
blade surface. However, the agreement, particularly that between the displacement thickness 
distributions, deteriorates in the vicinity of the trailing edge. This is due to the use of an 
O-rnesh around the blades in the Navier-Stokes analysis, which is not well-suited for predict- 
ing flows over thin or wedge-shaped trailing edge geometries. In such cases the lines of the 
O-mesh become severely skewed in the vicinity of the trailing edge, introducing inaccuracies 
into the numerical solution. Also, the displacement thickness [cf. (4.13)] should be evaluated 
by integrating along lines that are normal to the body surface. However, near the trailing edge 
the skewed “radial” lines of the O-mesh deviate significantly from lines normal to the body 
surface, producing questionable results for 6 in the vicinity of a blade trailing edge. Both the 
Navier-Stokes and the IVI analyses predict separation ( f w < 0) near the trailing edge, and 
give almost identical predictions for the location of the separation point (f w = 0); as indicated 
in Figure 4.4c. 

To test the robustness of the SFLOW-IVI analysis, additional calculations were carried 
out for viscous flows at Re = 10 6 . The inlet Mach number was held at M_ 0 0 = 0.3, but a wide 
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Figure 4.3: Inviscid ( ) and IVI, at Re = 10 s ( ) and Re = 10 6 ( ), 

solutions for the EGV cascade: (a) pressure coefficient; (b) displacement thickness; (c) surface 
shear stress. 

range of inlet flow angles, i.e., 36 deg < fl_oo < 54 deg was considered The transition point 
locations were held fixed at t/t T e = 0.01 for all values of fLoo- This location is the same as 
that reported earlier for the baseline (O-oo = 40 deg) calculation. The results are shown in 
Figure 4.5. Here, the predicted total pressure loss parameter, u = ( p T- * - ^r+oo )l( p T-«, - 
P-oo), where Pt is the total pressure, exit flow angle, fl +00 , and suction-surface separation 
point location, x sep , are plotted as functions of inlet flow angle fl-oo- At fl_ oc = 54 deg, 
the viscous flow is approaching stall, with the separation region spanning approximately 35 
percent of chord. Above 54 deg the IVI calculations did not converge due to a numerical 
instability. This is consistent with the known stability properties of the semi-inverse IVI 
iteration procedure when applied to flows with large-scale separations [Wig8l]. 
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Figure 4.4: Comparison of IVI ( ) and Navier-Stokes ( ) solutions for the EGV 

cascade at Re = 10 6 : (a) pressure coefficient; (b) displacement thickness; (c) surface shear 
stress. 

The total pressure loss parameter and the exit flow angle are plotted versus fl_oo in Fig- 
ures 4.5a and 4.5b, respectively. There is a range of inlet flow angles over which the loss remains 
relatively low, but u increases rapidly as the inlet flow angle is increased above 50 deg. The 
latter behavior corresponds to a significant increase in the extent of the separation region 
with increasing Cl- for Cl-^ > 50 deg, as can be seen from the results for x sep shown in Fig- 
ure 4.5c. A striking similarity exists between the variations in — fl +00 and x aep with. fl_oo, as is 
apparent from the results shown in Figure 4.5. The streamwise growth of the separation bub- 
ble, with increasing fl_oo, is accompanied by an increase in the suction-surface displacement 
thickness in the vicinity of the trailing edge. This produces a thickened displacement body 
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Figure 4.5: SFLOW-IVI predictions for the EGV cascade operating over a range of inlet flow 
angles: (a) loss parameter; (b) exit flow angle; (c) separation point location. 

(i.e., the profile made up of the viscous displacement thickness superimposed on the actual 
blade), thereby reducing the effective camber of the blade, and hence, the blade loading. As 
a direct consequence, there is a reduction in the turning of the flow, i.e., an increase in fl +O0 . 

The predicted streamline patterns indicating the size of the trailing-edge separation bubble 
for O_oo = 36, 45 and 54 deg are shown in Figure 4.6. For 36 < fi_oo < 45 deg, the separation 
bubble grows slowly, whereas a much more rapid growth occurs between 45 and 54 deg [see 
Figure 4.5c]. The “decambering” effect produced by the growth of the separation bubble is 
also indicated by the results in Figure 4.6. The kinks that appear in the streamlines near 
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Figure 4.6: Trailing-edge streamline patterns for the EGV cascade: (a) ft_oo = 36 deg; (b) 
Cl - oo = 45 deg; (c) 0_ oo = 54 deg. 

the trailing edge require some explanation. Since the blade trailing edge is wedge shaped, 
the surface coordinate line formed by the blade surface and reference wake streamline has a 
geometric singularity of “kink” at the trailing edge. This singularity influences the solution 
throughout the trailing-edge region as shown by the streamline plots in Figure 4.6. Because 
this singular behavior is highly localized, its effect on the solution for the overall flow field 
appears to be negligible. 

High-Speed Compressor Cascade (HSC) 

The HSC cascade or Tenth Standard Configuration, consists of modified NACA 5506 
airfoils. It has a blade spacing of unity and a stagger angle of 45 deg. We consider a high- 
subsonic inlet operating condition, i.e., M - <*, = 0.7 and Q_oo = 55 deg, and viscous flow's at 
Reynolds numbers of 10 s and 10 6 . Instantaneous transition is assumed to occur at ten and at 
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Figure 4.7: Streamline H-mesh for the HSC cascade. 

one percent of the surface arc length for the flows at Re = 10 s and 10 6 , respectively, on both 
the suction and pressure surfaces of each blade. The IVI analysis for the flow at Re = 10 5 was 
found to be sensitive to the specified transition location. In particular, the iterative solutions 
would not converge, if transition was specified to occur at one percent of arc length, whereas 
if transition was assumed to occur further downstream, the iterations converged. The cascade 
along with a streamline H-mesh, which has, for the sake of clarity, a lower grid point density 
than that used for the actual calculations, are shown in Figure 4.7. The results of the inviscid 
and viscous calculations are presented in Figures 4.8 and 4.9. 

The predicted pressure and displacement thickness distributions along the blade surfaces 
and the wake are shown in Figures 4.8a and 4.8b, respectively, for pure inviscid flow and for 
the viscous flows at Re = 10 5 and Re = 10 6 . The behavior of both of these quantities is similar 
to that observed for the EGV. The surface shear-stress distributions, shown in Figure 4.8c, 


76 



0 0.25 0.50 0.75 1.0 


x 


Figure 4.8: Inviscid ( ) and IVI, at Re = 10 5 ( ) and Re — 10 6 ( ), 

solutions for the 10th Standard Cascade: (a) pressure coefficient; (b) displacement thickness; 
(c) surface shear stress. 

indicate that the streamwise extents of the suction-surface separation bubbles are smaller than 
those predicted for the EGV cascade, decreasing from approximately 20 percent to about 8 
percent of chord as the Reynolds number is increased from 10 s to 10 6 . The kinks in the shear- 
stress distributions for the Re — 10^ case are associated with the instantaneous transition 
that occurs at t/tx£ = 0.1 on both the suction and pressure surfaces of each blade. 

The surface pressure coefficient, displacement thickness and shear-stress distributions de- 
termined for the flow at Re = 10 6 using SFLOW-IVI are compared with those obtained using 
the Navier-Stokes analysis of [DDE92] in Figure 4.9. The agreement is excellent except in the 
immediate vicinity of the trailing edge. Again, the differences between the two solutions are 
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Figure 4.9: Comparison of IVI ( ) and Navier-Stokes ( ) solutions for the 10th 

Standard Cascade at Re = 10 6 : (a) pressure coefficient; (b) displacement thickness; (c) 
surface shear stress. 

attributed to the use of different meshes (H- and 0-) around the blade. The two analyses give 
almost identical predictions for the location of the separation point. 

Turbine Cascade 

The turbine cascade is a modified version of the Fourth Standard Configuration described 
in the study of Fransson and Suter [FS83]. The blade geometry is shown in Figure 4.10 and was 
obtained by modifying the original blunt trailing-edge geometry to produce a wedge-shaped 
trailing edge, while retaining the original chord length, as discussed in [HV93]. As for the 
compressor solutions discussed above, the mesh shown in the figure has fewer grid lines than 
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Figure 4.10: The turbine cascade and streamline H-mesh. 


were used in the actual calculation. The streamline mesh employed for the turbine calculation 
was obtained from an inviscid solution calculated on a composite mesh (see, e.g., Figure 3.9). 
The composite mesh solution capability [UV91] is available in both the SFLOW and LINFLO 
analyses. 

The blade spacing and stagger angle for the turbine cascade are 0.76 and 56.6 deg, re- 
spectively, and the inlet Mach number M_oo and flow angle O_oo are 0.205 and 45 deg. The 
value of M-oo has been adjusted from the experimentally measured value of 0.190 to improve 
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Figure 4.11: Results for turbine cascade: (a) comparison of predicted and measured Mach 

number distributions: ( ) IVI; ( ) inviscid; symbols: experiment; (b) predicted 

displacement thickness distribution; (c) predicted surface shear-stress distribution. 

the agreement with the measured pressure distribution. The calculation was carried out at 
a Reynolds number of 5 x 10 5 with instantaneous transition occurring at 10 percent of the 
surface axe length downstream from the leading-edge stagnation point along both the suction 
and pressure surfaces of the blades. A converged solution could not be obtained for the turbine 
if the location of transition was specified to be too close to the leading edge. This is consistent 
with the behavior observed for the HSC cascade for the flow at Re = 10 5 . 

The IVI solution was obtained in 12 global inviscid/ viscid iterations. The viscous subit- 
eration procedure described in § 4.2 was very effective for this case, reducing the CPU time 
needed to converge the calculation from 1371 seconds without subiteration (requiring 115 
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global iterations) to 224 seconds (in 12 global iterations), using four viscous subiterations 
(i.e., Nv = 4) during each global IVI iteration. 

The computed and measured blade surface Mach number distributions are shown in Fig- 
ure 4.11a. Viscous effects produce a nearly uniform decrease in the suction surface Mach 
number distribution aft of x 0.4, whereas the pressure surface Mach number distribu- 
tion is almost unaffected. The agreement between the IVI solution and the experimental 
data is reasonable; the disagreement in the trailing-edge region can be attributed to the geo- 
metric modification mentioned above. It is difficult to draw definitive conclusions regarding 
the compaxison between the predictions and the data because the solution for this case is 
particularly sensitive to the inviscid mesh used. The predicted displacement thickess and sur- 
face skin-friction coefficient distributions are shown in Figures 4.11b and 4.11c, respectively. 
No separation was predicted for the flow conditions considered, however, the suction-surface 
viscous-layer is close to separation at the trailing edge. 

Timing Study and Convergence Behavior 

Because the development of an efficient analysis has been a major objective of this analyti- 
cal investigation, a timing study was conducted for the three cascade configurations examined 
herein. This study provides both a measure of the computational effort currently required to 
obtain solutions using SFLOW— IVI, and benchmarks against which future efforts to improve 
efficiency can be compared. The results are summarized in Table 1. In addition to the CPU 
time tc, the relaxation factor w, and the number of global iterations N G required to converge 
solutions using a tolerance level e of 0.001 are given in Table 1. The execution times were 
determined using a nearly optimal value of w, as determined by a trial and error procedure. 

The calculations were carried out on an HP-Apollo 720 workstation where SFLOW-IVI 
has been compiled using an optimizing preprocessor. No attempt has been made to “tune” 
the code to take advantage of special features of the optimizer. The times given in Table 1 
are CPU times for the portion of the calculation associated with the IVI iteration loop. Any 
overhead associated with the initialization of the data structure, the generation of the mesh 
and the calculation of the initial inviscid solution is not included. However, this overhead 
amounts to a small percentage of the overall CPU time required by the SFLOW-IVI analysis. 
Note that each of the solutions was obtained in less than five minutes. Recent calculations on 
an IBM RS /6000 3CT workstation indicate a factor of four (4) reduction in the computing 
times reported in Table 1 . 

It is difficult to make direct comparisons with Navier-Stokes CPU times since these can 
vary considerably, even by orders of magnitude, depending on numerous factors, including 
whether the code is a research or a design code, the number of grid points, the grid stretching, 
the convergence tolerance, and so on. An estimate based on a Navier-Stokes analysis that is 
currently used in design indicates that the present IVI analysis requires one to two orders of 
magnitude less CPU time to produce similar results. 

The convergence behavior of two parameters of interest to compressor blade designers was 
examined to determine if a different measure of convergence than that given by equation (4.28) 
would be more appropriate for engineering applications. For the two compressor cascades, 
the total pressure loss parameter u5 and exit flow angle fl+oo were monitored during the 
IVI iteration procedure. We have found that the values of uJ and ft +0O could be considered 
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Table 1. Summary of SFLOW-IVI CPU times, t c , for different cascade configurations. 


Configuration 

LJ 

N g 

tc (secs) 

EGV, Re = 10 6 

1.20 

24 

197 

EGV, Re = 10 5 

0.85 

38 

277 

HSC, Te = 10 6 

1.20 

27 

203 

HSC, Re = 10 s 

0.80 

40 

296 

Turbine 

0.55 

12 

224 




Figure 4.12: Convergence history for the EGV cascade at Re = 10 6 and = 40 deg: (a) 
total pressure loss parameter; (b) exit flow angle. 


converged at a significantly lower iteration count than was needed to satisfy the convergence 
criterion (e = 0.001); typically about one-third fewer iterations than are shown in Table 4.3 
are needed. Thus, even greater efficiency could be achieved in many cases by measuring 
convergence by the degree to which the parameters of interest have approached “asymptotic” 
values. This is demonstrated by the results presented in Figure 4.12, which show the behavior 
of & and f2 +oe , respectively, as functions of the iteration count for the EGV cascade operating 
at Re = 10 6 and U_ oc = 40 deg. This behavior is typical of that observed for all of the 
cases studied herein. The solution for the case illustrated in Figure 4.12 converged to within 
e = 0.001 in 24 iterations while the asymptotic value, indicated by the dashed horizontal line, 
was determined by converging the solution to c = 0.0001, for which 41 iterations were needed. 
For engineering purposes, this solution could be considered to be converged after about 15 
iterations, for which tc ~ 120 seconds. 
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Discussion 


Existing nonlinear inviscid and inverse viscous-layer analyses have been extended and cou- 
pled to provide a strong inviscid/viscid interaction analysis (SFLOW-IVI) for two-dimensional, 
steady, subsonic cascade flows. This IVI solution procedure can be used to predict the ef- 
fects of local strong interactions, including trailing-edge/near-wake interactions and small- 
to moderate-scale viscous-layer separations, on cascade performance. The present analysis is 
restricted to subsonic flows, but it can be extended to treat transonic flows in the future. 

The SFLOW-IVI analysis has proven to be both efficient and robust. Converged solutions 
for each of the baseline configurations examined were obtained in less than five CPU minutes 
on an HP- Apollo 720 Workstation. Even lower CPU times could be obtained by basing conver- 
gence on the global quantities of interest to an engine designer. Robustness was demonstrated 
via application to a wide range of inlet conditions, including cases with large-scale separation, 
spanning up to 35 percent of blade chord. 

A number of issues still need to be addressed in order to improve the accuracy of the 
SFLOW-IVI analysis and to expand its range of applicability. Among them are the inclu- 
sion of quasi- three-dimensional (i.e., streamtube contraction and radius change) effects, the 
incorporation of better models for transition and turbulence, and the addition of a procedure 
for updating the location of the wake streamline during the global iteration process. In addi- 
tion, the overall utility of this SFLOW-IVI analysis for design-system applications needs to 
be explored through further testing and validation. Finally, as this effort continues, a steady 
transonic capability should be developed and the focus should turn increasingly towards the 
development of a strong inviscid/viscid interaction capability for subsonic and transonic un- 
steady flows. 

4.4 Numerical Examples: Effects of Strong Steady Inviscid/Viscid Interaction 
on Unsteady Response 

We have applied the SFLOW-IVI and LINFLO analyses in an effort to estimate the effects 
of steady viscous displacement on the unsteady aerodynamic response of a cascade undergoing 
prescribed blade motions. The calculation procedure is as follows. A strong inviscid/viscid 
interaction solution is first determined by applying the SFLOW-IVI analysis using an H-type 
mesh for the flow in the outer or inviscid region. The blade and wake displacement-thickness 
distributions resulting from this calculation are then used in the surface conditions (4.2) and 
(4.4), and the SFLOW analysis is applied on a composite-mesh to determine an accurate 
solution for the inviscid component of the high Reynolds number flow. This composite-mesh 
SFLOW solution is then used to provide the steady background flow information needed to 
determine a linearized inviscid unsteady solution, also on a composite mesh, using LINFLO. 

In this approach, the SFLOW inviscid solution accounts for the effects of steady viscous- 
displacement via the imposition of the appropriate blade and wake boundary conditions. The 
linearized inviscid flow is determined as a solution of the field equation (3.9), with v, = 0, and 
the blade and wake surface conditions (3.11) and (3.12). Thus, steady viscous displacement 
effects are incorporated into the linearized inviscid analysis through the steady potential $, 
which is the potential for the outer inviscid flow of a strong inviscid/viscid interaction calcu- 
lation. Unsteady viscous-displacement effects are not taken into account by this procedure. 
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Figure 4.13: Surface Mach-number distributions for inviscid ( ) and viscous, at Re = 10 s 

( ) and Re — 10 6 ( ), flows through the 10th Standard Cascade operating at 

M_ «, = 0.70, Q_oo = 55 deg. 

The foregoing approach has been applied to predict steady and unsteady flows through 
the 10th Standard or HSC cascade operating at an inlet Mach number of 0.7 and inlet flow 
angle of 55 deg. The blade-surface, Mach number distributions, as determined using SFLOW 
on a composite mesh, for pure inviscid flow and for viscous flows are Re = 10 6 and Re = 10 s 
are shown in Figure 4.13. The C p , 6 and f w distributions for the^ flows, as determined by 
an H-mesh IVI solution, are shown in Figure 4.8. The predict, exit Mach number, exit 
flow angle, and lift force for three steady flows are 0.447, 40.2 deg, and 0.348 for Re — > oo; 
0.470, 41.5 deg, and 0.321 for Re = 10 6 ; and 0.488, 42.2 deg, and 0.303 for Re = 10 s . The 
viscous flows at Re = 10 6 and Re = 10 s separate from the suction surface at = 0.927 and 
£sep = 0.808, respectively. The trailing-edge streamlines for the flow at Re = 10 5 are depicted 
in Figure 4.14. 

Unsteady response predictions, i.e., predictions for the global and local works per cycle 
[Ver89a, Ver93], are shown in Figures 4.15-4.17 for blades undergoing single-degree-of-freedom 
torsional, with a = (1,0), and bending, with h y = (1,0), vibrations. The torsional vibrations 
occur about the blade midchords. The global work per cycle, Wc , is the work done by the 
airstream on a given blade over one cycle of its motion. Therefore, a prescribed blade motion 
is stable, neutrally stable, or unstable according to whether the (global) work per cycle is 
less than, equal to, or greater than zero, respectively. The local work per cycle or pressure- 
displacement function, wc , describes the distribution of the work per cycle over a blade 
surface. For subsonic flows in which the blades are undergoing small-amplitude rigid-body 
motions these quantities are determined from the relations 

W c = (j/) w c (TB)dfB , (4.29) 
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Figure 4.14: Trailing-edge streamlines for 10th Standard Cascade operating at A/_oo — 0.70, 
Q-oo = 55 deg, and Re = 10 s . 

and 

wc = — 7r/m{PB<*RjB * r B + PbR-b ' n s} + • • • • (4.30) 

Here R B = h + a x R P , Rp is the distance from the moving pitching axis to a point on the 
moving blade surface, Pq is the steady pressure at a point on the mean blade surface, p& is 
the complex amplitude of the unsteady pressure at the corresponding point on the moving 
blade surface, and the superscript * denotes the complex conjugate. 

The local work-per-cycle or pressure displacement function distributions for unit frequency 
(uj = 1) in- (a = 0) and out-of-phase (a = 180 deg) torsional blade motions about midchord 
are shown in Figure 4.15. For the in-phase torsional motion there are small differences between 
the inviscid and the viscous solutions, and no clear trend in the behavior of wc with increasing 
Reynolds number. The normalized global works per cycle |a| -2 Wc for the three flows at a = 0 
are -0.897, -0.777 and -0.784 for Re = 10 5 , Re = 10 6 and Re -» oo, respectively. For the 
out-of-phase torsional motion the local work per cycle decreases with Reynolds number on 
both the blade pressure and suction surfaces over approximately the first thirty (30) percent 
of blade chord, indicating that the stability margin for the out-of-phase torsional blade motion 
increases with decreasing Reynolds number. Here, |a| -2 Wc = —1.732, —1.578, and —1.252 
for the flows at Re = 10 5 , Re = 10 6 and Re-* oo, respectively. 

Similar results for unit frequency, in- and out-of-phase, bending vibrations are shown in 
Figure 4.16. For the in-phase bending motion there are changes in the local work per cycle 
with Reynolds number on both the blade suction and pressure surfaces in the vicinity of the 
leading edge, and along the suction surface, where wc increases with decreasing Reynolds 
number over an interval extending from 15% to 60% of blade chord. The local work per cycle 
distributions for the two viscous flows at Re = 10 5 and Re = 10 6 are almost identical. The 
normalized global works per cycle, |/ij,| -2 Vkc, for these flows are —3.488 for Re —* oo, —3.270 
for Re = 10 6 and -3.261 for Re = 10 5 . Thus, for the in-phase bending vibration, the stability 
margin decreases slightly with decreasing Reynolds number. For the out-of-phase bending 
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vibration, viscous displacement effects cause a small increase in wc on the pressure surface 
over the first twenty (20) percent of blade chord and a large increase over most of the suction 
surface. As a result of the latter, viscous effects, which increase with decreasing Re , are 
strongly destabilizing for the out-of-phase bending motion, but this type of blade motion has 
a substantial stability margin. The global works per cycle, \h y \~ 2 Wc , for the flows analyzed 
in Figure 4.16b are: = —10.14 for Re —> oo, —8.812 for Re = 10 6 and —8.667 for Re = 10 s . 

Global work per cycle, Wc, predictions for flc vs at Re = 10 5 , Re = 10 6 and Re — » oo 
through the 10th Standard Cascade are shown in Figure 4.17. Here the blades are under- 
going pure torsional and pure bending vibrations at unit frequency, and the response pre- 
dictions are given over an entire 360 deg range of interblade phase angle. For the inviscid 
flows, the interblade phase-angles at which acoustic resonances occur are crZ^ = —26.9 deg, 
<r±oo = H7.1 deg, = —31.8 and cr +00 = 59.8 deg, where the subscripts refer to the far 
upstream (— oo) and far downstream (-(-oo) regions and the superscripts refer to the direction 
of wave propagation. The viscous flows have the same upstream resonant conditions, but, 
since the steady exit flows vary with Reynolds number, slightly different downstream condi- 
tions. The results in Figure 4.17a indicate that steady viscous effects have a negligible impact 
on the torsional stability margins for superresonant blade motions, i.e., the motions ocurring 
between the lowest and highest resonant interblade phase angles. Also, viscous effects are sta- 
bilizing for the subresonant torsional blade motions. Viscous effects are destabilizing for both 
superresonant and subresonant bending vibrations (Figure 4.17b), but such motions tend to 
have substantial stability margins. 

We must emphasize that the foregoing parametric studies are very limited and much 
more detailed studies will have to be conducted in order to understand the effects of viscous 
displacement on blade flutter characteristics. Also, these results account, only partially, for 
steady viscous displacement effects and neglect unsteady viscous effects. In particular, the 
contributions from the right-hand-sides of the linearized surface boundary conditions (2.56)- 
(2.58) have not been taken into account. That said, if the analysis described above is found to 
provide useful approximations to the effects of viscous displacement on the unsteady pressure 
responses of blade rows, then convenient and very efficient viscous unsteady aerodynamic 
models could be applied to predict a wide range of blade-row aeroelastic and aeroacoustic 
phenomena. 
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Figure 4.15: Effect of steady viscous displacement on pressure- displacement function response 
for subsonic flow (M_ «> = 0.7, fl_oo = 55 deg) through the 10th Standard Cascade undergoing 
torsional blade vibrations about midchord at u> = 1: (a) in-phase (<x = 0 deg) blade motion; 
(b) out-of-phase (cr = 180 deg) blade motion. 
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Figure 4.16: Effect of steady viscous displacement on pressure-displacement function response 
for subsonic flow (Af_oo = 0.7, fl_oo = 55 deg) through the 10th Standard Cascade undergoing 
bending blade vibrations at u = 1: (a) in-phase (a = 0 deg) blade motion; (b) out-of-phase 
(<r = 180 deg) blade motion. 
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Figure 4.17: Effect of steady viscous displacement on work per cycle responses for subsonic 
flow (M_oo = 0.7, Cl - oo = 55 deg) through the 10th Standard Cascade undergoing prescribed 
blade vibrations at unit frequency (u> = 1): (a) torsional vibrations about midchord (b) 
bending vibrations. 
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5. The Unsteady Viscous Layer Analysis: UNSVIS 


The viscous layer analysis, UNSVIS, applies to the unsteady compressible flows in blade 
boundary layers and wakes. It was originally developed [DBJ + 89, PVK91] to determine the ef- 
fects of unsteadiness on turbine blade heat transfer. In this analysis, the equations that govern 
the flows in viscous layers are transformed using an unsteady turbulent generalization of the 
Levy- Lees transformation. The transformed equations are then solved using a finite difference 
technique in which the solution proceeds by marching in time and in the streamwise direction. 
Both laminar and turbulent flows can be studied, the latter using algebraic turbulence and 
transition models. Laminar solutions for a flat plate have been shown to approach classical 
asymptotic results for both high and low-frequency unsteady motions; turbulent flat-plate re- 
sults, to be in qualitative agreement with earlier predictions and experimental measurements 
[PVK91]. The numerical technique has also been applied to the stator and rotor of a low- 
speed turbine stage to determine unsteady effects on surface heating [DBJ + 89]. The results 
compare reasonably well with measured heat transfer data and indicate that nonlinear effects 
have minimal impact on the mean and unsteady components of the flow. 

Under the present contract, the UNSVIS analysis has been extended so that the unsteady 
viscous effects in the vicinity of leading-edge stagnation points and in blade wakes can be 
predicted. The stagnation point analysis also provides the “initial” upstream information 
needed to advance or march a numerical viscous-layer calculation downstream along the blade 
and wake surfaces. In addition, a wake analysis, used previously in UNSVIS, has been extended 
so that the changes or jumps in the inviscid velocity that occur across vortex-sheets unsteady 
wakes could be accommodated. 

The UNSVIS analysis is described in detail below and demonstrated via application to 
two cascade configurations. In particular, this nonlinear unsteady viscous-layer analysis has 
been applied to study the viscous-layer responses of an unstaggered flat-plate cascade and a 
turbine cascade to external pressure excitations. The flat-plate results indicate that viscous- 
layers respond linearly to upstream pressure excitations over a broad range of excitation 
amplitudes and temporal frequencies, but the responses to downstream pressure excitations 
at high amplitude and/or high temporal frequency have significant nonlinear content. The 
coupled LINFLO/UNSVIS analysis has been applied to a turbine cascade subjected to a pres- 
sure excitation from upstream to demonstrate the current, weak, inviscid/viscid interaction 
solution capability for unsteady flows on a realistic cascade configuration. 

5.1 Governing Equations 

For flows through blade rows, occurring at high Reynolds number, the viscous region con- 
sists of thin boundary layers that lie along the upper and lower surfaces of each blade and 
thin wakes which extend downstream from the blade trailing edges, as indicated in Figure 2.3. 
The equations that govern such flows have been derived in § 2.3 with respect to moving, 
reference, shear-layer surfaces S m , each of which is contained entirely within a viscous layer. 
These surfaces are taken to coincide with the suction and pressure surfaces of the blades and 
to lie entirely within the viscous wake. The curvatures of a reference shear layer surface are 
assumed to be 0( 1) and, at this point, we assume that the amplitude, |72j, of the unsteady 
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blade motion and the local time derivative d/dt\ x are of 0(1). Subsequently, we will restrict 
our consideration to blades undergoing small (i.e., 0(e)) amplitude motions. The exact loca- 
tions of the wake shear-layer surfaces (lines) are unknown a priori ; however, to within lowest 
order, wake boundary conditions can be referred to any arbitrary surface emanating from the 
trailing edge and lying within the actual viscous wake [Vel80]. For steady flow calculations, 
the reference wake surfaces are taken to be the aft stagnation stream surfaces, as determined 
from an inviscid solution. 

It is convenient to use forms of the viscous-layer equations that contain the relative fluid 
velocity and relative total enthalpy, i.e., 

U = V-£- and H r = H + U 2 /2, (5.1) 


respectively, as dependent flow variables. Here 7?.(x,f) = x(x,t) - x is prescribed such that 
72.(x, t) = 7?-b(x, t) for x € B. Thus, the two-dimensional flows in the thin boundary layers 
that lie along the upper and lower surfaces of each blade and in the thin wake that extends 
downstream from the blade trailing edge are governed approximately by the following system 
of field equations, cf. (2.28)— (2.30) : 


d£ 

dt 


and 


BHr 

dt 


dUr 

dt 


+ (U-V)#J 
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= Si(: 

M) , 


(5.2) 

’ + s- (Ra) ‘ 
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= S 2 (x,t) , 

(5.3) 


(? - (jv)-*; 

f )<' 


- M tK 

= 5 3 (x,<) . 










(5.4) 


Here, r and n measure distances along and normal to the moving shear layer surface, and 
U r and U n are the fluid velocity components in the positive r- and n-directions, respectively, 
relative to the moving shear-layer surface. A distinction should be made between the inde- 
pendent variables that describe the flows in the upper and lower surface boundary layers and 
in the wakes, e.g., by attaching subscripts to r and n. However, as a convenience, we are 
neglecting to make this distinction explicitly. 

The source terms on the right-hand-sides of (5.2)-(5.4) are given by 


Si = — V x • (pit) , S 2 = —p 


Dn 


Dt 


T + (n-V x )U r 


and 


S 3 = -p 


D_ 
L Dt 


(t/ T R T + K;/2) + (-R-V„)ft t ] , 


(5.5) 


where 1Z, = d'R,/dt\ x . If the blade motions are of small-amplitude, i.e., on the order of the 
displacement thickness, I, then to within the order of the viscous-layer approximation, the 
right-hand sides of (5.2)-(5.4) can be regarded as negligible. 
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The symbols p, U r , U n , P, Hr, p and k in (5.2)-(5.4) refer to ensemble (or Reynolds) 
averaged values of the fluid dynamic variables; v' r , v' n and h' T axe the values associated with 
random turbulent fluctuations; and the overbar indicates a turbulent correlation, which must 
be determined empirically. As a consequence of the high Re and surface curvature assump- 
tions, the pressure in the thin viscous layers is a function only of r and t, and Hr « T + U^/2. 
Also, the pressure and the flow properties at the edge of each viscous layer are equal to the in- 
viscid values of these variables at the reference shear layer surface. If the interactions between 
the flows in the viscous layer and the external inviscid stream axe weak, then the pressure 
within a viscous layer and the flow properties at its edge are determined by an inviscid solution 
for these flow variables. The latter is determined subject to surface conditions that account 
for the effects of viscous displacement. 

In addition to the foregoing field equations, the equation of state for a perfect gas, i.e., 
(2.11), empirical laws relating the molecular viscosity and the thermal conductivity to the 
temperature, e.g., (2.10), and equations relating the turbulent correlations v' T v' n — u' T u' n and 
h'j v ' n = h'Ru' n to the ense mble-averaged flow quantities, axe also required. The turbulent 
correlations v' r v' n and h' T v' n are related to gradients of the mean-flow variables, using Prandtl’s 
mixing-length hypothesis, cf., (2.32). 

Initial and Boundary Conditions 

The foregoing system of field equations, (5.2)-(5.4) is parabolic in time and in the stream- 
wise or r-direction. Therefore, the streamwise component of the relative velocity, U T , and the 
relative total enthalpy, Hr, must be known for all time at some upstream streamwise loca- 
tion, and these variables, along with the normal velocity, U n , must be known throughout the 
solution domain at some initial time. Also, conditions on the fluid properties at the edge(s) 
of the viscous layer, i.e., 

Ur - Ur,e(T,t) and Hr -+ Hr^t,!) for n -» ±oo , (5.6) 

where the limits +00 and —00 refer to the edges of the upper (+) and lower (— ) surface 
boundary layers and the upper and lower edges of a wake, a no-slip condition and either a 
prescribed temperature or heat flux condition at a solid blade surface, i.e., 

U = 0 and Hr — HR fW (r,t) or ^ — —k 1 (Re-Pr)~ 1 Q w for n = 0 , r < pte ; 

( 5 - 7 ) 

and a condition on the fluid velocity normal to a reference wake streamline, i.e., 

U • n = 0 for n = 0, r> tte , (5.8) 

must be enforced. Here the subscripts w and e denote the values of the fluid properties at a 
solid wall and at the edge of the viscous layer, respectively, and the subscript TE refers to the 
airfoil trailing-edge point. The relative fluid velocity, U T , e , and relative total enthalpy, H R>e , at 
the edges of the viscous layers are determined by the inviscid solution along the blade surfaces 
and the reference wake surfaces. 
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Turbulence and Transition Models 


The models used here and in [PVK91] to simulate the effects of turbulence and tran- 
sition on the flow in the viscous layer are the algebraic eddy-viscosity model proposed by 
Cebeci and Smith [CS74], the transition length correlation model proposed by Dhawan and 
Narashima [DN58], and the wake turbulence model proposed by Chang et al [CBCW86]. 
Also, since flows in turbomachines are known to be characterized by high freestream turbu- 
lence levels, a simple modification, developed by Yuhas [Yuh81], has been incorporated into 
the turbulence model to account for the effects of freestream turbulence on the viscous layer. 
These models are easy to implement, and are known to be reasonably accurate for steady 
flows with mild pressure gradients; however, they have only been developed for steady flows. 
We have modified them here for application to unsteady flows simply by replacing steady flow 
variables by their time-dependent counterparts. Thus, the ability of the models given below 
to accurately represent turbulence and transition in unsteady flows is not known; therefore, 
any resulting unsteady flow predictions must be interpreted with some caution. 

The Cebeci-Smith algebraic model divides a solid-surface boundary layer into inner and 
outer regions, where I = I, and I = t 0 , respectively. The inner model is applied from the wall 
out to the point at which t, = t 0 ; the outer model, from this point to the edge of the boundary 
layer. For unsteady flows over moving blades, we set the eddy viscosity in the inner region 
according to 

1 Qjj I 

ti = 7 Tr ^(0.41n) 2 [l - exp(-n/A T )] 2 Re 


dn 


(5.9) 
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(5.10) 


and 7 Xr is a longitudinal intermittency factor which models transitional flow. 
The eddy viscosity in the outer region is given by 
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The parameter 7r has been defined to account for low momentum thickness Reynolds number. 
The Clauser constant, xo, is usually set equal to 0.0168, but following Yuhas [Yuh81] we set 

Xo = 0.0168(1.0 + 18.4(^P)t» + 99.6( ; ^V)_ 0O ] , (5.13) 


where (v' • is the mean-square of the inflow freestream turbulence velocity, to account 

for the effects of freestream turbulence on the development of the turbulent boundary layer. 

Blade boundary layers can contain transitional regions of significant extent. Therefore, 
a transition model is needed to predict such flows. Several models have been proposed for 
steady transitional flows, based upon correlations with experimental data. One of these is the 
Dhawan-Narashima forced transition model [DN58] in which the intermittency factor in the 
transition region is given by 


iTr = 1 - eX P 


-4.6513 


(— y 

\T 2 - T X ) 


, 7i < T < r 2 , 


(5.14) 


Here t x and r 2 are the streamwise locations at the beginning and end of the transition re- 
gion, respectively. These locations can be specified as functions of time to model unsteady 
flows. The intermittency factor is set equal to 0 and 1 in laminar and fully turbulent regions, 
respectively. Therefore if r 2 is set equal to t x the transition is regarded as instantaneous. 

The eddy viscosity in the wake is based on the model developed by Chang et al [CBCW86]. 
At each streamwise station in the wake equation (5.11) is evaluated twice: once for — oo < 
n < rimin, where n m in is the location at which the streamwise velocity reaches a minimum 
value, and once for n^n < n < oo. The maximum of these two values, 6^, is then used to 
set the eddy viscosity, eyv, i.e., we set, 


CyV — Cmax ^mtx ^max,TE] CXp 


T ~ TTE 

20$tE - 


(5.15) 


In this equation r — tte is the distance measured along the reference wake surface from 
the trailing edge point, and <*>te is the sum of the upper and lower surface boundary- layer 
thicknesses at the trailing edge. The boundary-layer thickness is defined to be the normal 
distance from a blade surface to the point at which V T /V T<e = 0.995. 


5.2 Transformed Equations 

For laminar flows, the Levy- Lees transformation [see Blottner [Blo70]] provides indepen- 
dent variables that effectively capture the growth of the viscous layer with increasing stream- 
wise distance. In addition, the transformed equations reduce to similarity equations at an 
airfoil leading edge. These features facilitate the determination of a numerical solution. Thus, 
an extension of the Levy-Lees transformation, in which the laminar edge viscosity is replaced 
by an effective turbulent viscosity [VWV82], is applied herein to accommodate unsteady tur- 
bulent flows. Note that the dependent variables defined below differ from those used in § 4 
for steady flows. Also, in the present implementation, we do not restrict our consideration to 
adiabatic flows at unit Prandtl number. 
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We define new independent variables £ and rj according to 


( = f Qdr 

Jo 


and 77 = 


j7 T ,e\/Re 


I pdn , 
Jo 


(5.16) 


where Q = p e U T , e p e [ 1 + (e//2 ) e ] = pU T , e Pe«,e and e and ji axe the turbulent eddy and molecular 
viscosities in the steady background flow. In [PVK91] the temporal mean values of the fluid 
properties, p e , U T , e , and /x e , at the edge of the viscous layer were used to define Q. However, 
in the present discussion, we regard the unsteady inviscid flow as a small perturbation of a 
nonlinear mean flow. Therefore, we can replace the temporal mean values of these variables 
by their values in the steady background flow to avoid introducing additional nomenclature. 
The coordinate £ in (5.16) is a function of r alone; the coordinate 77 , by contrast, is a function 
of r, n, and t. 

We also introduce the new dependent variables 


T = 
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2£/_l_d77 
Q Ur.e 


dt] 



and 



(5.17) 


where Q = p e U T ,eP e [ 1 + (e//i)e]. After substituting Eqs. (5.16) and (5.17) into the viscous-layer 
field equations and performing some algebra, we find that 
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(5.20) 


95 



where the following parameters have been introduced: 


Qti = 



«3 
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a 4 
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2 i{QU T , e H Rfi ) 


dt 
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Ul,T ,-' a , 
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2{Q(Qa 1 U T , e )- 1 ^ , 

(5.21) 

ffil 1 + lPr/(^Pr T )] 

^ e A e [l + (e/£) e ]fh- 


Special attention must be exercised in applying the foregoing transformation to unsteady 
wakes. In the present analysis all quantities on the upper (i.e., n > 0) and lower (n < 0) sides 
of a wake are referenced to their respective edge conditions, i.e., we set T = for n > 0 

and T = U T /U~ e for n < 0. In general, the inviscid streamwise relative velocity, U T , e , and 
relative total enthalpy, Hr te , will be discontinuous across unsteady wakes. As a consequence, 
the dependent variables T and Q axe discontinuous across reference (n = 0) wake streamlines. 
These discontinuities must be taken into account in develop ng a numerical solution procedure, 
to ensure that the physical variables ( U T , U n and Hr) ai continuous across viscous wakes. 
Also, the upper- and lower-surface values of the independent variable £ will generally differ 
at a blade trailing edge. In the present effort, the upper-surface ^-variable is used to continue 
the viscous calculation into the wake, i.e., we set 


£vv — £te + / Q + dr , 
J tte 

However, two different ^-coordinates 


T > T te 


(5.22) 



t > tte , 


n < 0 , 


(5.23) 


axe used for the wake calculation. Thus, at the trailing edge, the lower-surface boundary-layer 
solution must be interpolated onto the wake rj - mesh before continuing the viscous calculation 
into the wake. 


Boundary and Initial Conditions on T and Q 

In addition to satisfying the field equations (5.18) — (5.20), the dependent variables T and 
Q must be prescribed functions of £ and 77 at some initial time, and prescribed functions of rj 
and t at initial or upstream streamwise locations on the upper and lower surfaces of each blade. 
Also, in terms of the transformed variables, the following conditions must be imposed at the 
edge(s) of the viscous layer, on the blade surface, and along the reference wake streamline, 
cf. (5.6), (5.7) and (5.8): 

J- — * 1 and Q — * 1 for |t/| — » 00 , (5.24) 
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T = 0, V = 0, and Q — G w {£,t) = Hr^/Hr^ or 


(5.25) 


and 


(dg/drj) w = 


-(2tRe)±%PrQ w {t,t) 

P e Ur,eH R,e 


for T ) = 0, £ < £tE , 


V = 0 for rj = 0, £ > £te , 


(5.26) 


respectively. 

For the calculations reported herein, the condition at the initial time is the viscous-layer 
solution in the absence of unsteady excitation, i.e., the solution for the steady background 
flow. The calculated transients resulting from this approximation have been found to die 
out with increasing time. The upstream profiles, required at each time step of the unsteady 
viscous-layer calculation, are obtained from a similarity analysis of the flow in the stagnation 
region, as described below. 


5.3 Solution Procedure 


Stagnation Region 

To determine solutions for the flows in the viscous layers, we first seek a local, similarity 
solution for the flow in a leading-edge stagnation region. This solution can then be used 
to set the instantaneous upstream profiles, •7 r (£/\f?,t) and that are required to 

initiate a space-marching, numerical, viscous-layer calculation. The overall strategy is to 
develop a similarity solution for a simple flow configuration and then, to adapt this result for 
application to the unsteady cascade problem. The present analysis for unsteady, compressible 
flow is modeled after the incompressible analyses developed independently by Rott [Rot56] and 
Glauert [Gla56]. Their analyses provided an exact solution to the Navier-Stokes equations for 
incompressible flow. An exact result can not be determined for compressible flow, however, 
because it is necessary to neglect the dissipation term in the energy equation to obtain a 
similarity solution. Fortunately, frictional heat dissipation in the stagnation region is of limited 
importance for a wide range of practical flows. 

We consider two-dimensional, compressible flow around a flat plate, which is oriented 
normal to the stream direction and undergoes a harmonic motion at velocity iu>Rr exp(iwf) 
in its own plane. As a mathematical convenience, we use complex-variable representations 
for the unsteady flow properties, but it is to be understood that the real parts of the various 
complex parameters represent the actual physical variables. The flow in the viscous layer 
will be determined in terms of space-fixed Cartesian coordinates, rs and ns. The resulting 
stagnation-region solution can then be expressed in terms of plate-fixed coordinates as a 
prerequisite to its implementation into a complete unsteady viscous-layer calculation. The 
rs and n s axes are directed along and normal to the plate, respectively, and the coordinate 
origin lies at the point at which the dividing inviscid streamline impinges on the plate. The 
fluid velocity components in the rs- and ns-directions are and V ns , respectively. 

The inviscid flow is steady in the space-fixed frame, with velocity and pressure gradient at 
the edge of the viscous layer (i.e., at n s = 0 in the inviscid region) given by 

V e = K s , e e Ts = cr s e Ts , (5.27) 
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and 


(5.28) 


8P dV Ts , t 2 

’ ■* Pe'' r s< e Qts ^ TS 7 "®^ Ts ’ 

The constant c in these equations is determined from the behavior of the inviscid solution in 
the vicinity of the mean stagnation-point location (7$ = 0). The total enthalpy at the edge of 
the viscous layer is constant, i.e., Hr, e = T e + e /2 = //_<*>. 

The flow in the viscous layer is described by the field equations (5.2)-(5.4) or, in this 
case, (2.28)-(2.30), the edge conditions (2.36) and the surface conditions (2.34). The surface 
velocity components are 


K s (rs,ns,t) = iuRrexp(iujt) and V ns = 0 for ns = 0 


(5.29) 


We assume that the flow in the stagnation region is laminar and that the heat generated by vis- 
cous dissipation is negligible, i.e., pd 2 V Ts !dn\ « 0. The streamwise velocity and temperature 
can be expressed in the forms 

K s ( ts,7M) = CT S f(r}) + itx>Rrg(T]) exp(turf) and f(T S ,T),t ) = f e (T S ,t)d(r)) , (5.30) 

where crsf'(r)) and iu>Rrg(r))exp(iujt) are the steady and unsteady components of the fluid 
velocity in the viscous layer and 


,=(^Y7 n w 

V PePe ) J 0 


(5.31) 


Note that dP/dns 0 and V^ s <C V rs ] therefore, in the viscous layer, p t j p ~ T j T e ~ 9 and 

Hr^f+V 2 / 2 . 

After substituting the foregoing results into the field equations that govern the flow in the 
viscous layer and performing the necessary algebra, we obtain an expression for the normal 
component of the fluid velocity, V^ s , and the following set of ordinary differential equations 
for /, g and §. 

, (5.32) 

<7 + f)S~fs' -(!«')' = 0 , (5.33) 

and 

77+/«=°> (5-34) 

where / = pp,/(p e p e ). The conditions on /, g and 0 at the plate surface (77 = 0) and at the 
edge of the viscous layer (77 — ► 00) are 

/(0) = f(0) = 0 , g( 0) = 1 and 0(0) = 0 or 0(0) = Tjf e ; (5.35) 

and 

/'(oo) = 1 , 0(00) = 0 and 0(oo) = 1 . (5.36) 

Note that for steady (g = 0), incompressible (0=1) flow, we recover the classical stagnation- 
region problem studied by Hiemenz. Also, equation (5.34) is identical in form to the energy 
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equation that describes steady, compressible, stagnation-point flow [Whi74]. Finally, the 
present analysis can be readily extended to consider more general unsteady plate motions, 

e-g., 

N 

Vr s = CT S f'(ri) + ^2 iu) n r Tin g n (r})exp(iu n t) . (5.37) 

n=l 

In this case solutions for N uncoupled equations for the g n , n_= 1, 2, 3, . . . , N, are required. 

The set of nonlinear ordinary differential equations for /, 9 and g can be solved using an 
implicit finite-difference technique similar to that described in the next sub-section and used 
herein to obtain solutions to the full unsteady viscous equations. The functions / and 9 must 
be determined simultaneously, but they are independent of g, which can be determined once 
/ is known. 

To provide the upstream profiles required for a complete viscous-layer calculation it is 
convenient to express the stagnation-region solution in terms of body- or plate-fixed coordi- 
nates. Thus, we consider plate-fixed Cartesian coordinates (t, n), where the r- and n - axes are 
parallel to the spaced-fixed ts- and ns-axes, respectively, and the mean position of the r, n 
coordinate origin coincides with the origin of the spaced-fixed Ts,ns-frame. The streamwise 
positions and velocities of a fluid particle in the two coordinate frames are related by 

rs = r + Rr exp(iu;f) (5.38) 

find 

V TS = U T + iujRr exp(zu><) . (5.39) 

After combining (5.38) and (5.39), we find that 

U T = crf'iv) + b Kv) exp (iut) , (5.40) 

where 

b = cFLr{ 1 — iu>/c) and (1 — iu>/c) \(g) = + {iu/c)[g(r]) — 1] . (5-41) 

At the plate surface U T = 0, and at the edge of the viscous layer 

Ur(v -*■ °o) = Ur,e = ct + bexp(iu>t) . (5.42) 

The relative total enthalpy H R = U*/2 + f is given by 

Hr = (cr/') 2 /2 -|- cbTf'\exp(iu)t) + 6 2 A 2 exp(2iurf)/2 + tj) . (5.43) 

The locations, i.e., r = ±|t/|, at which the stagnation-region solution can be applied as an 
initial condition for an unsteady viscous- layer calculation are determined by two criteria. First, 
the stagnation-point motion must be contained within the interval [— |r/|, |tj|] and second, 
to avoid stability problems in the subsequent boundary-layer calculation, those points at 
which flow reversals, i.e., ( dU r /dn) w < 0, occur must also be contained within this interval. 
Reverse flow is associated with the lag in the response of the low-momentum fluid near the 
wall to changes in the velocity at the edge of the viscous layer. The extent of the interval 
over which the stagnation point moves, i.e., [ — |tw|, Ita^I], i s determined from the inviscid 
velocity distribution (5.42) and is given by [-c _1 |6|,c|6|j. The extent of the reverse-flow 
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interval, i.e., [— I r /?L kfllli is determined by the maximum and minimum values of r for which 
(. dUr/drj ) w = 0. We find that 

< c . 1|6| i + iuc-'mirw 

1 — iu)/c 

Thus, |r/| > max(|rA/|, | t/j [ ) . For our application to unsteady cascade flows, we have found 
that setting |r/| = 2c _1 |6| leads to reasonable results in the stagnation region and allows us 
to continue the viscous-layer calculation along the blades and their wakes. 

In applying the foregoing stagnation- region analysis to unsteady cascade flows, we use the 
inviscid velocity distribution along a moving blade surface, i.e., 

V e • r = [V e + v e exp(zwf) + . . . ] • r = [V e + (u e + R) exp(iu;<) + ...]• r, (5.45) 

where V e (r) and v e = u e + R are the nonlinear steady velocity and the complex amplitude 
of the linearized unsteady velocity, respectively, to determine the parameters c and 6 in (5.40) 
and (5.43). In particular, these parameters are determined by matching the analytical velocity 
profile (5.40) to the calculated inviscid velocities at ±|tj|. Thus, with u specified as an input 
quantity, the functions / 6 and g determined by an implicit finite- difference procedure and 
the parameters b and c determined by the inviscid solution in the vicinity of the stagnation 
point, the streamwise velocity and total enthalpy profiles needed to initiate the unsteady 
viscous-layer calculation can be specified. 

Viscous-Layer Numerical Solution Procedure 

The transformed viscous-layer equations (5.18)-(5.20) are solved by applying finite- difference 
approximations for the various partial derivatives that appear, to convert tnis system of dif- 
ferential equations into a system of algebraic equations [VBHA91]. Streamwise and tempo- 
ral derivatives are approximated by first-order accurate backward difference expressions and 
normal derivatives, by second-order accurate central difference expressions The streamwise 
momentum and energy equations contain nonlinear terms, which must be li -ized to ob i 

solutions a 'ach point in space and time. This is accomplished by using i wton itera i 

procedure which the initial guess for the profiles comes from the solution at the pre> s 

streamwise (r) station or time step. The equations are solved by using a local iteration to 
remove the linearization error, and by repeating this iteration until the values of th r *’ow 
variables converge to within a specified tolerance level. The finite-difference approxi; n 

results in a block tri diagonal system of linear algebraic equations at each step, which is d 

using the Thomas block inversion algorithm. 

The analyses for surface boundary layers and wakes are identical, except for the implemen- 
tation of different boundary conditions. In addition to replacing the surface no-slip condition 
with a zero normal velocity condition at the reference wake position, the possibility of jumps 
in the relative tangential velocity and the relative total enthalpy across a viscous wake must be 
taken into account. As discussed previously, the transformed governing equations are written 
in terms of two different sets of variables, each applying on one side of the wake, and scaled to 
the edge conditions for that side. For example, on the upper side of the wake, = Ur!U+ 

and Q+ = Hr/ H% e \ on the lower side, T = and Q~ = H R /H Re . The discretized 
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equations on either side of the reference wake streamline are written in terms of that side’s 
variables. However, at the reference wake line, (77 = 0), variables from both sides of the wake 
are used in the momentum and energy equations, due to the use of central-difference approx- 
imations for the ^-derivatives. For consistency, the equations must be written in terms of a 
single set of variables. This is accomplished by writing the equations in terms of the upper 
surface quantities. Thus, whenever a lower surface variable appears in the equations at 77 = 0, 
it is rewritten in terms of the upper surface edge conditions. For example, the variable T~ is 
written in terms of upper surface variables as T~ • (U~ e /U+ e ) = U~ /U+ e . 

The finite-difference approximation is implemented on a grid that is nonuniform in both 
the £- and 77 - directions, but uniform in time. The grid distribution in the streamwise direction 
is chosen so as to cluster points near the blade leading and trailing edges, with each mesh 
interval being set equal to a constant times the previous one, according to the relation 

A6+i = Kt&Z, , 1 = 1, 2 ,..., L (5.46) 

where / is the streamwise mesh point index (/ = 1 at the blade leading edge and increases with 
distance downstream), K ^ > 1 from the leading edge of the blade to midchord, and K(_ < 1 
from midchord to the trailing edge. The grid is then stretched (i.e., /Q > 1) aft of the trailing 
edge, and generally extends between one and two blade chords downstream of the trailing 
edge. 

A stretched grid is also employed in the normal direction, with clustering near the blade 
surface to capture the large velocity gradients that occur, and to ensure that, for turbulent 
flow, there are enough points in the near- wall region to adequately resolve the laminar sublayer. 
Geometric stretching is again employed, with 


A? 7 n+1 = K„ A Vn , n = 1, 2, ... ,N (5.47) 

where n is a normal mesh point index (n = 1 at the blade surface or reference wake streamline 
and increases with increasing distance into the flow). Since the resolution requirements are 
a function of the solution, the normal grid stretching, K v , and the spacing at the wall, A rft, 
must be determined for each case. For turbulent flow, this is accomplished by monitoring the 
value of Y + at the first mesh point off the wall (i.e., l£ 2 ), where Y + is a Reynolds number 
based on the friction velocity, (f w /p w ) 1 ^ and the normal distance from the blade surface. For 
accuracy, T 2 + should be close to one, which can be achieved by adjusting Arft. The stretching 
parameter K v is then adjusted to place the outer edge of the mesh, tj - t] N , far enough from 
the surface to allow the edge conditions to be approached asymptotically. The number of 
points used in the normal direction is chosen to allow the flow over the entire viscous layer to 
be resolved accurately at all streamwise stations — generally, between 50 and 100 points are 
sufficient. 

The constant value of the time-step used for the temporal discretization is chosen based 
on the frequency of the imposed unsteady excitation and the number of time-steps prescribed 
for each period of oscillation, generally between 20 and 50. More complete details on the 
finite-difference analysis, used herein, for viscous-layer flows are given in [VBHA91]. 
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5.4 Numerical Results: Unsteady Flows Driven by Acoustic Disturbances 

Response predictions will be given below for two-dimensional blade rows operating at 
subsonic Mach numbers to demonstrate the present capabilities of the unsteady viscous-layer 
analysis, UNSVIS. In particular, this analysis will be applied to predict the viscous-layer 
responses of flat-plate and turbine cascades to acoustic excitations. For these applications, 
the blades are assumed to be stationary in the cascade frame of reference. Hence, R, = 0, 
U = V and Hr = H T . 


Unstaggered Flat- Plate Cascade 

As a benchmark case to demonstrate the UNSVIS analysis, we consider unsteady flows 
through an unstaggered flat-plate cascade. Here, the uniform inlet velocity V_oo is aligned 
with the mean positions of the blade chord lines. Hence, the inviscid steady background flow 
is uniform with velocity, pressure and total enthalpy given by V = e^, P = (7 M 2 ) -1 , and 
Ht = 1/2 + (7 — 1 ) -1 M -2 , respectively. The unsteadiness is excited by incident pressure 
disturbances, that carry energy toward the blade row from either far upstream (—00) or far 
downstream (+00). The acoustic excitations travel in the axial (or chordwise) direction, i.e., 
k v = K y = a jG = 0, and are therefore, described by a velocity potential of the form [cf. (3.20) 
and (3.21)]. 

<£( x , t) = —Re{ [*(/<*, TO o + w)] -1 p/,Too expf^/c^oof + wf)]} , (5.48) 

where u>, = k x and p/, TO o are the temporal frequency, axial wave number and complex 

amplitude of the unsteady pressure excitation, respectively. The axial wave number is related 
to the temporal frequency by [cf. (3.24)] 


i ± M ' (5.49) 

Note that for |/c^ i+00 | > |k^,_oo|, because a downstream acoustic excitation travels upstream 
against the freestream flow direction, whereas an excitation from upstream travels in the 
direction of the freestream flow. 

The flat-plate cascade produces no response to the pressure excitation described by (5.48); 
therefore, this equation describes the first-harmonic component of the inviscid unsteady flow. 
The time-dependent velocity, V, pressure, P, and total enthalpy, Hj , in this flow are given 
by 

V = V + v + ... = [1 ± MRe{p I ^ <x> exp[i(K tToo ( + ut)]}}e i + ... , 

P = P + p+... = (7-^ 2 ) _1 + Re{p/, Too exp[t(K 4>Too ^+a;<)]}-f-... , (5.50) 

and 

Ht = H t + hr + • •• = Ht + (1 ± M)Re{p/ >Too exp[i(/c iiToo ^ + u>f)]} + ... , 

where v = V<^,p = — D4>jDt , and Jit = — D4>/Dt-\- V^-e^ are the first-harmonic components 
of the unsteady velocity, pressure and total enthalpy, respectively. The real and imaginary 
components of the complex amplitude of the first-harmonic unsteady pressure acting along 
each blade and wake surface (0 < x - £ < 2) for a flat-plate cascade operating at M = 0.5 
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Figure 5.1: Unsteady pressure distributions along a blade and its wake for an unstaggered 
flat-plate cascade (0 = 0 = 0 deg, = 0.5 and G = 1) subjected to a unit-amplitude, 

|p/,-oo| = 1, pressure excitation from upstream with k,,_oo = crG~ x = 0. 



and subjected to upstream acoustic disturbances with p/,_oo = (1,0) and u> = 1, 2, 5 and 10 
are shown in Figure 5.1; similar results for downstream acoustic disturbances, are shown in 
Figure 5.2. 

The inviscid results (5.50) have been used to provide the pressure and edge-condition 
information for a series of unsteady viscous-layer calculations. The latter will allow us to 
examine the effects of the amplitude, temporal frequency and axial wave number of a pressure 
excitation on the behavior of various viscous-layer response quantities; in particular, the dis- 
placement thickness, I(x,<), surface shear stress f w = (Re)~ l (fidV T /dn) w and wake centerline 
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Figure 5.2: Unsteady pressure distributions along a blade and its wake for an unstaggered 
flat-plate cascade (with fi = 0 = Odeg, M_ <*, = 0.5 and G = 1) subjected to a unit-amplitude, 
|p/ 1+ oo | = 1, pressure excitation from downstream with K ,, t+00 = aG~ l = 0. 


velocity Vq. We can decompose each of these quantities into Fourier series; e.g., 


= / [1 -pVr/(p e V T ,e)]dn= 6 n (x)exp(inut) , 

J 0 

n=— oo 


(5.51) 


where 


w r2ir/w ~ 

M x ) = 2 ~ J Q 6(x,t)exp(-inurt)dt , n = 0,±1,±2, ... , 


(5.52) 


to examine the behavior of their Fourier components. In this way we can gain insight into 
the relative importance of nonlinear effects on the flow in a viscous layer and, therefore, into 
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whether or not a linearized viscous analysis could be applied to provide efficient and meaningful 
unsteady viscous-layer solutions. Note that the lower limit on the integral in (5.51) must be 
changed to — oo, if the wake displacement thickness is to be determined. 

A series of results were determined for acoustic excitations travelling towards the blade row 
from upstream or downstream. In each case the steady Mach number is 0.5 and the Reynolds 
number, Re, is 10 6 . The excitations occur at the frequencies, w, listed above, and at complex 
amplitudes, pj, ±00 , of (0.1,0), (0.3,0), (0.5,0) and (0.75,0). The pressure, P = in the 

steady background flow is 2.857. The viscous-layer numerical calculations were initiated at 
(x,t) = (0.01,0). Laminar similarity solutions were used for 0 < x < 0.1 and imposed as the 
initial conditions in x and t for the numerical calculations. The viscous flows were assumed 
to undergo an instantaneous transition from laminar to turbulent flow at x = 0.02. 

The boundary layer calculations were carried out using 51 uniformly stretched (with K v = 
1.10) points across each boundary layer, with Arj = 0.0175 at the blade surface. In the wake 
the grid consisted of 101 points across the viscous layer stretched in the same manner as for the 
surface boundary layers. A total of 25 uniform time steps were used per temporal period of the 
unsteady excitation. Two different axial or streamwise mesh distributions were used — one 
for the excitations coming from upstream; the other, for excitations coming from downstream. 
In each case the streamwise distribution was selected so that there were at least 20 axial mesh 
lines per wave length, 2i t/k^oo, for the highest frequency considered, i.e., u = 10. 

For disturbances originating upstream, a variably spaced streamwise mesh was used with 
points clustered near the blade leading and trailing edges. The minimum streamwise spacing 
on the blade was Ax « 0.0177 at the blade edges, and the maximum was Ax « 0.0611 
near midchord. The stretching used in the wake was identical to that used for the forward 
portion of the blade, with Ax « 0.0177 in the first wake interval and monotonically increasing 
to Ax rs 0.11 at the downstream boundary of the computational domain (x = 2.0). The 
resulting grid had 29 points on the blade surface and 20 points along the wake. For the 
disturbances originating downstream a nearly uniform grid was employed, with Ax « 0.025 
on the blade and 0.025 < Ax < 0.030 along the wake, where the grid was mildly stretched 
in the flow direction to distribute the points throughout the interval x € (1,2]. The resulting 
grid had 40 points along the blade and 36 points along the wake. 

First, we consider the solution for a pressure disturbance from upstream with pj oo = 
(0.5, 0), u> = 5 and a = 0 deg. The temporal means, <$ 0 , r wfi and Vg >0 , and the magnitudes, |<5„|, 

|r w , n | and |% n |, of the first two harmonics (n = 1, 2) of the displacement thickness, 6, surface 
shear stress, t w , and wake centerline velocity, Vg, as determined by the unsteady viscous-layer 
solution, are presented in Figure 5.3 along with the corresponding steady (|p/,-oo| = 0) results, 
which were also determined using the UN S VIS code. The steady displacement thickness and 
surface shear stress are given by 

/•oo 

S(x) = / [1 - pV T /(p e V T , e )}dn (5.53) 

and t w = (R ,e)~ l (jidV T /dn) w , where p and V T are the density and streamwise velocity, respec- 
tively, in the steady background flow. The differences between the steady and the temporal 
mean values of the unsteady viscous quantities, and the amplitudes of the higher (n > 2) har- 
monic unsteady quantities provide a measure of the relative importance of nonlinear effects 
on the unsteady flow in the viscous layer. 
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Figure 5.3: Temporal mean and Fourier magnitudes of the displacement thickness, I, wall 
shear stress, t w , and wake centerline velocity, V<g, for turbulent flow through an unstaggered 
flat-plate cascade (f 1 = 0 = Odeg, G = 1, M = 0.5 and Re = 10 6 ) subjected to an incident 
pressure disturbance from upstream with pi-oo = (0.5,0), u> = 5 and a = Odeg. 

The results in Figure 5.3 indicate that nonlinear effects axe relatively small for the un- 
steady flow driven by the prescribed upstream pressure excitation. However, similar results 
in Figure 5.4 for an acoustic excitation from downstream, with p/ +OG = (0.5,0), u> = 5 and 
<r = 0, indicate the nonlinear content in the viscous-layer displacement-thickness response 
to be quite significant. An unexpected result of the latter calculation is the increase in the 
predicted time-mean of the unsteady displacement thickness with increasing distance along 
the wake (i.e., as x — > 2.0). In an attempt to determine whether this effect is physical or nu- 
merical in origin, an unsteady viscous solution was calculated using a grid with twice as many 
uniformly distributed streamwise points. It was found that, although slightly less pronounced, 
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Figure 5.4: Temporal mean and Fourier magnitudes of the displacement thickness, <5, wall 
shear stress, t w , and wake centerline velocity, Vg, for turbulent flow through an unstaggered 
flat-plate cascade (SI = 0 = Odeg, G = 1, M = 0.5 and Re = 10 6 ) subjected to an incident 
pressure disturbance from downstream with p/, +o0 = (0.5,0), u = 5 and a — Odeg. 

a similar behavior was also present in the fine grid calculation. This unexpected behavior of 
the viscous layer solution for an upstream traveling acoustic disturbance will require some 
further study to resolve. 

The behaviors of the zeroth- and first-harmonic components of the viscous parameters 
are illustrated in Figures 5.5 and 5.6 for upstream acoustic excitations with pi-oo = (0.5,0), 
a — o, and oj = 1,2,5 and 10. Corresponding results for downstream acoustic excitations 
are given in Figures 5.7 and 5.8. For the upstream excitations, the time-mean values and the 
magnitudes of the first-harmonic viscous-layer response quantities (Figure 5.5) are not affected 
significantly by changes in excitation frequency. In addition, the first-harmonic displacement 
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Figure 5.5: Fourier amplitudes of 6, t w and Vq for an unstaggered flat-plate cascade (ft = 
0 = Odeg, G = 1, M = 0.5 and Re = 10 6 ) subjected to pressure excitations from upstream 
with p/,_oo = (0.5, 0) and a = 0 deg. 

thickness is essentially out-of-phase with the first-order streamwise velocity fluctuation, v r<e , 
at the edge of the viscous layer (Figure 5.6), and the first-harmonic wall shear stress and wake 
centerline velocity are essentially in phase with this edge-velocity fluctuation. 

The corresponding results in Figures 5.7 and 5.8 for acoustic excitations from downstream 
show that frequency has a much greater impact on the temporal mean and the first-harmonic 
amplitudes of the viscous-layer response quantities. In particular, the time-mean displacement 
thicknesses vary significantly with u>, indicating that nonlinear effects become increasingly 
important with increasing frequency. As shown in Figure 5.8, the first-harmonic displacement 
thickness lags the fluctuating component of the streamwise velocity at the edge of the viscous- 
layer by a phase angle ranging from approximately —140 to —180 deg; the first-harmonic 


108 





X 


Figure 5.6: Relative phase angles of first-harmonic displacement thickness, wall shear stress, 
and wake centerline velocity for an unstaggered flat-plate cascade (fl = 0 = Odeg, G = 1, 
M = 0.5 and Re = 10 6 ) subjected to pressure excitations from upstream with p/,-oo = (0-5, 0) 
and a = Odeg. 

wall shear stress leads this fluctuating velocity by an angle ranging from 0 to 40 deg; and 
the minimum wake-streamwise velocity is approximately in-phase with the wake-edge velocity 
fluctuation. 

As a final illustration, we consider the response of the viscous layer to acoustic excitations 
at different amplitudes. Numerical results are presented in Figures 5.9 and 5.10 for excitations 
at uj = 5, a = 0 and p/, TOO = (0.1,0), (0.3,0), (0.5,0) and (0.75,0). For disturbances from 

upstream (Figure 5.9) the effects of unsteadiness on the time-mean values of 8 and f w are 
small, except, perhaps, for the highest disturbance amplitude, |pj,_oo| = 0.75, considered. The 
effects of disturbance amplitude on the time-mean wake centerline velocity are negligible. The 


109 





10.0 


Re 1 / 2 ! 

5.0 


0.0 




Figure 5.7: Fourier amplitudes of 6 v and Vg for an unstaggered flat-plate cascade (ft = 
0 = Odeg, G = 1, M = 0.5 and Re = i0 6 ) subjected to pressure excitations from downstream 
with p/,+oo = (0.5, 0) and a = Odeg. 


results for the acoustic excitations from downstream in Figure 5.10 once ag indicate that 
such disturbances produce much stronger nonlinear viscous-layer response^ articularly in 
the displacement thickness. For both the upstream and downstream acoustic excitations, the 
phase angles (not shown) of the first- harmonic viscous quantities relative to the edge velocity 
fluctuation are almost independent of the excitation amplitude. 

The foregoing results indicate that the viscous-layer response parameters 6, f w and Vg be- 
have essentially in a linear manner for acoustic excitations originating upstream of the blade 
row. For excitations coming from downstream, nonlinear effects can be significant for high 
temporal frequencies (say u> > 5) and/or high excitation amplitudes (|p/ >+00 | > 0.5). For 
M = 0.5, the magnitude of the axial wave number for a pressure excitation from downstream 
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Figure 5.8: Relative phase angles of first-harmonic displacement thickness, wall shear stress, 
and wake centerline velocity for an unstaggered flat-plate cascade (fl = 0 = Odeg, G = 1, 
M = 0.5 and Re = 10 6 ) subjected to pressure excitations from downstream with p/,+oo = 
(0.5,0) and a = Odeg. 

is three times that for an excitation from upstream. This produces a corresponding ratio in 
the magnitudes of the pressure gradients that drive the viscous-layer solution, and is perhaps 
responsible for the relatively strong nonlinear content in the responses of the viscous layer to 
downstream disturbances. Another factor is the direction of propagation of the acoustic exci- 
tation. We have performed a series of calculations to isolate these effects and have found that 
the latter, i.e., direction of propagation, to be the dominant factor in promoting a nonlinear 
response of the viscous layer. It should also be noted that Navier-Stokes calculations have 
been performed recently for compressor cascades excited by acoustic disturbances [DV94]. 
The Navier-Stokes results indicate that the unsteady surface pressure responses to acoustic 
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Figure 5.9: Fourier amplitudes of 8, t w and Vg vs amplitude, Ip/.-ool, for an acoustic excitation 
from upstream, with u = 5 and a = Odeg, interacting with an unstaggered flat-plate cascade 
(SI = 0 = Odeg, G = 1, M = 0.5 and Re = 10 6 ). 

excitations from upstream are linear over a wide range of excitation amplitudes, but acoustic 
excitations from downstream produce responses with significant nonlinear content. 

The spurious predictions for the behavior of the displacement thickness in the far wake, 
i.e., increasing with increasing x , for high frequency or high amplitude disturbances from 
downstream is a somewhat disturbing aspect of the present unsteady viscous-layer solutions. 
A series of laminar calculations were performed, and the results showed a similar behavior. 
Thus, the increase in displacement thickness in the far wake ( x — ► 2) cannot be attributed to 
the turbulence model used in the present study. 
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Figure 5.10: Fourier amplitudes of S, f w and Vg vs amplitude, |p/ >+00 |, for an acoustic excita- 
tion from downstream with, ui = 5 and o — Odeg, interacting with an unstaggered flat-plate 
cascade (f 1 = 0 = Odeg, G = 1, M = 0.5 and Re = 10 6 ). 


Turbine Cascade 

We proceed to consider an unsteady flow through a “real” blade cascade to demonstrate the 
coupled LINFLO/UNSVIS weak, inviscid/viscid interaction, solution capability. In particular, 
we consider a flow at a Reynolds number, Re, of 5.0 x 10 5 , an inlet Mach number, M_ <», of 
0.19 and an inlet flow angle, tl-oo, of 45 deg, through the turbine cascade known as the 4th 
Standard Configuration. The unsteadiness is excited by an acoustic or pressure disturbance 
from upstream, characterized by pi-oo = (0.35,0), /c^.-oo = &/G = 0 and u> — 1.0. The 
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Figure 5.11: Contours of the in-phase component (real part) of the unsteady pressure for the 
turbine cascade subjected to an acoustic excitation from upstream with pi-oo = (0.35,0), 
u = 1 and <j = 0. 


inviscid steady Mach number field and steady Mach number distribution along a blade surface 
are shown in Figure 3.2. For the unsteady flow being considered, the unsteady pressure at 
a given point in the field is essentially in phase with the acoustic excitation. The in-phase 
component (real part) of the unsteady pressure field, as determined using LINFLO, is depicted 
in Figure 5.11. The acoustic excitation at w = 1, <7 = 0 produces a superresonant (1,1) 
acoustic response. Therefore, the pressure contours in Figure 5.11 indicate much different 
pressure amplitudes at the far upstream and downstream boundaries than that associated 
with the acoustic excitation alone. 

The temporal mean value and the upper and lower bounds (envelope) for the inviscid 
surface velocity are shown in Figure 5.12. The amplitude of the first-harmonic component of 
the surface velocity, which is proportional to |p/,_oo|, is relatively small for this case. But, for 
pressure excitations of higher amplitude, the viscous-layer separates from the suction surface 
of a blade, just upstream of the trailing edge, thereby precluding the continuation of the direct 
unsteady viscous-layer calculation into the wake. 
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Figure 5.12: Temporal mean values and upper and lower bounds of the inviscid surface 
(viscous-layer edge) velocity for the turbine cascade subjected to an upstream pressure exci- 
tation with pi - oo = (0.35,0), u> = 1 and a = 0. 

Since the turbine blade has a blunt leading edge, the unsteady stagnation-point analysis, 
described in § 5.3, was applied at each time step to provide the upstream velocity profiles 
needed to advance the viscous-layer solutions along the blade suction and pressure surfaces. 
For this case, the motion of the stagnation point is confined to a small interval of length 
2|r/| ~ 0.0037. A series of velocity profiles associated with the periodic flow within the 
stagnation region at four different times, i.e., t — 7r/2, 7r, 3 tt/ 2 and 2tt , are shown in Figures 5.13 
and 5.14. Here the instantaneous streamwise- velocity profiles are presented in the blade fixed 
frame of reference, and the abscissa on each plot refers to the location on the blade surface 
at which the velocity profile is determined. The interval shown in Figure 5.13 is centered 
about the mean location of the stagnation point; that in Figure 5.14, about the instantaneous 
location, i.e., the location at which U T>e = 0 in the blade-fixed frame of reference. Note 
that the velocity profiles are shown over a much narrower interval, [— 0.02|r/|, 0.02 |t/|], in 
Figure 5.14, where the profiles indicate that reverse flow occurs in the immediate vicinity of 
the instantaneous stagnation point location. 

The viscous-layer calculation for the reference turbine blade and its wake was performed 
assuming that instantaneous transition from laminar to turbulent flow occurs at r /tje = 0.05 
on both the pressure and suction surfaces of each blade. Here, r is the distance along the blade 
surface measured from the mean, leading-edge, stagnation point location, and the subscript 
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Figure 5.13: Velocity profiles in the neighborhood of a mean stagnation point location for the 
turbine cascade subjected to an upstream acoustic excitation with p/,_oo = (0.35, 0), u = 1 
and cr = 0: (a) t = 7r/2; (b) t = tt; (c) t = 37r/2; (d) t = 2 tt. 



Figure 5.14: Velocity profiles in the neighborhood of an instantaneous stagnation point lo- 
cation for the turbine cascade subjected to an upstream acoustic excitation with pi~ 0 o = 
(0.35,0), u> = 1 and a = 0: (a) through (d) as in Fig. 5.13. 
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TE refers to the values of r at the blade trailing edge. The grid used in this calculation 
had 77 points along the blade surface and 54 points along the reference wake line. It was 
stretched with At rj 0.002 at the furthest upstream point and At as 0.0001 at the trailing 
edge. The largest value of At on the blade, i.e., 0.052, occurs near midchord. The streamwise 
intervals grow aft of the trailing edge from At 0.0001 to approximately 0.083 one chord 
length downstream of the blade row. The viscous-layer calculation was carried out using 71 
uniformly stretched points across each boundary layer, with K n = 1.045 and A 7? = 0.04 at 
the blade surface. The wake grid consisted of 141 points across the viscous layer stretched 
in the same manner as on the blade surface. A total of 40 uniform time steps were used per 
temporal period of the unsteady excitation. 

Results of the unsteady viscous-layer calculation are shown in Figures 5.15 through 5.17. 
Temporal mean values and upper and lower bounds for the displacement thickness and wall 
shear stress along the upper and lower surfaces of the reference (m = 0) turbine blade are 
shown in Figure 5.15; corresponding results for the wake displacement thickness, 6w, and 
minimum streamwise velocity, V T , win in Figure 5.16. Here, the upper and lower bounds of a 
viscous-layer response quantity, say the displacement thickness, are defined by 

~± 00 

5 =<o±2£|«„|. (5.54) 

n= 1 

Wake velocity profiles at four different instants of time are depicted in Figure 5.17. The 
unsteady response of the viscous layer is essentially linear for this example, i.e., the temporal 
mean and the steady viscous solutions are almost identical, and the Fourier amplitudes of the 
higher (n > 2) harmonic components of 6, f w and K min are negligible. 

The foregoing results demonstrate the new capabilities that have been added to the 
UNSVIS code under Contract NAS3-25425, i.e., an unsteady stagnation region analysis and 
an unsteady wake analysis. These results also demonstrate the present weak, inviscid/ viscid 
interaction, prediction capability that results from a sequential coupling of a linearized, in- 
viscid, unsteady solution, determined using LINFLO, and a nonlinear unsteady viscous-layer 
solution. Unfortunately, since boundary-layer separation usually occurs in realistic configu- 
rations and a weak interaction analysis breaks down in such cases, the present analysis has 
only a limited range of application. This is particularly true for compressor cascades where, 
because of adverse mean pressure gradients, local separations almost always occur. There 
is, therefore, an important need to develop a simultaneous coupling (or strong interaction) 
solution procedure for unsteady cascade flows. 
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Figure 5.15: Temporal mean values and upper and lower bounds for the displacement thickness 
and surface shear stress along a turbine blade surface for an unsteady flow excited by an 
upstream pressure excitation with p/_<x> = (0.35,0), u> = 1 and a = 0. 
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Figure 5.16: Temporal mean values and upper and lower bounds for the displacement thickness 
and minimum wake streamwise velocity along a turbine wake for an unsteady flow excited by 
an upstream pressure excitation with pi-oo = (0.35, 0), u> = 1 and a — 0. 





6. Concluding Remarks 


Under the present effort, we have contributed to the development of efficient, and rea- 
sonably comprehensive, unsteady fluid dynamic analyses that can be used in turbomachinery 
aeroelastic and aeroacoustic design studies. In particular, we have invoked the assumptions 
of high-Reynolds-number, “attached” flow and small-amplitude unsteady excitation to de- 
velop asymptotic unsteady aerodynamic models that apply to realistic cascade configurations 
and mean operating conditions. Based on the high Reynolds number assumption, an invis- 
cid/viscid interaction model has been formulated for unsteady flows through two-dimensional 
cascades. In addition, based on the small- disturbance assumption, a linearized analysis has 
been formulated for the unsteady flow in the inviscid region, but at this point, the unsteady 
flows in viscous layers, i.e., boundary layers and wakes, are still regarded as nonlinear. To 
further expedite the unsteady flow predictions, we have assumed that the nonlinear steady 
background flow in the inviscid region is isentropic and irrotational — an assumption that 
leads to considerable simplifications in the equations that describe the behavior of inviscid 
unsteady perturbations. The foregoing assumptions can lead to very efficient predictions of 
the unsteady pressure responses of realistic cascades to prescribed structural (blade) motions 
and external unsteady aerodynamic (entropic, vortical, and acoustic) disturbances. 

To provide a strong inviscid/ viscid interaction analysis for unsteady cascade flows, several 
component analyses must be constructed, along with methods for coupling these components 
into an overall solution procedure. Here, the component analyses include a full potential 
analysis to predict the isentropic and irrotational steady background flow, a linearized inviscid 
analysis to predict the behaviors of unsteady entropic, vortical and acoustic perturbations, 
and a nonlinear unsteady viscous layer analyses, which, in the future, might be replaced by 
nonlinear steady and linearized unsteady analyses, to predict the flows in boundary layers 
and wakes. In the present report, we have described and demonstrated the linearized inviscid 
analysis LINFLO, the steady, strong, inviscid/viscid interaction analysis SFLOW-IVI and 
the nonlinear unsteady viscous-layer analysis UNSVIS. The steady full-potential analysis, 
SFLOW, developed at NASA Lewis under a related research program [HV93, HV94], serves 
as the inviscid component of the steady inviscid/viscid interaction analysis, and is also used 
to provide the steady background flow information needed for a LINFLO linearized inviscid 
analysis. The SFLOW-IVI analysis entails the iterative coupling of SFLOW to a steady, 
inverse, viscous-layer analysis. 

The LINFLO analysis describes unsteady perturbations of a potential steady background 
flow. It applies to unsteady flows excited by prescribed blade vibrations and external entropic, 
vortical and acoustic disturbances. Applications of LINFLO to blade flutter and blade-row 
aeroacoustic response predictions are described in [Ver89a, UV91, Ver93] and [KV93b, KV94], 
respectively. Under the present contract, LINFLO was extended to predict the unsteady flows 
excited by entropic and vortical gusts. Because the steady background flow is assumed to be 
potential, closed form solutions can be determined for the entropy and vorticity or rotational 
velocity fluctuations. Consequently, numerical field methods are only required to determine 
the unsteady potential, and hence, the unsteady pressure. The former is governed by an 
inhomogeneous wave equation in which the source term depends upon the rotational velocity 
field. Since only a single partial differential equation must be solved numerically, the LINFLO 
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analysis provides very efficient unsteady aerodynamic response predictions. Moreover, such 
predictions have been shown to be in good agreement with those based upon time-accurate, 
nonlinear, Euler and Navier-Stokes solutions for subsonic unsteady flows excited by vortical 
and acoustic excitations [DV94] and for subsonic and transonic unsteady flows excited by 
blade vibrations [AV94]. 

In future work, the LINFLO vortical gust response prediction capability should be ex- 
tended to transonic flows. Also, the convected potential developed in [AG89] and used herein, 
cf. (3.33), should be modified to remove the indeterminacy in the unsteady velocity, associ- 
ated with entropic gusts, at blade and wake surfaces. This would lead to improved resolutions 
of unsteady flows excited by entropic disturbances. Finally, since we have experienced some 
difficulties in providing accurate unsteady pressure responses for thick, highly-cambered, tur- 
bine blades, particularly the pressure responses to high frequency vortical gusts, work should 
be directed towards improving LINFLO for application to vortical gust/turbine blade- row 
interactions. 

The SFLOW-IVI analysis describes high Reynolds number, steady flows, containing re- 
gions of strong inviscid interaction, including local separations, through turbomachine cas- 
cades. Here, the flow in the outer or inviscid region is governed by the full potential equation; 
that in the inner or viscous layer regions, by Prandtl’s viscous layer equations. Inviscid and 
viscid solutions are iteratively matched using a semi-inverse global iteration procedure until a 
converged result for the complete flow field is determined. The SFLOW analysis [HV93, HV94], 
with modified surface conditions to account for the effects of viscous displacement, is used 
to determine the inviscid component of the flow, and an inverse viscous layer analysis, to 
determine the flows in viscous-layer regions. 

As part of the present effort, the SFLOW-IVI analysis has been applied to predict steady 
viscous flows through compressor and turbine cascades, including flows with extensive trailing- 
edge separations. Results, for selected cases, particularly those for surface shear stress, have 
been shown to be in good agreement with corresponding Navier-Stokes predictions. The 
inviscid components of SFLOW-IVI solutions have also been used to provide the steady 
background flow information for a linearized inviscid analysis, in an effort to determine the 
effects of steady viscous displacement on blade f Iter margins. For the limited applications 
considered here, steady viscous displacement effects tend to be stabilizing for torsional blade 
motions, but destabilizing for bending vibrations. 

As it stands, the SFLOW-IVI analysis is efficient and robust, but it is restricted to sub- 
sonic, adiabatic flows at unit Prandtl number. In future work the viscous-component of the 
SFLOW-IVI analysis could be easily extended so that heat conducting flows at arbitrary 
Prandtl numbers can be considered. More importantly, the SFLOW-IVI analysis should be 
extended for application to transonic flows. A steady, transonic, inviscid/ viscid interaction 
analysis would be of value for steady-state design applications, and could be used in conjunc- 
tion with LINFLO to determine the impact of steady viscous displacements on unsteady shock 
loads. Recent Navier-Stokes studies [AV94] have indicated that, for cascades operating near 
design, the effects of viscous- displacement on unsteady aerodynamic response are significant 
only in the vicinity of a shock, where such effects tend to weaken the shock and shift its mean 
position. In this case, a steady analysis that accounted for such effects could be coupled to 
LINFLO to provide useful unsteady aerodynamic response information for transonic flows. 
This type of approach has been applied successfully in external aerodynamics [Edw93]. Its 
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use for internal flows could lead to very efficient unsteady aerodynamic response predictions 
for viscous transonic unsteady flows. 

The nonlinear unsteady viscous-layer analysis and code, UNSVIS, has been developed to 
predict the unsteady flows in blade boundary layers and wakes. This analysis can be applied 
to predict the viscous-layer responses that arise from imposed inviscid conditions at the blade 
and wake surfaces. At present only a direct viscous-layer calculation and a weak or sequential 
coupling of inviscid and viscous-layer solutions has been considered. The development of a 
strong inviscid/viscid interaction analysis involving a simultaneous coupling of the inviscid and 
viscous solutions should be considered in future work. Under the present effort, the existing 
unsteady viscous layer analysis has been extended, by incorporating a similarity analysis for 
the flow in the vicinity of a moving stagnation point and by accounting for the jumps in the 
inviscid flow variables across vortex-sheet unsteady wakes, and applied to predict unsteady 
cascade flows excited by prescribed acoustic excitations. 

The coupled linearized inviscid LINFLO and nonlinear viscous-layer UNSVIS analyses 
have been demonstrated via applications to an unstaggered flat-plate cascade subjected to 
acoustic excitations from upstream and downstream, and to a turbine cascade interacting with 
an upstream acoustic excitation. The flat-plate example is, perhaps, the simplest unsteady 
cascade problem that can be analyzed, both because of its geometric simplicity and because 
the unsteady pressure is nonsingular at the flat-plate leading edges. The numerical results 
indicate that the viscous layer responds linearly, for the most part, to acoustic excitations from 
upstream, but significant nonlinear response components occur for downstream excitations 
at high temporal frequency and/or high amplitude, which travel against the mainstream 
flow velocity. A similar conclusion, based on Navier-Stokes solutions for acoustic excitations 
interacting with a low-speed compressor cascade has also been reported in [DV94]. The 
numerical results for the turbine demonstrate the present weak inviscid/viscid interaction 
solution capability for a realistic cascade configuration. 

Because of boundary-layer separation, the range of application of a weak unsteady in- 
viscid/viscid interaction analysis is severely limited. For example, the mean pressure rise 
produced by a compressor blade row typically causes boundary-layer separations near blade 
leading edges, thereby precluding a continuation of a direct viscous-layer calculation along 
the blades and into their wakes. Thus, the development of a strong inviscid/viscid interaction 
analysis for unsteady flows will be needed so that the effects of steady and unsteady viscous 
displacements on the unsteady aerodynamic responses of blade rows can be predicted. The 
responses of viscous layers to prescribed blade vibrations and entropic and vortical excita- 
tions are other issues that require further study. In particular, as an important next step, 
the UNSVIS analysis should be applied to unsteady cascade flows excited by prescribed blade 
motions. The viscous layer responses in such flows should be studied to understand the rela- 
tive importance of nonlinearities and the potential impact of unsteady viscous displacement 
on unsteady pressure response. 
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